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The  Jovian  decametric  radiation  is  known  to  be  excited, 
in  part,  by  the  satellite  lo.  However,  much  of  the  radiation 
is  apparently  not  related  to  lo  and  the  provocation  for  this 
component  is  not  understood.  This  study  examines  the  non-Io 
component  by  the  systematic  removal  of  all  information 
associated  with  the  "favored"  orbital  positions  of  lo  at  18, 
22,  26,  and  27  MHz. 

Three  specific  examinations  have  been  made.   The  first 
has  considered  the  possibility  of  a  correlation  between  the 
activity  of  the  radio  sources  and  density  waves  in  the  solar 
wind.  A  byproduct  of  this  aspect  of  the  study  is  the  develop- 
ment of  two  models  for  solar  wind  velocity  stream  inter- 
actions.  These  models  predict  the  development  of  large  density 
waves,  mainly  past  1  AU,  in  the  solar  wind.   The  waves  are  the 
result  of  interactions  (collisions)  between  high  speed  and  low 


speed  particle  streams  in  the  solar  wind. 

A  positive  correlation  between  these  proposed  waves  and 
the  source  A  and  source  C  emission  is  found  in  the  two  study 
periods,  1967-68  and  1974.  During  the  same  periods  there  is 
evidence  of  a  negative  correlation  with  source  B. 

The  second  part  of  the  study  has  examined  the  activity 
occurrence  spectra  for  evidence  of  periodicity.  A  27  day 
period  was  found  in  the  1974  data  which  may  be  associated 
with  the  solar  correlation.   In  the  1967-68  and  1974  appari- 
tions a  prominent  peak  was  found  at  6.35  days  in  both  source 
A  and  C.   This  peak  is  coincident  with  an  alias  of  a  principal 
feature  in  the  spectrum  of  a  joint  lo-Europa  exciter.  No  other 
satellite  effect  can  cause  this  feature. 

Another  spectral  peak  was  found  at  2.7  days  in  both 
apparitions  in  source  A.   This  maximum  is  not  explainable  at 
this  time. 

The  third  examination  has  considered  the  question  of 
randomness.   The  object  was  to  determine  whether  it  was 
possible  that  the  non-Io  component  was  the  result  of  some, 
possibly  intrinsic,  random  process.   In  this  study  two  char- 
acteristics have  been  examined.   The  first  was  the  length 
of  the  storms  themselves.   Employing  the  frequency  of  length 
distribution  of  the  storms  it  is  shown  that  storm  lengths  are 
not  randomly  distributed. 

The  other  topic  was  the  question  of  the  likelihood  that 


the  6.4  day  period  features  seen  in  the  Chree  analysis  of  the 
radio  activity  could  be  the  result  of  random  influences.   This 
effect  was  shown  to  be  non-random  with  a  confidence  exceeding 
99.5% 

The  study  concludes  that  below  the  strata  of  events  known 
to  be  associated  with  lo,  there  are  secondary  and  tertiary 
exciters  which  influence  the  non-random  emission  in  the  18  to 
27  MHz  regions. 


CHAPTER  1 
INTRODUCTION 

Jovian  Decametric  Radiation 

Sources 

The  process  by  which  the  intermittent  "storms"  of 
Jovian  decametric  emissions  are  generated  has  been  the  subject 
of  considerable  study  since  the  radiation  was  discovered  by 
Burke  and  Franklin  (1955). 

An  examination  of  pre-discovery  data  by  Shain  (1956) 
revealed  that  the  radiation  was  not  uniformly  distributed 
with  respect  to  the  Jovian  central  meridian  longitude  (CML) . 
He  found,  rather,  that  the  signals  originated  from  specific 
"sources"  which  were  apparently  fixed  to  Jupiter's  high 
latitude  System  II  longitude  scheme.   Thus,  these  source  loca- 
tions rotated  with  the  planet,  though  further  studies  by  many 
observers  showed  a  slight  drift  with  respect  to  System  II. 
The  adoption  of  a  radio  longitude  scheme.  System  III,  largely 
removed  this  small  discrepancy. 

If  one  plots  the  probability  of  activity  in  the  vicinity 
of  20  MHz  versus  CML,  the  result  is  a  histogram  such  as  Figure 
1.  The  features  labeled  A,  B,  and  C  represent  the  three  major 
source  regions.  The  nomenclature,  source  "A",  "B",  or  "C"  was 
adopted  by  the  Florida  group  (Carr  et  al. ,  1961)  early  in  their 
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3 
studies  arid  will  be  used  in  this  presentation.   However, 
sources  A,  B,  and  C  correspond  respectively  to  sources  2, 
1,  and  3  defined  by  Douglas  (1964)  and  the  main,  early,  and 
third  sources  of  Dulk  ( 1965a, b) . 

lo  Effect 

Subsequently,  Bigg  (1964)  determined  that  there  was  a 
component  of  the  decametric  emission  which  was  strongly  corre- 
lated with  the  phase  from  superior  geocentric  conjunction  of 
the  inner  most  Galilean  satellite,  lo.  More  recent  studies 
(e.g.,  Register,  1968;  Miller  et  al.,  1972)  have  shown  that 
all  three  of  the  classical  decametric  sources  contain  both  an 
lo  related  component  and  non-Io  related  component.   Source  A, 
in  particular,  and  to  a  lesser  extent,  source  C,  contain 
significant  non-Io  components  in  the  vicinity  of  20  MHz.   This 
situation,  depicted  by  Leacock  (1971)  in  Figure  2,  seems  to 
suggest  that  the  position  of  lo  enhances  a  pre-existing  radia- 
tion process.   More  recently,  Desch  (1976)  has  shown  that  source 
B  also  has  a  significant  non-Io  component  at  lower  flux  levels. 

The  data  presented  in  Figure  2  correspond  to  the  eight- 
year  period  1962-70,  including  the  1967-68  apparition  which  is 
of  most  interest  in  the  present  study.   Thus,  the  definitions 
of  the  sources  in  terms  of  central  meridian  longitude  (CML)  and 
the  "favored"  lo  ranges  from  superior  geocentric  conjunction 
(ISC)  are  drawn  by  inspection  from  this  plot.   Though  there  are 
small  variations  with  radio  frequency  evident  in  Leacock 's 
diagram,  the  values  in  Table  1  have  been  adopted  for  the  18  to 
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Figure  2.   Maps  of  the  occurrence  of  Jovian 

decametric  events  as  a  function  of 
CML  and  ISC.   (After  Leacock,  1971.) 


27  MHz  region  for. the  purposes  of  this  investigation. 

Table  1 
18  -  27  MHz  Source  Definitions 


Source 


CML,  System  III 
1957.0  (degrees) 


ISC  (degrees) 


200  -  290 

75  -  200 

285  -  40 


220  -  270 

60  -  130 

225  -  260 


Research  Tactics 

Most  of  the  work  that  has  been  devoted  to  the  decametric 
radio  emission  from  Jupiter  can  be  lumped  into  one  of  two  cate- 
gories.  The  first  of  these  categories  is  the  study  of  the 
intrinsic  properties  of  the  emissions  themselves.  Among  these 
properties  are  the  wave  polarization,  energy  flux,  dynamic 
spectra,  etc.   The  second  category  is  the  study  of  relation- 
ships between  some  property,  such  as  the  occurrence  of  the 
radiation,  and  other  events  in  the  outside  world.   Investiga- 
tions of  this  kind  ultimately  take  the  form  of  some  kind  of 
cross -correlation . 

This  characterization  of  research  tactics  betrays  the 
implicit  assumption  that  the  production  of  the  intermittent 
radiation  is  (at  least)  a  two-step  process.  The  first  step 
consists  of  having  the  necessary  physical  conditions  for  a 


radiation  mechanism  to  operate.   For  example,  these  might  be 
the  existence  of  a  magnetic  field  and  a  source  of  electrons 
that  might  conspire  to  produce  gyromagnetic  or  synchroton 
radiation. 

The  second  step  is  the  intermittent  availability  of  an 
"exciter,"  some  event  which  completes  the  chain  of  conditions 
required  to  set  off  the  radiation  mechanism.   For  example, 
recent  work  by  Krausche  et  al.  (1973) ,  Ellis  (1973) ,  and 
Krausche  (1975)  seems  to  suggest  that,  when  the  lo  related 
source  B  radiation  is  produced,  it  comes  from  electrons 
emitting  doppler  shifted  gyromagnetic  radiation  as  they  move 
upward  along  the  planet's  magnetic  flux  lines.   The  role  of 
the  exciter  in  this  process  is  played  by  lo.  According  to 
Goldreich  and  Lynden-Bell  (1969),  lo  energizes  magnetospheric 
electrons  while  driving  them  down  to  the  ionosphere  by  acting 
as  a  unipolar  electrical  generator.   When  lo  is  within  the 
favored  range  with  respect  to  a  source,  radio  emission  is 
stimulated.   Smith  (1975)  has  recently  suggested  that  subse- 
quent diffusion  of  electrons  back  into  the  void  swept  out  by 
the  lo  flux  tube  plays  a  role  in  the  non-Io  emission. 

The  Present  Study 

Underlying  Processes 

The  question  of  the  lo  related  emission  is  far  from 
solved  since,  among  other  things,  it  is  not  clear  why  the 
localized  sources  themselves  exist.   This  points  out  a  very 


fundamental  problem,  since  a  large  amount  of  the  decametric 
emission  does  not  correlate  with  lo,  but  almost  all  of  it 
correlates  with  the  defined  source  locations. 

Ignoring  the  question  of  what  is  special  about  the  source 
longitudes,  it  is  apparent  that  these  longitudes  are  somehow 
associated  with  locations  where  conditions  are  conducive  to  the 
emission  of  signals  receivable  at  the  earth.   If  one  considers 
only  the  non-Io  component,  a  natural  question  is  whether  this 
component  is  also  triggered,  at  least  in  part,  by  some  exciter 
which  also  has  origins  beyond  the  planet  itself.   Alternatively, 
one  should  also  ask  whether  there  are  unknown  systematic  influ- 
ences, such  as  beaming  effects,  which  affect  terrestrial  reception. 

In  view  of  the  success  of  the  lo  relationship  it  is 
reasonable  to  examine  the  remaining  three  Galilean  satellites. 
This  has  been  done  by  many  workers  with  no  sure  success; 
however,  a  more  thorough  discussion  of  this  point  is  reserved 
for  a  later  chapter. 

At  this  point,  a  matter  of  semantics  requires  discussion. 
The  principal  thrust  of  this  study  is  directed  at  the  question  of 
the  existance  of  non-Ioian  Jovian  exciters  or  other  reception 
influences.   In  this  regard,  a  distinction  between  the  term 
"mechanism"  and  "process"  is  made.   The  expression  "mechanism" 
will  be  reserved  to  indicate  the  immediate  source  of  the  emission 
(e.g.,  synchrotron,  Cherenkov,  etc.).   On  the  other  hand,  the 
term  "process"  will  refer  to  the  joint  entity  of  exciter  and 
radiation  mechanism,  though  this  will  frequently  be  treated  as 
a  "black  box"  operation. 


Randomness 

The  objective,  then,  is  to  examine  the  process  itself, 
and,  if  possible,  to  find  other  exciters.   In  order  to  do  this, 
it  is  important  to  examine  the  question  of  whether  the  non-Io 
component  is  a  random  process.   If  it  is  shown  not  to  be  random, 
then  the  existence  of  a  systematic  exciter  is  possible.   There 
are  statistical  tools  available  to  examine  the  randomness  of  a 
function  and  these  will  be  employed  to  deal  with  this  question. 

Jupiter-Solar  Correlation 

Historically,  this  investigation  began  as  an  effort  to 
see  if  the  non-Io  component  of  the  various  sources  was  related 
to  solar  wind  influences.   As  such,  considerable  attention  has 
been  given  to  the  sun,  solar  magnetic  sectors,  the  solar  wind 
and  a  model  for  solar  wind  velocity  stream  interactions. 

Subsequently,  it  was  recognized  that  in  order  to  establish 
any  such  correlation  one  must  also  consider  more  fundamental 
questions,  such  as  that  of  randomness.   The  result  was  a  broader 
assault  on  the  general  problem.   Even  so,  a  preoccupation  with 
periodic  processes  will  become  evident  in  what  follows.   This 
interest  is  motivated,  in  part,  by  the  prejudice  that  the  solar 
rotation  may  have  an  influence.   But,  the  concern  for  periodici- 
ties is  also  the  result  of  realizing  that  if  there  are  periodic 
exciters,  they  would  be  the  easiest  to  separate  from  the  noise 
in  the  data  due  to  their  inherent  redundancy. 


Methods  of  Study 

In  approaching  these  problems  an  examination  of  both 
Jovian  and  solar  wind  data  will  be  undertaken.   The  treatment 
of  these  data  will  consist  of  techniques  such  as  autocorrelation, 
occurrence  spectral  analysis  and  the  development  of  various 
distributions.   Thus,  each  class  of  data  will  be  viewed  by 
itself  to  discover  any  clues  that  may  be  there  to  be  found. 
Then,  the  two  sets  of  data  will  be  studied  with  respect  to  one 
another,  principally  by  way  of  Chree  (1912)  analysis  followed 
by  cross  correlation. 

There  are  already  a  number  of  known  or  strongly  suspected 
periodic  influences  on  the  Jovian  emissions,  six  to  be  exact. 
These  will  be  dealt  with  in  the  next  chapter,  but  the  point 
here  is  that  the  various  techniques  outlined  above  can  be 
strongly  influenced  by  periodic  processes.   Consequently,  each 
treatment  of  the  data  must  either  recognize  these  effects  in 
the  reduced  output,  or  be  contrived  in  such  a  way  that  the 
data  reduction  is  blind  to  the  known  effects.  As  will  be  shown, 
this  is  not  a  trivial  problem.   Indeed,  modern  statistics  is 
just  beginning  to  deal  with  the  problem  of  irregular  time 
series.   Specifically,  the  data  received  at  the  earth  are 
riddled  with  "holes"  caused  by  the  inability  of  a  station  to 
observe  a  given  source  continuously.   Thus,  the  table  of,  let 
us  say,  daily  numerical  measures  of  activity  has  the  potential 
for  significant  error  in  that  it  cannot  reflect  activity  that 
occurred  during  that  part  of  the  day  when  the  source  could  not 
be  seen.   Nevertheless,  these  problems  can  be  managed,  provided 


10 


one  does  not  lose  sight  of  their  existence. 

Choice  of  Apparition 

Another  related  problem  is  that  of  whether  or  not  the  pre- 
sumed, yet-to-be-discovered,  exciter  is  truly  periodic,  lasting 
for  ail  time,  or  whether  it  is  almost  periodic  and  exists  for 
a  time  less  than  that  of  the  study. 

Since  the  examination  of  a  possible  Jupiter-solar  relation- 
ship is  a  major  part  of  this  work,  the  above  problem  has  special 
significance.   In  particular,  the  solar  phenomenon  under  study 
is  the  control  of  the  solar  wind  apparently  exerted  by  the  solar 
magnetic  sectors.   These  sectors  can  be  relatively  stable  for 
a  time  on  the  order  of  months — but  not  much  longer.  With  this 
in  mind,  it  has  been  necessary  to  confine  the  study  to  single 
Jovian  apparitions. 

Specifically,  it  was  decided  to  restrict  the  study  to 
apparitions  when 

1.  The  solar  magnetic  sector  structure  was  simple 
and  stable. 

2.  Good  solar  wind  data  were  available. 

3.  Jovian  decametric  data  were  available. 

The  1967-68  apparition,  blessed  with  a  stable  two  sector 
sun,  was  the  first  to  meet  all  these  criteria  and  hence  has 
received  the  bulk  of  the  attention.   Very  recently  Svalgaard 
(1976)  has  published  sector  structure  data  which  showed  that 
the  1974  apparition  also  satisfies  these  criteria  and  this 
period  has  also  been  studied. 
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Source  of  Data 

The  Jovian  decametric  data  employed  in  this  research 
are  those  of  the  University  of  Florida  Radio  Observatory. 
Solar  magnetic  sector  structure  data  have  been  provided  by 
Wilcox  (1970,  1971)  and  Svalgaard  (1975).   The  solar  wind  data 
for  1967-68  were  gathered  by  Pioneer  VI  and  VII  in  solar  orbit 
near  1  AU,   essentially  in  the  ecliptic  plane,  and  published 
in  Solar-Geophysical  Data*.   The  solar  wind  velocity  data  for 
1974  are  from  Bame  et  al.  (1976). 


U.  S.  Department  of  Commerce,  Environmental  Sciences 
Services  Administration  periodical:   Government  Printing 
Office,  Washington,  D.  C. 


CHAPTER  2 
OBSERVATIONAL  PROBLEMS  -  THE  STROBES 

Periodic  Modulation 

As  mentioned  in  the  preceding  chapter,  there  are  six 
certain  or  nearly  certain  periodic  influences  on  the  observation 
of  Jovian  decametric  radiation.   These  influences  arise  from  the 
rotations  of  the  earth  and  Jupiter,  the  Jovian  synodic  and 
sidereal  revolutions,  the  earth's  tropical  or  seasonal  revolu- 
tion, and  finally,  the  revolution  of  lo. 

Rotation  of  Jupiter 

Since  the  three  classical  sources  are  of  limited  angular 
extent,  their  rotation  with  the  planet  presents  them  to  the 
earth  only  once  every  ten  hours  or  so.   The  lAU  definition  of 
the  sidereal  radio  rotation  period,  the  System  III  period,  is 

Vi    m    e 

9  55  29  .37,  though  later  work  by  Smith  et  al.  (1965),  Gulkis 
and  Carr  (1966)  and  by  Duncan  (1971)  indicates  that  the  true 
mean  radio  period  is  more  clearly  9'^55"'29^.  70  or  some  0.33 
seconds  longer.   However,  near  opposition  the  motion  of  the 
earth  can  reduce  the  synodic  radio  period  to  as  low  as  9^54"^49^. 

The  longitudinal  extent  of  each  of  the  sources  is  on  the 
order  of  90°,  and  therefore  the  presentation  of  a  source  to 
the  earth  takes  place  for  a  period  of  about  2.5  hours  during 


12 


13 


each  System  III  rotation. 

Rotation  of  Earth 

Since  the  earth  rotates  on  its  axis  with  a  sidereal 
period  of  about  23  56^^04^,  a  given  station  fixed  on  the  earth 
near  the  equator  can,  at  most,  observe  Jupiter  for  about  12 
hours  out  of  every  24.   Thus,  under  ideal  circumstances,  the 
odds  of  being  able  to  view  each  source  during  an  observation 
period  are  good;  still,  about  every  other  Jovian  rotation  goes 
unobserved. 

It  might  be  added  at  this  point  that,  of  course,  due  to 
the  orbital  motions  of  the  earth  and  Jupiter,  the  rotation 
period  of  the  earth  with  respect  to  Jupiter  is  not  generally 
the  sidereal  period  of  the  earth.   It  is  the  sidereal  value 
only  when  the  earth  is  at  greatest  eastern  or  western  elonga- 
tion as  viewed  from  Jupiter.  When  the  earth  is  at  inferior 
conjunction  (i.e.,  Jupiter  is  at  opposition)  the  period  is 
shortened  by  almost  one  minute  to  23^55"'08^,  while,  when  the 
earth  and  Jupiter  are  at  superior  conjunction  with  respect  to 
each  other,  the  period  is  lengthened  about  39  seconds  to 
23^56"'43^. 

Jovian  Synodic  Year 

The  above  variation  of  the  earth's  Jovicentric  rotation 
period,  and  its  counterpart  in  the  Jovian  rotation  period, 
are  both  related  to  Jupiter's  synodic  year.   However,  they  are 
higher  order  effects  or  pertubations  associated  with  the  rota- 
tion periods.   There  is  yet  another  factor  which  is  far  more 
potent. 
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Unlike  many  other  areas  of  radio  astronomy  which  use 
higher  frequencies,  work  at  decametric  wavelengths  is  strongly 
influenced  by  fluctuations  in  terrestial  ionospheric  trans- 
parency.  Ionization  levels  in  the  daytime  F  layer  and 
sporadically  in  the  nighttime  E  layer  may  reach  levels  which 
render  the  ionosphere  opaque  by  reflecting  decametric  radiation, 
impinging  on  the  earth,  away  into  space.   This  problem  is  com- 
pounded by  the  fact  that  opacity,  or  if  one  wishes,  reflectivity, 
is  proportional  to  the  secant  of  the  zenith  angle  for  an  earth- 
bound  observer.   Consequently,  as  the  level  of  ionization  builds 
up,  such  an  observer  will  generally  find  both  natural  and  man- 
made  radio  sources  from  beyond  the  horizon,  reflected  by  the 
ionosphere  at  large  zenith  angles  and  picked  up  by  the  antenna 
sidelobes,  producing  significant  interference.   This  terrestrial 
interference  will  frequently  terminate  observation  attempts 
even  before  smaller  zenith  angle  opacity  precludes  reception 
of  extraterrestrial  signals. 

It  is  possible  to  minimize  the  ionospheric  effects  either 
by  observing  at  higher  frequencies  where  the  transparency  is 
greater  or  by  employing  a  narrow-beam  antenna  with  strong  side- 
lobe  suppression.   For  various  reasons  it  is  generally  either 
undesirable  or  not  practical  to  take  these  steps,  though  the 
University  of  Florida's  640  dipole  array  at  26.3  MHz  is  an 
exception.   As  a  result,  most  Jovian  decametric  observations 
are  subjected  to  a  strong  diurnal  effect.  As  Jupiter  moves 
into  the  daytime  sky  during  its  synodic  cycle,  observation 
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becomes  progressively  more  limited,  usually  ceasing  altogether 
near  conjunction.   After  conjunction  the  situation  improves, 
reaching  its  peak  at  opposition.   The  net  effect  is  the  impo- 
sition of  a  modulation  with  a  period  equal  to  the  Jovian 
synodic  year. 

Jovian  Sidereal  Year 

Jupiter's  spin  axis  is  inclined  3° 07'  to  the  normal  of 
the  plane  of  its  orbit.   The  planet's  orbit,  in  turn,  is  in- 
clined 1**18'  to  the  ecliptic  plane.  These  two  effects  combine 
to  cause  the  Jovicentric  declination  of  the  earth,  D„,  to  vary 
about  +4"  over  a  Jovian  sidereal  year.   The  radiation  centroid 
and  abundance  for  Source  A  apparently  oscillate  with  a  period 
near  the  sidereal  period  of  11.86  years.   This  effect  has  been 
explained  by  Gulkis  and  Carr  (1966)  and  by  Donivan  and  Carr 
(1969)  as  the  result  of  the  source  A  emission  being  confined 
to  a  relatively  narrow  sheet.   Thus,  as  Jupiter  moves  about 
its  orbit  the  variation  of  D^,  directs  more  or  less  of  the  beam 
towards  the  earth.   The  last  peak  in  activity  due  to  this 
effect  was  about  1964.3  (Carr,  1972).   Presumably  the  following 
maximum  will  be  found  to  have  occurred  around  1976.2. 

Terrestrial  Tropical  Year 

The  sporadic  nighttime  E  layer  ionization  mentioned  above 
is  dependent  on  the  earth's  tropical  year.   Sporadic  E  levels 
tend  to  peak  during  the  summer  months,  causing  not  only  a 
diminution  of  extraterrestial  signals,  but  even  in  the  presence 
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of  locally  adequate  penetration,  ionization  clouds  near  the 
horizon  often  cause  electrical  interference — principally  light- 
ning— from  beyond  the  horizon  to  skip  into  the  side  lobes  of 
the  antenna.   The  E  layer  has  a  minor  peak  in  the  winter  months 
but  is  not  too  bothersome.   Not  only  is  the  probability  of 
opacity  much  less  but  the  lightning  problem  is  largely  a 
summertime  phenomenon  as  well. 

lo  Effect 

The  role  of  the  satellite  lo  in  the  Jovian  decametric 
process  has  been  outlined  in  Chapter  1.   lo  is  particularly 
effective  in  producing  source  B  activity,  at  least  that  which 
can  be  observed  from  the  earth. 

lo  revolves  about  Jupiter  with  a  period  of  about  42^28"^39^, 
while  Jupiter  rotates  under  it  with  a  synodic  period  of  12^57"'03^. 
There  is  evidence  to  suggest  that  the  principal  period  for  the 
lo  related  activity  is  the  13-hour  period.   Duncan  (1965)  re- 
ported an  occurrence  spectral  peak  at  12*^57"*08^.   Further, 
superposed  epoch  analysis  of  activity  using  this  period  revealed 
a  peak  when  lo  was  near  Jovian  longitude  200° — the  plane  thought 
to  contain  the  Jovian  magnetic  pole.   Duncan  points  out,  however, 
that  in  order  for  emission  to  be  beamed  at  the  earth,  CML  110° 
(source  B)  or  CML  250°  (source  A)  must  also  be  aimed  at  the 
earth.   He  then  suggests  that,  since  the  same  relative  angular 
separation  of  lo  and  sources  A  and  B  occurs  at  the  13 -hour  in- 
tervals, every  13  hours  the  source  emission  probability  is 
enhanced  irrespective  of  geocentric  or  heliocentric  central 
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meridian  longitude.   If  this  were  fact  it  would  certainly 
have  a  strong  effect  on  models  for  the  radiation  mechanism. 
A  review  of  Duncan's  analysis  by  this  author  does  not  reveal 
a  compelling  case  for  the  reality  of  the  13-hour  period  it- 
self, though  the  rest  of  his  conclusions  seem  valid. 

Kaiser  and  Alexander  (1973)  have  done  a  similar  study 
employing  17  years  of  data.   They  find  almost  six  times  as 
much  power  in  the  13-hour  peak  as  in  the  42-hour  period. 
However,  in  an  earlier  analysis  Bigg  (1966)  concluded  the 
opposite;  that  is,  that  the  13-hour  period  was  merely  a  mathe- 
matic  artifact  and  that  the  42-hour  period  was  dominant. 
Suffice  it  to  say  at  this  point  that  the  lo  effect  can  only 
influence  the  radiation  received  at  the  earth  when  the  source 
is  aimed  at  the  earth  and  lo  is  near  longitude  200".   This 
is  not  in  dispute  and  the  periodic  motion  of  lo  about  Jupiter 
can  be  considered  mathematically  independent,  as  can  also  the 
Jovian  rotation. 

Other  Modulations 

There  are  a  number  of  other  sources  of  unwanted  influences 
that  might,  in  effect,  contaminate  the  data.   Periodic  emissions 
of  data,  as  might  occur  from  the  observing  habits  of  regularly 
scheduled  observers,  will  cause  no  serious  problems  due  to  the 
nature  of  the  data  reduction  procedure.   Aperiodic  effects, 
particularly  random  ones,  will  also  have  little  effect.   An 
examination  of  the  data  has  not  revealed  any  substantive 
problems.   An  open  mind,  nevertheless,  must  be  maintained  in 
this  regard  and  cross  checks  on  the  data  will  be  made  and 


18 
evaluated  in  what  follows. 

Radiation  Process  Model 

The  Observation  Function 

The  preceding  sections  have  been  written  to  underscore 
the  fact  that  whether  or  not  one  observes  Jovian  decametric 
radiation  depends  strongly  on  several  known  periodic  influences. 
Even  though  some  of  these  influences  may  be  intrinsic  to  the 
generation  of  some  of  the  emissions,  as  in  the  lo  effect,  they 
can  all  be  treated  as  observational  artifacts.   They  are  known 
to  exist  and  it  is  not  wished  to  study  them  here.   Rather,  if 
they  can  be  removed,  what  remains  is  that  which  is  to  be 
studied  here. 

Since  this  examination  has  been  contrived  to  cover  single 
apparitions,  it  is  possible  to  ignore  several  of  the  known 
periodicities  (specifically,  those  effects  which  have  periods 
in  excess  of  one  apparition,  or  about  200  days).  This  elimi- 
nates synodic,  tropical  and  sidereal  year  influences  and 
leaves  the  rotation  of  Jupiter  and  the  earth,  and  the  revolution 
of  lo. 

It  is  now  wished  to  develop  a  mathematical  expression 
as  a  model  for  the  observations  themselves.   This  expression 
can  take  the  form  of  a  time  series  vdiich,  in  its  simplest  form, 
has  a  value  of  one  when  the  source  is  being  observed  and  a 
value  of  zero  when  it  is  not.   This  function  will  be  periodic, 
containing  components  corresponding  to  the  Jovian,  terrestrial 
and  loian  periods.   If  we  assume  that  these  three  components 
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are  independent  we  may  define  four  functions,  each  having  a 
value  of  one  when  the  condition  is  true  (e.g.,  the  Jovian  source 
is  aimed  at  the  earth)  and  a  value  of  zero  otherwise.   For  a 
given  source  S,  these  functions  are, 

Jg(t):  Is  the  source  aimed  at  the  earth? 

E  (t) :  Can  the  station  "see"  Jupiter? 

Ig{t):  Is  lo  in  the  favored  range? 

Ng(t):  Is  lo  out  of  the  favored  range? 

Of  course  Ng(t)  is  just  1  -  Ig(t).  The  composite  function  repre- 
senting whether  or  not  a  specific  source  is  being  observed  then 
is 

0  (t)  =  Jg(t)E_{t)I  (t)  for  lo  related  observations 

Og(t)  =  J  {t)E_{t)N_(t)  for  non-Io  related  observations, 

The  composite  observation  function  0(t)  will  have  a  value  of 
one  only  when  all  the  necessary  conditions  are  met,  that  is, 
the  source  S  is  aimed  at  the  earth,  the  earth  station  is  aimed 
at  Jupiter  and  listening,  and  lo  is  in  the  appropriate  ISC  range. 

The  Activity  Function 

Though  0(t)  expresses  whether  or  not  the  observer  has 
the  opportunity  to  observe  radiation,  it  in  no  way  guarantees 
that  such  emissions  will  be  seen.   For  example,  during  the 
1967-68  apparition  source  A  was  active  at  18  MHz  only  7.7% 
of  the  time  during  non-Io  observations  and  32%  of  the  time 
during  lo  observations  using  a  relatively  low  sensitivity  Yagi 
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array.  A  more  complete  view  of  activity  statistics  is  dis- 
cussed at  the  end  of  this  chapter,  including  a  comparison  of 
activity  indices  of  low  sensitivity  systems  versus  high  sensi- 
tivity systems.   The  principal  point  here  is  that  there  is  an 
activity  function,  analogous  to  the  observation  function,  which 
expresses  whether  or  not  emissions  are  detected.   Actually, 
there  are  two  activity  functions,  one  of  this  is  global,  that 
is,  an  intrinsic  function  which  is  defined  as  having  a  value 
of  one  whenever  a  given  source  emits,  with  lo  in  the  appropriate 
range,  regardless  of  whether  the  source  is  aimed  at  the  earth 
or  whether  the  earth  station  can  see  the  source,  and  a  value  of 
zero  otherwise.   The  other  activity  function  represents  merely 
the  observed  activity,  that  which  occurred  when  Oj,{t)  was  equal 
to  one. 

If  one  represents  the  global  or  intrinsic  activity  func- 
tion as  Ag(t)  and  the  observed  activity  function  as  A'(t),  they 
are  related  by 

A^(t)  =  Ag{t)Og(t) 

which  for  non-Io  activity  is 

A^(t)  =  Ag(t)Jg(t)Eg(t)Ng(t). 

While  the  intrinsic  function  A(t)  above  represents  the  entire 
activity  produced  by  the  underlying  processes,  the  earth-bound 
observer  is  permitted  only  a  tantalizing  view  of  this  process 
as  if  it  were  illuminated  by  a  strobe  lamp  driven  by  a  voltage 
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containing  the  constituent  frequencies  of  0(t). 

In  principle,  then,  given  adequate  data  one  might  hope 
to  find  the  intrinsic  function  A  (t)  by  the  simple  expediency 


A' (t) 
A,(t)  =    ^ 


S        0(t) 

Effect  on  Jovian  Data 

The  Raw  Time  Series 

The  previous  section  has  indicated  the  direction  of 
this  program:  the  development  of  procedures  that  might  permit 
the  separation  of  the  "true"  activity  from  the  artifact  modu- 
lated activity.   Such  procedures,  however,  are  constrained  to 
dealing  with  the  actual  raw  data. 

While  the  concept  of  Ag{t)  expressed  above  is  simple 
enough,  the  fact  is  that  the  various  observational  functions 
are  not  known  analytically.   It  is  true  that  in  principle  they 
can  be,  but  for  practical  purposes  it  is  far  more  important  to 
deal  with  the  actual  functions  as  they  really  evolved.  As  such, 
they  are  dependent  on  many  transient  influences:  the  occurrence 
of  thunderstorms,  the  inability  of  an  observer  to  stand  his 
watch,  the  presence  of  hunters  with  walkie-talkies,  arid  other 
earth-bound  stations  on  nearby  frequencies  providing  inter- 
ference, etc. 

The  observed  activity  time  series  will  be  developed  from 
what  the  Florida  group  calls  an  "intermediate  deck."  This  is 
a  catalog  of  Jovian  decametric  radio  observations  which  reports 
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the  CML  at  the  beginning  and  end  of  each  observation,  and  at 
the  beginning  and  end  of  any  observed  activity,  the  correspond- 
ing position  of  lo  (ISC)  and  so  forth.   All  these  angular 
measures  are  generally  reported  in  blocks  or  "zones"  of  5°. 
Thus,  0*>-5°  CML  corresponds  to  zone  1,  while  zone  72  represents 
355°-360°  CML. 

Since  the  intermediate  deck  reports  observational  data 
independently  of  activity  data,  it  will  also  be  possible  to 
develop  the  corresponding  raw  observation  time  series.   The 
observation  series  is  somewhat  model  dependent.  Nevertheless, 
as  previously  stated,  the  probability  of  observation  is  assumed 
to  be  constantly  one  or  constantly  zero.   This  is  a  simplifi- 
cation, since  the  distribution  of  activity  within  a  given  source, 
such  as  is  shown  in  Figure  1,  suggests  that  Jc(t)  is  more  nearly 
gaussian.   In  all  probability,  the  incidence  of  reception  at 
an  earth-bound  observatory  is  highest  at  transit  and  lowest 
near  the  horizon,  with  a  corresponding  effect  on  E(t).   The 
precise  shape  of  all  of  the  component  functions  is  not  known, 
however,  and  the  rectangular  approximation  seems  a  good  com- 
promise at  this  point.   At  worst,  it  somewhat  exaggerates  the 
weight  of  events  that  occur  near  the  limits  of  one  of  the 
component  functions . 

At  this  point,  we  have  two  raw  time  series  whose  function 
values  are  either  zero  or  one.   The  problem  then  is  to  try  to 
compensate  the  activity  function  for  the  observational  effects. 
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Compensation  Procedure 

It  has  been  stated  earlier  that,  in  principle,  the  under- 
lying activity  function,  Ag,  could  be  known  from 


0(t) 

However,  as  0{t)  is  often  zero,  the  mathematical  expression  is 
undefined  over  much  of  its  range.   Physically,  this  means  that 
one  cannot  say  whether  or  not  there  is  an  underlying  activity 
when  the  source  was  not  observed.  On  the  other  hand,  A„(t) 
does  have  some  real  values.  Whenever  0(t)  is  one,  the  expres- 
sion yields  meaningful  numbers. 

Due  to  the  desire  to  explore  specifically  the  possibility 
of  solar  correlation  and  the  fact  that  solar  wind  data  were 
generally  confined  to  one  sample  per  day,  the  decision  was  made 
to  adjust  the  available  values  of  A  (t)  to  one  value  per  day. 
Thus,  the  intrinsic  activity  function  was  chosen  to  be 
estimated  by 

11 

where  A' (t).  and  0(t) .  are  expressed  in  5°  zones.   The  activity 
function,  then,  is  nothing  more  than  a  daily  activity  index 
expressed  as  the  ratio  of  the  number  of  source  zones  seen 
storming  to  the  number  of  source  zones  observed  each  day. 


24 


A  table  of  values  for  A^lt)  will  have  a  real  value  for 
most  days,  since  each  source  is  usually  observed,  at  least  in 
part,  every  day.   Nevertheless,  as  often  as  one  day  in  four  for 
a  source  may  be  skipped. 

Influence  of  the  Compensation 

The  question  naturally  arises,  then,  as  to  what  effect 
these  covert  and  overt  indeterminacies  will  have  on  any  attempt 
to  use  Ag(t)  to  study  properties  of  the  presumed  underlying 
phenomena.   Since  the  time  series  will  become  the  input  to 
autocorrelation  and  power  spectrum  routines,  it  was  indeed 
disappointing  to  discover  that  the  bulk  of  current  statistical 
theory  was  predicated  on  the  existence  of  evenly  spaced,  un- 
interrupted time  series. 

In  the  present  case,  the  situation  is  further  complicated 
by  the  fact  that  the  "samples"  in  the  time  series  are  themselves 
composed  of  samples.   That  is,  the  daily  values  of  A(t)  are  calcu- 
lated from  a  sample  of  the  total  activity  for  that  day.   The 
minimum  valid  sample  size  for  a  given  observation  is  one  5° 
zone,  while  the  maximum  sample  size  may  exceed  20  zones. 

The  number  of  zones  sampled  on  a  given  day,  as  expected, 
will  be  shown  to  be  a  somewhat  complicated  periodic  function  of 
time.  As  a  result,  the  sample  size  from  which  the  daily 
values  of  A(t)  are  derived  is  variable.   This  means,  among 
other  things,  that  the  weight  of  each  daily  number  is  not  con- 
stant.  Thus,  from  time  to  time,  when  there  are  only  a  few 
samples,  a  value  of  A(t)  will  appear  which  is  inordinately 
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high  or  low,  and  yet  will  be  considered  on  a  even  par  with 
other  values  drawn  from  larger  samples  more  likely  to  approach 
the  population  mean. 

A  subsequent  chapter  contains  a  discussion  of  the  various 
statistical  properties  of  the  data,  and  it  will  be  shown  that 
about  90%  of  the  non-zero  values  for  0(t)  contained  5  or  more 
zones,  while  only  7%  of  the  observed  storms  lasted  longer  than 
five  zones.   Consequently,  the  errors  introduced  by  varying 
sample  size  are  not  expected  to  be  serious. 

Having  calculated  a  table  of  the  daily  values  of  A(t), 
some  of  which  are  indeterminate,  the  question  then  arises  as 
to  what  effect  these  indeterminacies  will  have  on,  say,  the 
power  spectrum  of  the  time  series.   Jones  (1962)  has  studied 
this  problem  for  data  with  "regularly  missed"  observations. 
In  his  work  he  assumed  that  data  were  collected  at  equally  spaced 
intervals,  but  that  periodically  observations  had  to  be  skipped 
for  whatever  reason. 

In  making  the  spectral  analysis  he  used  the  common  tech- 
nique of  first  calculating  the  autocorrelation  function  and 
then  taking  its  Fourier  transform.   However,  in  the  autocorre- 
lation he  included  only  pair  products  where  both  terms  existed 
and  simply  excluded  those  which  did  not  contain  both  terms. 
He  then  showed  that,  using  this  technique,  the  error  in 
estimating  the  spectral  density  was  not  greatly  increased. 
Figure  3  shows  his  plot  of  the  upper  bound  of  the  increase  in 
covariance  for  a  repeated  sequence  of  a  observations  followed 
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0.2  - 


Figure  3.   Effect  of  missing  data  on  power  spectrum  esti- 
mates,  a  is  the  upper  limit  for  the  standard 
deviation  of  the  spectriam  of  a  time  series 
missing  one  point  every  a  points.   a»j  is  the 
standard  deviation  for  N  actual  and  missed 
points,  while  Cj,^  is  the  standard  deviation 
for  M  actual  points.   (After  Jones,  1962.) 
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by  one  skipped  observation.   It  shows,  for  example,  that  if 
every  fourth  observation  is  missed  the  standard  deviation  of 
the  spectral  density  will  increase  by  at  most  25%. 

Jovian  Data  Reduction 

Zone  Filters 

The  Jovian  data  reduction  process  in  this  study  begins 
by  acting  on  the  Florida  intermediate  deck,  in  either  card  or 
magnetic  tape  form  with  a  computer  program  designated  as  Jupzone 
(see  Appendix) .   This  routine  consists  of  a  series  of  filters 
which,  first,  selects  only  intermediate  deck  records  which  are 
within  the  desired  apparition.   Second,  it  rejects  any  of  those 
records  which  do  not  represent  the  observations  of  the  chosen 
Jovian  source.  Third,  it  sorts  the  accepted  records  into  four 
categories: 

1.  lo  related  activity, 

2.  lo  related  observations, 

3.  Non-Io  related  activity, 

4 .  Non-Io  related  observations . 

Finally,  the  program  reduces  the  records  in  each  of  these  cate- 
gories to  a  single  entry  per  day — the  total  number  of  5°  zones 
in  that  category  for  that  day. 

The  output  from  the  Jupzone  program  consists  of  the 
(abbreviated)  Julian  date  and  a  list  of  the  four  zone  counts 
for  each  day,  for  each  specified  source,  during  the  given 
apparition.   This  information  is  produced  in  card  form  and 
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serves  as  the  input  to  a  variety  of  other  analysis  routines 
which  will  be  described  later.   Table  2  shows  a  printed  example 
of  a  segment  of  this  output. 

The  Reduced  Time  Series 

It  should  be  clear  that  the  time  series  generated  by 
Jupzone  can  easily  be  used  to  produce  the  function  A  (t)  for 
both  the  non-Io  and  lo  related  components.   While  much  of  the 
ultimate  data  analysis  does  deal  directly  with  A  (t) ,  some 
routines  require  the  separate  listing  of  activity  and  obser- 
vation.  Consequently,  the  series  are  left  in  the  separated 
form  at  this  level,  enabling  subsequent  programs  to  choose 
whether  or  not  to  calculate  Ag(t).  For  example,  by  leaving 
the  series  separate  it  was  possible  to  suppress  the  effect  of 
sample  size  biasing  in  reporting  the  activity  index  for  an 
entire  apparition,  Ag(t),  in  the  following  section.   In  this 
process,  separate  sums  of  the  daily  zone  counts  were  made 
over  the  entire  apparition , and  then  their  ratio  was  taken, 
rather  than  taking  the  average  of  the  daily  ratios. 

Gross  Statistics 

Radio  Frequencies 

During  the  1967-68  apparition  the  reception  frequencies 
employed  in  this  study  were  18,  22,  and  27  MHz.   While  the  Uni- 
versity of  Florida  operates  receiving  systems  outside  this 
range,  the  source  definitions  in  the  18  to  27  MHz  region  are 
essentially  identical.  At  the  same  time,  data  are  more 
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Table  2 
Example  of  Filtered  Jupiter  Data 
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consistently  available  in  this  range  than  at  lower  frequencies 
where  the  various  ionospheric  effects  become  increasingly  more 
difficult. 

For  the  study  of  the  1974  apparition  the  same  choices 
in  frequency  were  made.   However,  data  were  also  available 
from  the  640  dipole  array  at  26.3  MHz.   This  filled  aperture 
array  has  been  described  extensively  by  Desch  (1976)  and  offers 
a  sensitivity  some  three  orders  of  magnitude  greater  than  con- 
ventional Yagi  antennas. 

Activity  Indices 

Table  3  reports  the  values  of  the  activity  function, 
A_(t),  where  the  summation  is  over  the  entire  apparition  rather 
than  on  a  once-a-day  basis.   The  Table  shows  that  in  1967-68 
and  in  1974  with  the  standard  low  sensitivity  systems,  the 
values  of  the  lo  related  activity  indices  are  within  an  order 
of  magnitude  for  all  sources  at  a  given  frequency.   It  is  also 
evident  that,  for  these  systems,  there  is  a  general  inverse 
relationship  between  frequency  and  observed  activity. 

It  is  interesting  to  note  that  the  non-Io  component  is 
far  more  dependent  on  which  source  is  observed  than  is  the  lo 
related  component.   For  the  low  sensitivity  systems,  the  non-Io 
activity  indices  for  source  B  are  2  to  3  times  less  than  sources 
A  and  C  at  18  and  22  MHz. 

With  respect  to  the  comparison  of  1967-68  and  1974  it  is 
pointed  out  that  in  1974  non-Io  activity  was  observed  about 
twice  as  often  as  it  was  during  1967-68.   This  is  entirely 
consistent  with  the  D-  hypothesis  since  the  earlier  apparition 
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occurred  near  the  minimum  of  the  presumed  effect  while  the 
1974  apparition  was  more  nearly  at  the  expected  1976  maximum. 

Perhaps  the  most  dramatic  effect  shown  in  Table  2  is 
the  impact  of  the  high  sensitivity  system  at  26.3  MHz  with 
respect  to  the  established  low  sensitivity  patterns.   The 
"big  array"  shows  that  Jupiter  is  a  far  more  prolific  emitter 
than  the  conventional  receiving  systems  have  indicated.   Not 
only  is  the  non-Io  related  component  more  frequently  active, 
but  the  non-Io. related  activity  is  anywhere  from  3  to  30 
times  more  prominent. 


CHAPTER  3 
SOLAR  WIND 

Steady  Flow  Models 

Early  Clues 

Perhaps  the  earliest  recognized  clues  to  the  existence 
of  a  wind  of  particles  blowing  out  from  the  sun  were  corre- 
lations between  certain  geomagnetic  events  and  solar  phenomena. 
By  the  early  1930's,  two  such  classes  of  correlated  magnetic 
disturbances  were  identified.   Chapman  and  his  collaborators 
(Chapman,  1963)  established  the  connection  between  solar  flares 
and  the  sudden  commencement  (SC)  type  of  geomagnetic  storm. 
They  noted  that  many  large  solar  flares  were  followed  some 
12  to  36  hours  later  by  a  SC.   Since  it  was  clear  that  the 
assumed  exciter  was  either  electromagnetic  or  particulate, 
the  time  delay  was  due  either  to  some  kind  of  hysteresis 
effect  in  the  geomagnetic  reaction  to  an  exciter  propagating 
at  the  speed  of  light,  or  it  was  due  to  the  transit  time  of 
a  particulate  exciter  moving  at  1000  to  3000  km/s.   The 
hypothesis  which  evolved  held  that  a  cloud  of  ionized  particu- 
late ejecta  from  the  flare  compressed  the  geomagnetic  field 
sharply  when  the  two  collided  producing  the  rapid  onset  of 
the  observed  increase  in  the  horizontal  component  of  the 
field.   The  development  of  a  ring  current  due  to  the  subsequent 
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trapping  of  some  of  the  ejecta  accounted  for  the  depression 
in  the  field  which  followed  in  a  day  or  so. 

Taken  by  itself,  the  flare-SC  relationship  did  not  indi- 
cate whether  the  wind  was  a  continuous  process  or  just  a 
transient  one  associated  only  with  flares.   However,  Bartels 
(1932)  reached  the  conclusion  that  there  was  a  long-term 
process  in  his  study  of  "M-region  events".   These  events 
represented  a  second  class  of  geomagnetic  storms  which,  rather 
than  being  transient,  frequently  recurred  with  a  period  of 
about  27  days.  These  storms  showed  no  particular  association 
with  active  regions  of  the  sun,  but  Bartels  correctly  identi- 
fied their  period  as  being  due  to  persistant  streams  of  par- 
ticles emitted  at  low  latitudes  on  the  rotating  sun. 

Another  clue  turned  up  years  later  in  the  study  of 
cometary  tails.   It  was  noted  that  in  tails  containing  large 
concentrations  of  heavy  ions  such  as  N^^  and  CO-'*'  (rather 
than  dust)  radiation  pressure  was  inadequate  to  push  the  tails 
away  from  the  sum.   Nevertheless,  the  "heavy"  tails  pointed 
away  from  the  sun  just  as  did  the  lighter  dust  tails. 
Biermann  (1951)  suggested  that  the  "heavy"  tails  were  being 
pushed  away  by  collisions  between  the  ions  in  the  tail  and 
particles  continuously  flowing  out  from  sun  with  velocities 
of  500  to  1500  km/s. 

Hydrodynamic  Models 

Biermann 's  work  motivated  Parker  (1958)  to  study  the 
theoretical  feasibility  of  such  an  outflow.   He  recognized 
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that  for  any  reasonable  temperature  in  the  corona  the  thermal 
velocity  of  hydrogen  ions  was  no  more  than  about  one-half 
the  500  km/s  solar  escape  velocity  at  that  level.   It  was  not 
believed  that  even  the  high-end  tail  of  the  velocity  distri- 
bution could  account  for  the  number  of  particles  thought  to 
be  escaping.   Parker  approached  the  problem  in  two  ways.   His 
first  tack  was  to  examine  whether  or  not  hydrostatic  equi- 
librium was  possible.   In  doing  so  he  assumed  the  sun  was  a 
spherically  symmetric  gravitating  body  whose  surface — the 
outer  corona — was  at  a  given  temperature.  He  further  assumed 
that  above  this  surface  there  were  no  sources  or  sinks  of  heat, 
that  mass  was  conserved,  that  all  flow  was  strictly  radial  and 
that  there  was  no  magnetic  field.   He  then  invoked  the  baro- 
metric equation  of  hydrostatic  equilibrium: 

Ap(r)  +  g(r)  =  0 

where  p  is  the  thermal  pressure,  g  is  the  force  due  to  gravity 
and  r  is  the  radial  distance.   By  solving  this  relationship 
for  p(r)  he  was  able  to  show  that  the  pressure  at  infinity 
was  not  zero,  in  fact,  it  was  several  orders  of  magnitude 
higher  than  estimates  of  the  pressure  of  the  interstellar 
medium.   Hence,  he  concluded  that  there  must  be  an  effusion 
of  plasma — a  solar  wind. 
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Having  proved  the  failure  of  static  equilibrium,  his 
second  step  was  to  modify  the  barometric  equation  by  adding 
a  term  to  account  for  bulk  motion,  giving 

Vp(r)  +  g(r)  =  -VK(r) 

where  K  is  the  bulk  kinetic  energy.   He  proceeded  to  solve 
this  equation  for  the  solar  wind  bulk  velocity  u(r).   Remark- 
ably he  found  that  his  hydrodynamic  model  predicted  that  the 
velocity  was  a  monotonically  increasing  function  of  radial 
distance.   In  other  words,  despite  the  solar  graviational 
field,  the  wind  continued  to  accelerate  as  depicted  in 
Figure  4.   Parker  goes  on  to  show  that  mathematically  this 
unexpected  result  stems  from  the  fact  that  the  expansion 
occurs  in  three  dimensions  rather  than  one.   If  the  model  had 
assumed  an  infinite  flat  sun,  the  wind  velocity  would  decrease 
monotonically  with  distance. 

When  Parker's  paper  was  published  there  were  no  direct 
measurements  of  the  solar  wind.   The  earliest  direct  measure- 
ments would  not  occur  until  the  following  year,  1959,  when 
the  Russian  probe  Luna  1  was  flown.   There  were  other, 
essentially  contemporary,  theories  such  as  those  of  Chapman 
(1954)  and  Chamberlain  (1960)  which,  proceeding  from  different 
assumptions,  held  that  there  was  little  or  no  wind,  though 
interplanetary  space  was  filled  with  matter  more  or  less  in 
hydrostatic  equilibrium.   However,  once  the  direct  measurements 
were  obtained  it  was  clear  that  Parker's  assumptions  were  more 
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Figure  4.   Solar  wind  bulk  velocity  u(r).   (After 
Parker,  1958.) 
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nearly  correct.   With  the  launch  of  Mariner  2  for  Venus  in 
1962,  extended  measurements  became  available  (Neugebauer  and 
Snyder,  1966)  and  the  number  of  models  proliferated  rapidly. 

The  basic  requirements  to  devise  a  model  are  a  set  of 
conservation  equations  for  mass,  momentum,  and  energy  together 
with  appropriate  boundary  conditions.   Both  the  form  of  the 
equations  and  the  boundary  conditions  are  influenced  by  the 
assumption  of  some  kind  of  heating  mechanism  which  acts  on 
the  wind  somewhere  between  2  and  20  solar  radii.   This  process 
is  the  primary  energizer  for  the  wind. 

The  nature  of  this  heating  mechanism  was,  and  still  is, 
poorly  understood.   It  is  also  unfortunate  that  the  number  of 
terms  that  might  conceivably  be  included  in  the  conservation 
equations  is  large.   For  example,  the  energy  equation  might 
have  significant  terms  representing  gravity,  bulk  velocity, 
thermal  pressure,  magnetic  pressure,  thermal  conduction, 
thermal  radiation,  viscosity,  magnetohydrodynamic  radiation 
and  damping,  plasma  oscillation  radiation,  state  changes, 
species  interactions  and  so  forth.   A  great  deal  of  the  model 
building  that  has  been  done  has  involved  modifying  the  heat- 
ing mechanism  assumption  by  selecting  which  of  the  terms  to 
retain  in  the  conservation  equations  and  which  to  delete  as 
insignificant.   Parker,  for  example,  used  only  gravity, 
thermal  pressure  and  bulk  kinetic  energy  terms  in  his  energy 
equation;  Noble  and  Scarf  (1963)  added  terms  for  thermal 
conduction  and  viscosity.   Whang  and  Chang  (1965)  deleted 


39 

viscosity  and  tried  convection  of  kinetic  energy;^  Hartle 
and  Sturrock  (1968)  deleted  convection  and  argued  that  the 
plasma,  being  fully  ionized  and  essentially  collisionless, 
can  be  regarded  as  a  two  fluid  gas — protons  and  electrons — 
and  thus  each  species  can  have  a  different  temperature.   They 
made  the  appropriate  change  in  the  equations  and  produced 
another  model.   Subsequently,  Yeh  (1970)  developed  a  three 
fluid  model  by  adding  alpha  particles,  while  Cuperman  and 
Metzler  (1975)  considered  a  collisional  three  fluid  model 
with  provision  for  turbulence  heating. 

It  is  surprising  that  so  many  models  have  been  developed 
which  ignore  the  influence  of  the  magnetic  field.   One  would 
think  that  the  field  would  be  a  powerful  agent  in  the  lower 
corona  where  the  primary  heating  of  the  wind  plasma  occurs. 
Furthermore,  even  far  from  the  sun  heat  conduction  in  a  fully 
ionized  plasma  would  only  take  place  parallel  to  the  field. 
It  is  even  more  surprising  that  the  field-free  models  actually 
predict  the  various  wind  parameters  fairly  accurately,  though 
they  consistently  underestimate  the  average  bulk  velocity 
at  1  AU  by  20  to  30%.   The  results  of  a  number  of  steady  flow 
models  are  presented  in  Table  4,  together  with  various  experi- 
mental findings.   In  this  table  the  value  reported  for  Parker 
(1958)  corresponds  to  a  coronal  temperature  of  8  x  10  K°  at 
14  solar  radii. 

In  reviewing  the  various  models  it  is  difficult  to 
understand  why  most  modelers  appeared  content  with  bulk 
velocities  in  the  300-330  km/s  range.   Snyder  etal.  (1963), 
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Table  4 


Comparison  of  Theoretical  and  Experimental 
Solar  Wind  Parameters  at  1  AU 


u       n         T      B 
Model  (B=0)  (km/s)    (cm-3)    (105k°)    (y) 


Parker  (1958) 
Noble  &  Scarf  (1963) 
Whang  &  Chang  (1965) 
Hartle  &  Sturrock  (1968) 
Cuperman  s;  Metzler  (1975) 

Model  (Bt^O) 

Urch  (1969) 
Wolff  et  al.  (1971) 
Whang  (1971) 
Whang  (1972) 

Observation 


350 

352 

7 

2.8 

260 

8 

1.6 

250 

15 

3.5 

318 

10 

1.5 

371 

8 

4.39 

3.7 

303 

9 

2.0* 

8 

302 

5 

1.5 

5.8 

331 

6 

1.5* 

6.3 

Snyder  et  al.  (1963) 

504 

Wilcox  &  Ness  (1965) 

315 

9 

Hundhausen  (1970) 

325 

4       1 

Gosling  et  al.  (1971) 

426 

1.4^ 


''Electron  temperature 
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the  review  by  Ness  (1968)  and  Gosling  et  al.  (1971),  all 
clearly  show  that  the  mean  value  over  a  solar  rotation  was 
greater  than  400  km/s.   Most  modelers  emphasized  that  their 
models  were  for  the  quiet  sun,  but  1962  was  near  the  minimum 
of  Cycle  20  and  the  values  reported  by  Mariner  2  averaged  504 
km/s.   The  300  km/s  figure  is  more  nearly  the  lowest  observed 
value  published  after  1963. 

MHD  Models 

Formal  inclusion  of  the  magnetic  field  in  a  steady  flow 
model  seems  to  have  been  done  first  by  Weber  and  Davis  (1967) . 
The  one  fluid  inviscid  conduction  model  of  Urch  (1969)  was 
the  first  to  resolve  the  bulk  velocity  discrepancy,  but  it 
overestimated  the  proton  temperature  and  underestimated  the 
azimuthal  velocity,  both  by  an  order  of  magnitude.   Wolff 
et  al .  (1S71)  produced  a  viscous  two-fluid  conduction  MHD 
model  with  the  same  basic  discrepancies  as  Urch. 

In  an  effort  to  resolve  the  velocity  problem  in  his 
earlier  field-free  model  (Whang  and  Chang,  1965)  with  its 
assumed  strictly  radial  flow,  Whang  (1971)  proposed  an  MHD 
model  which  included  a  process. by  which  magnetic  field 
energy  was  continuously  converted  into  kinetic  energy  as  a 
by-product  of  the  expansion  process.   His  calculations  indi- 
cated that  this  process  would  increase  the  velocity  at  i  AU 
by  17%.   He  subsequently  extended  this  work  to  a  two  fluid 
model  (Whang,  1972).   However,  Barnes  (1974)  showed  that  in 
MHD  models  the  assumption  of  strictly  radial  flow  is  ho  longer 
valid.   Near  the  sun  the  intense  magnetic  field  attempts  to 
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impose  corotation  on  the  plasma.   The  resultant  azimuthal 
acceleration  imparts  an  azimuthal  velocity  component  which  can 
no  longer  be  ignored.   The  net  result  was  that  Whang's  radial 
velocity  increase  was  reduced  to  about  1%  at  1  AU. 

Recently  Nerney  and  Suess  (1975)  presented  a  model  which 
assumed  symmetry  about  the  sun's  spin  axis  only.   Their  model 
shows  a  weak  meridional  magnetic  pressure  gradient.   The  result- 
ing force  pushes  plasma  out  of  the  equatorial  plane  and  toward 
the  poles.   This  transport  process  imposes  a  -2.6  km/s  meridional 
velocity  component  on  the  solar  wind  near  1  AU.   They  suggest 
that  as  much  as  20%  of  the  equatorial  plasma  may  be  displaced 
from  the  solar  equatorial  plane  at  5  AU. 

Comments 

The  underestimation  of  the  bulk  velocity  of  the  solar 
wind  remains  a  feature  of  most  models.   Consequently,  the 
observed  bulk  kinetic  energy  is  on  the  order  of  60%  greater 
than  that  for  which  the  models  can  account.   Since  95%  or  more 
of  the  total  energy  flux  is  due  to  bulk  kinetic  energy,  a 
major  source  of  the  wind's  energy  is  still  not  appreciated. 

A  universal  feature  of  the  models  considered  is  the  call 
for  a  monotonically  increasing  bulk  velocity  with  radial  dis- 
tance.  Taken  at  face  value  this  would  imply  that  one  should 
expect  the  solar  wind  to  be  moving  faster  at  Jupiter  than  at 
the  earth.   This  increase  varies,  from  3  to  29%  depending  on 
the  model.   It  also  represents  an  extrapolation  of  models 
which  have  been  "tuned"  to  fit  known  conditions  in  the  corona 
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and  at  1  AU. 

Measurements  made  by  Mariner  2  failed  to  find  any 
decrease  in  wind  velocity  as  it  moved  in  toward  Venus 
(Snyder  et_al . ,  1963).   Likewise,  the  expected  increase  in 
velocity  was  not  seen  in  comparative  measurements  made  by 
Pioneers  X  and  XI  in  the  vicinity  of  4.5  AU  and  1.5  AU 
respectively  (Collard  amd  Wolfe,  1974).   However,  since 
the  models  have  steady  flow  as  a  basic  assumption  and  the 
actual  flow  is  quite  variable,  the  Mariner  and  Pioneer 
measurements  may  not  be  directly  comparable  to  velocities 
in  the  steady  flow  models.   At  the  very  least  it  seems  that 
there  should  be  some  kind  of  weighting  in  the  averaging 
process  to  enhance  the  importance  of  the  velocity  on  days 
when  the  particle  density  is  high  and  diminish  the  impor- 
tance of  the  velocity  on  days  when  the  particle  density  is 
low. 

The  Observed  Wind 

Rather  than  being  uniform  and  symmetrical,  the  solar 
wind  is  found  to  have  large  scale  structure  affecting  all  of 
its  measured  parameters.   At  times  this  structure  is  very 
regular  and  predictable,  while  on  other  occasions  it  appears 
chaotic. 

Bart el s '  Streams 

After  a  preliminary  review  of  the  first  4.5  months  of 
Mariner  2's  solar  wind  data, Snyder  et  al.  (1963)  reported 
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that  the  solar  wind  velocity  was  deeply  modulated,  containing 
four  or  five  major  peaks  every  27  days.   Furthermore,  the 
velocity  peaks  strongly  correlated  with  peaks  in  the  geo- 
magnetic index  K  .   Figure  5  shows  a  segment  of  their  data. 
The  inescapable  conclusion  was  that  Bartels'  picture 
of  the  recurrent  "M-region"  geomagnetic  disturbances  being 
due  to  the  existence  of  long-lived  high-flux  streams  of  solar 
particles  was  indeed  correct.   It  would  turn  out,  in  fact, 
that  geomagnetic  conditions  were  so  well  correlated  with  the 
solar  wind  that  they  could  be  used  to  extrapolate  back  to  as 
early  as  1932  to  produce  estimates  of  the  daily  orientation  of 
the  interplanetary  magnetic  field  (IMF)  as  has  been  done  by 
Svalgaard  (1976). 

Interplanetary  Magnetic  Sectors 

In  late  1963  the  IMP  satellite  was  launched  into  a 
very  eccentric  orbit  which  carried  it  into  the  interplanetary 
medium.  From  its  sensors  Wilcox  and  Ness  (1965)  were  able  to 
determine  that  the  magnetic  field  intensity,  particle  density 
and  particle  flux  were  also  strongly  correlated  with  the 
high  velocity  streams.   During  the  three-month  period  of  the 
satellite's  observations  they  found  that  there  were  four 
high-velocity  streams  in  each  27  day  solar  rotation.   As  the 
sun  rotated  three  of  these  streams  were  seen  at  intervals  of 
about  7.7  days  and  the  fourth  stream  3.9  days  later.   The 
width  of  the  three  wider  streams  was  about  4  days  followed 
by  3  to  4  days  of  slow  wind  speed. 
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6  10     20 
NOV. 


Figure  5.   Solar  wind  bulk  velocity  from  Mariner  2 

and  geomagnetic  index  ZK  .   (After  Snyder 
et  al.,  1963.)         P 
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The  observations  of  the  magnetic  field  showed  that 
during  the  entire  period  from  one  fast  stream  to  the  next, 
the  field  direction  was  either  directed  inward  or  outward 
from  the  sun.   However,  just  before  the  next  fast  stream, 
the  magnetic  field  direction  reversed  itself.   Consequently, 
the  plasma  was  organized  into  four  regions  or  "magnetic 
sectors"  of  alternately  inward  and  outward  directed  magnetic 
field.   The  boundary  between  two  successive  sectors  is  very 
abrupt.   The  field  reversal  at  the  boundary  was  seen  to  occur 
in  less  than  5  minutes.   Figure  6  shows  this  four  sector 
situation  as  it  was  observed  a  year  later  by  IMP  2  (Wilcox 
and  Colburn,  1969) . 

Extended  observations  have  shown  that  the  two  sector 
and  four  sector  configurations  are  those  most  frequently  seen 
(Svalgaard,  1976).   When  these  occur  the  pattern  usually  re- 
mains stable  for  several  months  to  a  year.  At  other  times 
the  pattern  may  be  very  confused,  having  as  many  as  twelve 
sectors  of  relatively  short  life  leading  to  rapidly  changing 
configurations. 

The  question  of  solar  cycle  effects  on  sectors  has 
been  examined  extensively  by  Wilcox  and  Colburn  (1969,  1970), 
Intriligator  (1974),  Vasyliunas  (1975)  and  others  with  mixed 
results.   Most  recently  Bame  et  al.  (1976)  have  concluded  that 
high  amplitude  streams  with  high  peak  velocities  and  broad 
angular  extent  are  more  probable  in  years  of  low  or  declining 
sunspot  numbers. 
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SOLAR  MAGNETIC  SECTORS 
OCT.  64- JAN.  65 


Figure  6.   Solar  magnetic  sectors  1964-65. 
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The  sharpness  of  the  sector  boundary  has  been  interpreted 
to  mean  that  the  magnetic  fields  within  a  sector  are  actually 
open  lines  rather  than  part  of  a  closed  loop  structure.   The 
field  itself  is  regarded  as  being  an  extension  of  closed  loop 
features  in  the  lower  corona.  As  field  lines  from  these  lower 
structures  extend  into  the  2  to  20  solar  radii  region  they  are 
swept  away  by  the  wind  as  it  is  rapidly  accelerated  outward. 
Newkirk  (1971)  has  identified  a  class  of  linear  north-south 
loop  patterns  in  computer  extrapolations  of  the  photospheric 
fields  into  the  corona.   These  features,  dubbed  "coronal 
arcades"  for  their  similarity  in  appearance  to  covered  walk- 
ways, are  thought  to  be  the  origin  of  the  sector  boundaries 
(Newkirk,  1972). 

The  regions  from  which  the  high  speed  streams  are  emitted 
may  be  "coronal  holes"  characterized  by  abnormally  low  densities 
and  temperature  (Krieger  et_al . ,  1973).   In  these  regions  the 
diverging  magnetic  field  is  believed  to  play  the  role  of  a 
de  Laval  nozzle  in  accelerating  the  wind  to  abnormally  high 
speed  (Adams  and  Sturrock,  1975). 

Velocity  Stream  Interactions 

Snowplow  Effect 

During  periods  of  stable  sector  structure  an  observer 
in  interplanetary  space  near  1  AU  sees  a  periodic  fluctuation  in 
the  solar  wind  velocity.   This  variation  will  have  a  frequency 
of  twice  the  number  of  sectors  since  each  sector  has  a  high 
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speed  and  a  slow  speed  component.   Figures  7  and  8  show  that 
the  day  of  the  sector  boundary  passage  is  at  or  near  a  velocity 
minimum,  followed  1  to  3  days  later  by  a  maximum—the  beginning 
of  the  high  speed  stream — which  is  sustained  for  a  few  days 
then  diminishes  to  a  minimum  at  the  next  boundary  passage. 

As  the  high  and  low  speed  components  propagate  through 
space  the  fast  stream  from  the  beginning  of  one  sector 
gradually  overtakes  the  slow  stream  from  the  end  of  the  pre- 
ceeding  sector.   Each  stream  consists  of  fully  ionized  plasma 
in  a  magnetic  field  and  consequently  one  stream  may  not  pass 
through  another,  though  Papadopoulos  (1973)  has  shown  that 
interpenetration  of  two  streams  may  occur  up  to  about  one  ioh 
gyroradius.   Eventually  there  is  a  "collision"  between  the  two 
streams  which  compresses  both  streams  adiabatically,  raising 
the  temperature,  magnetic  field  and  particle  density  (Burlaga 
et_al . ,  1971).   The  plots  in  Figure  7  show  the  density  build- 
up in  the  trough  just  before  the  velocity  wave.   This  process 
was  aptly  likened  to  the  buildup  of  snow  before  the  blade  of 
a  snowplow  by  Sonett  in  a  comment  on  a  paper  by  Davis  (1972). 

The  collision  process  is  inelastic  as  it  seems  the 
material  from  both  streams  becomes  part  of  a  joint  entity — the 
building  density  wave — share  momentum,  and  continue  at  a  speed 
intermediate  to  those  of  the  original  streams  with  the  sector 
boundary  trapped  in  the  wave.   In  time  the  original  velocity 
waves  will  be  damped  out  by  these  interactions  lea/ing  in- 
stead a  series  of  density  waves  moving  at  the  average  wind 
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Figure  7.   Solar  wind  velocity  and  particle  density 
versus  days  from  magnetic  sector  boundary, 
(After  Wilcox  and  ^7ess,  1965.) 
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velocity.   At  this  point  the  magnetic  fields  will  become 
perpendicular  to  the  radial  flow  of  the  density  waves. 
Siscoe  and  Finley  (1970)  have  applied  velocity  perturbations 
to  a  non-spherical ly  symmetric  steady  flow  wind  model  with 
essentially  this  result,  while  the  comparative  study  of  the 
Pioneer  X  and  XI  data  by  Collard  and  Wolfe  (1974)  shows  the 
strong  damping  effect  of  the  velocity  waves  near  4.5  AU  in 
1973. 

There  is  another  effect  of  the  stream  interaction  pro- 
cess which  may  be  a  factor  in  stimulating  Jovian  emissions — 
shock  waves.   The  bulk  wind  velocity  is  superalfvenic  and 
supermagnetosonic.   In  the  region  between  the  earth  and  Jupiter 
where  much  of  the  stream  interaction  occurs,  the  fast  magneto- 
sonic  wave  velocity  is  in  the  30  to  70  km/s  range.   It  will 
be  seen  shortly  that  the  velocity  difference  between  colliding 
streams  is  often  greater  than  the  MHD  wave  velocity.   As  a 
result  a  by-product  of  such  a  collision  is  the  production  of 
two  fast  mode  shocks,  one  moving  outward  and  other  other  moving 
inward  relative  to  the  collision  but  carried  outward  by  the 
bulk  v/ind  flow.   Experimental  evidence  of  this  effect  is  seen 
in  geomagnetic  disturbances  referred  to  as  sudden  impulse  (SI) 
pairs  (see  for  instance  Razdan  et  al . ,  1965;  Dessler  and  Fejer, 
1973;  and  Beard,  1965). 

Stream  Interaction  Models 

Since  one  of  the  objectives  of  this  research  is  to  com- 
pare the  arrival  of  the  density  wave  peaks  at  Jupiter  with 
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observed  Jovian  decametric  emissions,  a  model  describing  the 
process  of  the  formation  of  the  density  waves  is  required. 
Velocity  wave  data  at  1  AU  will  be  used  as  input  to  such  a 
model  to  predict  the  density  waves'  arrival  at  Jupiter.   In 
order  to  simplify  this  task  as  much  as  possible  the  two 
apparitions  to  be  studied  were  chosen  to  be  periods  with 
stable  two  sector  suns  as  shown  in  Figure  9.   The  benefit 
of  this  choice  is  two-fold.   First,  it  reduces  the  number  of 
waves  considered  to  a  minimum — two  per  solar  revolution. 
Second,  it  reduces  the  amount  of  interaction  to  a  minimum 
since  the  fast  and  slow  streams  start  out  as  widely  separated 
as  possible.   Figure  8  shows  the  velocity  wave  structure  at 
1  AU  during  these  periods. 

In  devising  a  suitable  model  for  the  stream  interaction 
it  must  be  remembered  that  only  large  scale  effects  are  of 
interest.   The  microscopic  details  of  the  interaction  appear 
to  be  very  complex  and  they  will  be  ignored  as  much  as  possible 
in  favor  of  predicting  the  macroscopic  development  of  the 
density  wave.   Intuitively,  the  smaller  the  amount  of  inter- 
action that  takes  place  before  the  wave-Jupiter  encounter  the 
less  important  the  microscopic  details  are. 

Input  Data 

Before  proceeding  with  the  development  of  the  models  a 
discussion  of  the  treatment  of  the  experimental  data  is  in 
order.   The  reduction  procedure  for  the  October  1967  -  May  1968 
apparition  and  that  of  the  April  1974  -  December  1974  apparitions 
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Figure  9.   Solar  magnetic  sectors  1967-68  and  1974, 
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differed  in  some  respects.   The  ultimate  goal,  however,  was 
the  same. 

Figure  10  shows  the  relative  positions  of  the  earth,  sun, 
and  Jupiter  for  some  arbitrary  time  before  opposition.   Since 
the  velocity  of  the  streams  is  radial,  the  data  of  interest  are 
the  velocities  as  a  function  of  time;  at  1  AU  on  the  line 
connecting  the  sun  with  Jupiter,   Since  in  general  no  obser- 
vations are  continuously  made  at  that  point,  it  is  necessary 
to  either  interpolate  or  extrapolate  what  data  are  available 
to  that  point. 

Despite  the  radial  motion  of  the  stream,  the  "garden 
hose"  effect  causes  the  position  of  the  point  of  intersection 
of  the  stream  on  the  1  AU  orbital  circle  to  rotate  with  the 
sun.  Consequently,  any  projection  of  measurements  made  off 
the  sun- Jupiter  radial  to  the  radial  involves  not  only  a 
projection  in  space  but  in  time  as  well.  The  intersection 
point  moves  in  the  posigrade  sense  about  13 •*  per  day.   Conse- 
quently, if  a  measurement  is  made  at  a  point  13°  behind  the 
radial  line,  that  value  becomes  an  estimator  of  the  velocity 
on  the  radial  the  following  day.   If  interpolation  between  two 
satellites  which  straddle  the  sun-Jupiter  radial  is  employed, 
one  must  assume  that  any  changes  that  occur  are  sufficiently 
gradual  that  no  inflections  lie  between  the  two  observations. 
If  extrapolation  from  a  single  point  is  used,  one  simply 
assumes  that  the  velocity  is  stable  over  the  period  of  time 
involved. 
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JUPITER 


Figure  10.      Sun,    earth  and  Jupiter  geometry. 
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Gosling  and  Bame  (1972)  have  shown  that  if  one  makes  a 
prudent  choice  of  apparition  the  stream  structure  one  full 
rotation  away — 27  days — may  autocorrelate  with  the  current 
structure  r    >0.1.      It  must  be  added  that  such  a  high 
correlation  is  the  exception  rather  than  the  rule  and  thus 
the  operative  phase  is  "prudent  choice."   The  validity  of 
assumptions  of  stability  is  strengthened  if  the  27-day  auto- 
correlation coefficient  is  reasonably  high  and  one  interpo- 
lates or  extrapolates  over  periods  significantly  less  than 
27  days. 

The  solar  wind  data  are  in  the  format  of  a  single 
average  value  per  day.   It  is  also  necessary  to  point  out 
that  since  the  Jovicentric  solar  sector  rotation  period  is 
the  period  of  real  interest  in  comparing  solar  and  Jovian 
events,  the  solar  rotation  periods  referred  to  in  the  balance 
of  this  study  will  be  the  solar  magnetic  sectors'  Jovicentric 
synodic  rotation  period.   This  period  is  nominally  25.7 
terrestrial  solar  days. 

1967-68  apparition.   The  velocity  information  during 
this  time  was  that  from  Pioneers  VI  and  VII.  These  satellites 
were  in  solar  orbit  very  near  1  AU.   They  were  separated  by 
about  155°  in  the  ecliptic  plane  such  that  they  straddled  the 
sun- Jupiter  radial  throughout  the  apparition.   Thus,  one 
satellite  or  the  other  was  never  more  than  5  rotation  days 
from  the  radial.   The  table  of  values  for  the  velocity  on  the 
radial  was  obtained  by  a  linear  interpolation  in  space  of  the 
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respective  estimators  for  each  day  at  the  radial  line.  This 
is  to  say,  the  value  of  the  velocity  at  each  satellite  was 
projected  in  time  to  the  radial.   The  value  chosen  for  the 
velocity  on  the  radial  for  a  given  day  was  the  average  of  the 
corresponding  measurements  from  each  satellite  linearly 
weighted  in  favor  of  the  nearer  satellite.   The  autocorrelation 
of  the  1967-68  period  is  shown  in  Figure  11  to  have  a  value  of 
0.44.   This  is  regarded  as  quite  satisfactory  for  the 
"interpolation. " 

1974  apparition.   The  solar  wind  velocity  data  for  this 
period  of  time  were  provided  by  the  IMP  7  and  8  satellites 
which  were  in  highly  eccentric  earth  orbits.   Unfortunately, 
this  situation  was  not  as  satisfactory  as  the  1957-68  period 
since  interpolation  was  not  possible.  The  only  alternative 
was  to  extrapolate  the  single  daily  average  to  the  sun-Jupiter 
radial  line.   Figure  12  shows  the  autocorrelation  function  for 
the  1974  period.   The  peak  in  the  function  after  one  solar 
rotation  is  0.81.   However,  the  1974  data  picture  may  be 
overly  optimistic.   At  this  writing  the  only  available  velocity 
data  for  the  1974  are  for  the  81-day  period  from  June  25  to 
September  14.   These  data  show  remarkable  repeatability  in 
the  velocity  structure.   However,  this  is  an  inadequate 
length  of  time  to  permit  proper  cross-correlation  with  Jovian 
data.   Svalgaard  (1976),  on  the  other  hand,  has  shown  that 
the  magnetic  sector  structure  was  quite  stable  throughout  most 
of  the  year.   In  order  to  satisfy  the  requirements  of  the 
cross-correlation  algorithm  the  June-September  data  have  been 
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Figure  11.   Solar  wind  velocity  autocorrelation  function, 
1  AU  1967-68. 
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Figure  12.   Solar  wind  velocity  autocorrelation  function, 
1  AU  1974. 
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used  to  simulate  the  81-day  periods  immediately  preceding 
and  following  the  true  June-September  period.   No  doubt  this 
procedure  has  introduced  some  errors  which  may  be  eliminated 
when  more  complete  solar  wind  data  become  available. 

Magnetic-Field-Free  Model 

The  simplest  model  for  the  stream  interaction  is  to 
assume  there  is  no  interaction.   In  the  absence  of  a  magnetic 
field  the  solar  wind  is  simply  a  collisionless  plasma.   Conse- 
quently, no  interaction  between  streams  occurs.  Nevertheless, 
this  model  predicts  the  buildup  of  a  large  scale  (length 
>0.2  AU)  density  wave  as  the  fast  stream  approaches  the  slow 
stream  leaving  a  rarefaction  behind.   This  development  of  a 
density  wave  from  a  velocity  wave  is  analogous  to  the  "bunching" 
process  in  a  klystron. 

The  field-free  model  is  applicable  in  situations  where 
the  initial  velocity  spread  between  adjacent  streams  and  their 
initial  separation  in  space  are  such  that  fast  streams  do  not 
often  have  the  opportunity  to  pass  completely  through  slow 
streams  before  reaching  Jupiter.   Since  the  model  permits  no 
stream  interaction,  at  some  point  in  space  all  of  the  fast 
streams  will  have  passed  through  the  slow  streams  and  the 
density  wave  will  decay.   In  nature,  this  decay  probably  does 
not  occur?  consequently,  a  basic  test  of  the  applicability  of 
the  field-free  model  will  be  to  ascertain  whether  the  model, 
given  the  input  velocities  for  a  given  apparition,  predicts 
such  a  decay  before  the  encounter  with  Jupiter,  (i.e.,  before 


62 

about  5.2  AU).   If  the  decay  does  not  occur,  the  model  can  be 
used.   The  precise  manner  of  performing  this  test  will  be 
detailed  shortly. 

A  FORTRAN  program  designated  "Snoplow"   (see  Appendix) 
was  developed  to  generate  the  modeled  values  for  wind  density 
and  flux  as  a  function  of  time  and  space.   Starting  with  the 
table  of  solar  wind  velocities  at  1  AU  on  the  sun-Jupiter 
radial,  this  program  extrapolates  these  values  back  in  time 
and  space  to  the  sun.   Then  the  model  releases  groups  of  par- 
ticles from  the  sun  with  constant  arbitrary  density  at  six- 
hour  intervals.   The  four  groups  for  a  given  day  are  assigned 
the  wind  velocity  extrapolated  from  the  earth  for  that  day. 
The  progress  of  each  group  is  tracked  for  forty  days  to  insure 
that  even  the  slowest  reaches  Jupiter  before  they  are  discarded. 

Once  a  day  the  groups  are  examined  to  determine  the  den- 
sity and  flux  in  space  at  increments  of  0.2  AU.   Since  part  of 
this  program's  objective  is  to  permit  a  view  of  the  wave  build- 
ing process,  the  decrease  in  density  caused  by  the  expansion  of 
the  wind  (i.e.,  the  1/r  dependence)  is  deliberately  omitted 
in  this  calculation.   As  a  result,  the  "density"  and  "flux" 
calculated  are  really  the  product  of  the  true  density  or  flux, 
the  square  of  the  radial  distance,  and  some  scale  factor. 
Clearly,  for  any  given  radial  distance  this  product  is  directly 
proportional  to  the  true  density  at  that  point.   Consequently, 
the  term  "density"  will  frequently  be  used  to  denote  the  product. 

The  net  result  of  Snoplow  is  two  two-dimensional  arrays, 
one  representing  the  density,  and  the  other  the  flux,  as 
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functions  of  time  and  space  out  to  6  AU.,  Figure  13  is  an 
example  of  such  an  array.   Prom  this  one  can  select  a  row 
corresponding  to  a  given  day  and  see  the  density  or  flux  as  a 
function  of  space,  or  one  can  select  a  column  corresponding 
to  a  radial  distance  and  see  the  density  or  flux  as  a  function 
of  time. 

Growth  test.   In  order  to  test  the  applicability  of  the 
field-free  model,  a  simple  auxiliary  program  examines  the 
density  array.   At  each  increment  of  radial  distance  the  pro- 
gram calculates  the  average  and  standard  deviation  of  the  den- 
sity across  the  apparition.   The  program  then  plots  the  standard 
deviation,  expressed  in  units  of  the  average  density,  for  each 
increment  in  radial  distance.   This  function,  then,  is  merely 
the  normalized  RMS  wave  amplitude  as  a  function  of  distance. 
If  a  density  wave  building  process  is  occurring,  the  amplitude 
will  monotonically  increase  with  distance.   Conversely,  if  a 
decay  is  occurring  the  amplitude  will  monotonically  decrease 
with  distance.   The  derivative  of  the  amplitude  function  may 
also  be  used  as  an  index  of  growth. 

Simple-Field  Model 

A  second  model  has  been  developed  to  handle  the  case 
where  the  growth  test  contraindicates  the  field-free  model. 
This  model  implicitly  assumes  the  presence  of  the  interplanetary 
magnetic  field,  but  does  not  attempt  to  account  for  the  fine 
detail  of  the  stream  interaction.   The  model  invokes  three 
assumptions: 
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1.  Groups  of  particles  cannot  pass  through  one 
another  (B  7^  0) , 

2.  The  stream  collisions  are  completely  inelastic, 

3.  Momentum  is  conserved. 

These  assumptions  have  been  incorporated  as  a  user  option  in 
the  program  Snoplow  outlined  in  the  previous  section. 

Since  Snoplow  treats  each  four  hour  group  of  particles 
as  independent  entities  rather  than  as  members  of  a  particular 
stream,  the  assumptions  are  applied  at  the  group  level. 
Consequently,  this  model  will  permit  the  wave  to  build  up  as 
an  infinitely  thin  sheet  of  mass.   This  is  physically  un- 
realistic because  both  plasma  and  magnetic  pressure  will 
limit  the  compression  of  the  density  front.   The  effect  is  off- 
set somewhat  by  the  coarse  scale  of  the  radial  distance  used 
in  the  density  and  flux  arrays.   This  insures  that  the  density 
wave  will  never  appear  smaller  than  0.2  AU. 

There  is  another,  probably  related,  effect  that  the  model 
does  not  account  for.   If  one  refers  back  to  Figure  7,  Wilcox 
and  Ness  show  that  the  peak  of  the  density  wave  front  actually 
precedes  the  peak  of  the  velocity  front  by  about  one  day  at  1  AU. 
The  forward  one-half  amplitude  point  of  the  density  front  pre- 
cedes the  corresponding  point  on  the  velocity  front  by  two  days. 
However,  the  model  will  place  the  density  front  at  the  collision 
point — the  leading  edge  of  the  velocity  wave.   There  are  no 
data  available  on  the  phase  relationship  between  velocity  and 
density  fronts  at  5.2  AU.   However,  if  the  1  AU  characteristic 
persists  out  to  Jupiter,  the  model  will  consistently 
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overestimate  the  arrival  date  of  the  density  fronts  at  Jupiter 
by  as  much  as  2  days. 

Model  Selection 

The  first  step  in  the  selection  of  the  appropriate  model 
was  to  invoke  the  field-free  model  for  both  apparitions.   The 
growth  test  routine  was  then  run  to  plot  the  density  wave 
growth  for  each  apparition. 

Figure  14  shows  the  amplitude  plot  for  1967-68.   It 
indicates  that  the  RMS  density  wave  amplitude  at  1  AU  was  0.47, 
which  is  quite  consistent  with  values  seen  at  the  earth.   The 
plot  also  shows  that  the  principal  evolution  of  the  density 
fronts  took  place  beyond  1  AU  and  that  the  growth  arrested  at 
about  3.4  AU.   There  is  a  plateau  which  is  essentially  con- 
stant through  5.8  AU  and  then  the  decay  process  apparently 
begins  in  earnest.   It  was  decided  that,  since  there  was  little 
decay  at  5.4  AU  (the  location  of  Jupiter  in  1967-68),  the  field- 
free  model  was  an  acceptable  choice  for  this  apparition. 

The  field-free  amplitude  plot  for  the  1974  apparition 
is  shown  in  Figure  15.   In  contrast  to  1967-68,  it  shows  a 
much  higher  growth  rate  with  the  amplitude  peaking  at  2.2  AU 
followed  by  an  immediate  decay.   This  characteristic  clearly 
shows  that  the  field-free  model  is  not  satisfactory  for  the 
1974  period. 

The  accelerated  growth  of  the  1974  density  waves,  as 
compared  to  the  1967-68  waves,  stems  directly  from  the  more 
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Figure  14.   Field-free  model  RMS  density  wave  amplitude, 
1967-68. 
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Figure  15.   Field-free  model  RMS  density  wave 

amplitude,  1974.   The  inflection  at 
2.2  AU  shows  that  the  model  fails 
beyond  that  point. 
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profound  peak-to-peak  amplitude  of  the  1974  velocity  wave. 
Referring  back  to  Figure  8  it  is  clear  that  value  of  the  peak- 
to-peak  velocity  averages  from  cycle  to  cycle  was  130  km/s  in 
1967-68,  but  the  value  in  1974  was  450  km/s.   Consequently, 
the  fast  streams  overtook  the  slow  ones  much  faster  in  1974. 

The  simple-field  model  was  then  applied  to  the  1974  data 
and  a  new  amplitude  plot  obtained.   The  latter  is  displayed  in 
Figure  16.   The  1  AU  value,  0.77,  is  higher  than  that  usually 
observed  but  not  substantially  different  from  the  1  AU  field- 
free  value.   The  rapid  growth  evident  in  the  field-free  model 
is  also  evident  here  in  the  simple-field  model.   The  1  AU 
amplitude  value  is  probably  reasonable  in  that  light,  though 
1974  observational  data  are  not  yet  available. 

The  growth  process  is  seen  to  diminish  as  it  approaches 
the  Jovian  orbit,,  suggesting  that  the  interaction  process  is 
nearing  completion.  This  is  an  expected  result  for  this  model. 
Once  a  fast  stream  has  completely  accumulated  the  preceding 
slow  stream,  the  only  further  interaction  which  can  take  place 
will  be  between  the  density  fronts  themselves.   Since  they  are 
widely  separated  and  moving  at  about  the  same  speed,  such 
collisions  would  be  very  infrequent. 

In  comparing  the  two  models  out  to  2.2  AU  before  the 
field-free  model  breaks  down,  the  two  both  provide  similar 
amplitude  values.   As  expected,  even  in  this  region,  the 
simple-field  model  shows  a  somewhat  higher  growth  rate.   The 
accumulated  amplitude  exceeds  that  of  the  field-free  model  by 
34%  at  2.2  AU.   Nevertheless,  it  is  comforting  to  see  that 
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Figure  16.   Simple-field  model  RMS  density  wave  amplitude 
1974. 
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where  they  are  both  valid,  these  two  basically  very  different 
models  do  not  perform  differently  in  any  essential  way.   It  is 
clear. though  that,  of  the  two,  the  simple-field  model  would  have 
to  be  used  for  the  197  4  data. 


CHAPTER  4 
SOLAR-JUPITER  CORRELATION 

Previous  Investigations 

The  matter  of  the  relationship  of  solar  and  geomagnetic 
events  was  discussed  at  some  length  in  the  preceding  chapter. 
By  the  time  Burke  and  Franklin  published  their  discovery  of 
the  Jovian  decametric  emission  in  1955>  the  solar-terrestrial 
correlations  were  already  well  recognized,  even  though  their 
precise  nature  was  not.   Since  solar  events — particularly 
solar  flares— were  seen  to  produce  such  dramatic  effects  as 
visible  and  radio  aurora,  auroral  absorptions  (PCAs)  and  pro- 
found geomagnetic  disturbances,  it  was  reasonable  for  early 
students  of  the  Jovian  radiation  to  consider  the  possibility 
of  solar-Jovian  relationships.  This  natural  flow  of  thought 
no  doubt  was  intensified  by  the  timing  of  the  discovery  just 
prior  to  the  International  Geophysical  Year,  the  peak  of 
solar  Cycle  19,  and  Parker's  theoretical  insights  on  the 
nature  of  the  solar  wind. 

Several  processes  were  proposed  to  explain  the  emissions 
such  as  Jovian  lightning,  plasma  oscillations  excited  by 
shocks  caused  by  Jovian  volcanoes,  and  atmospheric  chemical 
activity.   In  1958,  Carr  proposed  that  the  emissions  were 
triggered  by  energetic  particles  of  solar  origin  trapped  by 
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the  Jovian  magneto sphere.   At  about  the  same  time  Kraus  (1958) 
noticed  two  intense  Jovian  storms  following  two  large  solar 
flares  and  made  the  same  suggestion.   The  possibility  of  such 
a  connection  was  actively  pursued  by  a  number  of  investigators. 

What  early  observers  did  not  know  was  that  there  was  a 
very  powerful  and  periodic  exciter,  lo,  which  was  responsible 
for  much  of  the  activity  they  were  seeing.   The  influence  of 
the  lo  effect  was  clearly  evident  in  these  early  observations. 
Carr  et  al.  (1958)  reported  that  they  had  found  a  very  definite 
8-day  period  in  their  observed  activity.   (The  8-day  period — 7.71 
days  actually — is  nothing  but  the  alias  of  the  second  harmonic 
of  the  lo  revolution  period  with  the  24-hour  sampling  period 
imposed  by  the  rotation  of  the  earth.)   They  then  suggested 
that  the  period  was  due  to  tidal  action  induced  by  one  of 
Jupiter's  satellites.   Unfortunately  this  possibility  was 
not  explored  more  exhaustively  at  that  time.   The  discovery  of 
the  lo  effect  would  wait  for  Bigg  in  1964. 

The  presence  of  the  lo  modulation  would  cloud  efforts 
to  detect  solar  correlations.   There  was  yet  another  pitfall. 
Until  Snyder  et  al.  (1953)  published  their  Mariner  2  obser- 
vations, the  velocity  range  and  the  stream  structure  of  the 
solar  wind  were  not  known.   Consequently,  there  was  great 
uncertainty  as  to  the  appropriate  delay  time  between,  say, 
an  M-region  geomagnetic  storm  and  an  expected  Jovian  event. 
The  only  events  where  the  timing  might  be  inferred  were 
solar  flares  followed  by  SCs.   In  this  case  the  velocity  of 
the  wind  disturbance  was  measurable,  at  least  between  the 
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sun  and  the  earth.   Since  these  early  investigations  took 
place  near  the  peak  of  Cycle  19,  flares  were  in  relative 
abundance . 

Pursuing  the  search  for  a  solar  correlation  Carr 
et_al.  (1960)  presented  evidence  to  suggest  that  there  was 
an  enhancement  of  Jovian  activity  8  days  after  geomagnetic 
disturbances.   Warwick  (1960)  examined  solar  radio  continuum 
events  and  reported  increased  activity,  particularly  1  to  2 
days  after  the  cessation  of  the  continuum.   More  evidence 
for  an  8-day  delay  after  geomagnetic  events  was  introduced  by 
Carr  et  al.  (1961). 

On  the  other  hand.  Six  (1962)  examined  the  8-day  delay 
and  found  no  clear  correlation.   Lebo  (1964)  studied  the 
delay  from  flares  and  geomagnetic  storms  but  extended  the 
maximum  length  of  delay  considered  from  15  to  35  days.   This 
was  an  important  step  because  the  just  published  Mariner  2 
data  indicated  that  the  average  solar  wind  velocity  was  in 
the  500  km/s  range  and  the  instantaneous  values  were  often 
smaller.  Thus,  the  average  delay  between  the  earth  and 
Jupiter  was  15  days.   In  doing  so  he  found  indications  of 
peaks  in  Jovian  activity  with  delays  in  8  to  9  days  and  15 
to  21  days.   Lebo  also  suggested  that  future  studies  be 
expanded  to  cover  M-region  events  as  well. 

Armed  with  recognition  of  the  lo  effect  Sastry  (1968) 
examined  the  question  anew.   He  recognized  correctly  the 
complexities  imposed  by  the  "artificial"  periodicities  in 
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the  Jovian  data  on  efforts  to  find  time  dependent  effects, 
such  as  time-of-f light  delays.   However,  he  did  not  attempt 
to  remove  or  compensate  for  the  effects.   Sastry  reexamined 
Jovian  radio  activity  from  1958  through  1963  and  concluded 
that  the  probability  of  emission  is  higher  after  geomagnetic 
activity  than  before.   He  cautioned  that  the  correlation  might 
be  an  artifact  produced  by  the  periodicities. 

A  paper  by  Conseil  et  al.  (1971)  shows  a  tanta- 
lizing but  invalid  relationship  between  the  sign  of  the 
derivative  of  the  solar  wind  velocity  and  the  definition  of 
"lo-related"  activity.   They  have  considered  Jovian  activity 
and  solar  wind  velocity  measurements  of  Pioneers  VI  and  VII  in 
orbit  at  1  AU  during  1967-68 — the  same  period  as  one  of  the 
two  examined  in  this  study.   Unfortunately,  they  explicitly 
required  that,  "propagation  effects  [in  the  solar  wind]  are 
not  very  important."   It  must  be  clear  that  the  extent  of 
these  interactions  is  all- important.  However,  the  effect 
they  observed  may,  in  fact,  be  real,  but  for  the  wrong  reason. 

During  the  1967-68  apparition,  it  has  been  argued  here 
that  the  difference  between  the  velocities  of  high  and  low 
speed  streams  was  such  that  the  field-free  or  non-interaction 
model  adequately  describes  the  evolution  of  the  density  wave 
out  to  Jupiter.   Conseil 's  assumption  of  the  unimportance  of 
the  interaction  constitutes  the  invocation  of  the  field-free 
model.   Their  error  was  in  failing  to  consider  the  density 
wave. 
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Using  the  satellites'  wind  velocity  data,  Conseil  and 
his  colleagues  developed  a  table  of  solar  wind  velocities 
versus  dates  at  Jupiter.   They  then  examined  the  Jovian 
decametric  data  for  correlations  with  the  velocity  table. 
What  they  reported  was  that  the  "favored"  range  of  lo  phase 
(ISC)  for  source  A  and  source  C  storms  shifted  toward  lower 
angles  on  days  when  the  wind  velocity  was  increasing  and  toward 
higher  angles  when  the  velocity  was  decreasing.   This  effect 
is  shown  graphically  in  Figure  17. 

However,  it  seems  very  unlikely  that  the  velocity  deriva- 
tive at  1  AU  would  be  maintained  at  Jupiter.   The  derivative 
of  the  velocity  is  necessarily  far  more  sensitive  to  stream  ' 
interactions  than  is  the  velocity  itself.  But,  it  is  known 
that  a  positive  velocity  derivative  at  1  AU  corresponds  to  the 
leading  edge  of  the  fast  stream  which,  in  turn,  is  coincident 
with  the  peak  of  the  density  wave.   In  a  similar  manner,  a 
negative  velocity  derivative  corresponds  to  the  trough  of 
the  density  wave.  As  a  result,  one  may  interpret  their  re- 
sult as  actually  being  a  correlation  between  the  "favored" 
range  of  lo  and  the  density  wave,  not  the  velocity  wave  as 
they  propose. 

Kovalenko  (1971)  examined  the  relationship  between 
Jovian  emissions  and  interplanetary  shocks  caused  by  flares 
and  stream  collisions.   Using  the  "Forbush  decreases"  in 
cosmic-ray  intensity  as  an  indicator  of  a  shock  passage  at 
the  earth,  he  used  an  extrapolation  procedure  identical  to 
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Figure  17.   Decametric  storms  versus  ISC  for  sources 
A  and  C,  1967-68.   (After  Conseil  et  al . , 
1971.) 
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the  one  described  in  the  preceding  chapter  to  place  the  shock 
on  the  sun-Jupiter  radial.   He  concludes  that  for  the  period 
1961  to  1964  there  was  a  one-to-one  correspondence  between 
geomagnetic  events  and  Jovian  activity. 

The  principal  shortcoming  of  Kovalenko's  work  stems  from 
the  fact  that  he  had  no  solar  wind  velocity  data.   As  a  conse- 
quence he  was  forced  to  infer  the  velocity  at  1  AU  using  an 
algorithm  which  related  the  wind  velocity  to  the  magnitude 
of  the  shock  and  its  accompanying  geomagnetic  storm.   He 
also  failed  to  make  any  separation  of  lo  and  non-Io  related 
events.   Consequently,  the  7.7  day  lo  cycle  may  well  have  pro- 
duced the  "27-day"  cycle  he  reports  (4  x  7.7  =  28.3).   Had  the 
cycle  been  due  to  a  solar  influence  the  period  would  have  been 
the  Jovicentric  synodic  period,  25.7  days. 

In  1972  in  a  precursor  of  the  present  study,  the  author 
presented  an  apparent  correlation  between  the  non-lo  related 
component  of  source  A  and  the  projected  passage  of  sector 
boundaries  at  Jupiter  in  the  1964-65  apparition  (Kennedy, 
et_al. ,  1973) .    In  that  work  the  occurrence  of  sector  boundaries 
at  the  earth,  projected  to  the  sun-Jupiter  radial,  and  Jovian 
radio  data  were  each  subjected  to  a  superposed  epoch  analysis 
using  the  Jovicentric  25.7  day  period  as  the  base  period. 
These  two  plots  were  then  shifted  relative  to  each  other  to 
see  if  there  was  a  time  shift  corresponding  to  the  time-of- 
f light  delay  in  the  sector  boundary  arrival  at  Jupiter.   The 
best  match  between  the  two  plots  occurred  at  a  shift  of  13.5 
days.   This  delay  indicates  an  average  solar  wind  velocity  of 
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514  km/s  during  the  period  of  the  study  (the  average  for  the 
entire  year  1964  was  given  by  Gosling  et  al.  (1971)  as  4G5  km/s) 
Figure  18  shows  this  comparison. 

This  approach  was  regarded  as  simplistic  in  that  it 
demanded  that  the  sector  boundary  propagate  at  the  average 
wind  velocity.   At  the  time  it  seemed  likely  that  the  velocity 
of  the  sector  would  be  modified  by  the  stream  interactions — 
which  had  not  been  accounted  for.   In  retrospect,  however, 
the  procedure  may  have  been  better  than  it  appeared.  The 
sector  boundary  does  become  trapped  in  the  density  wave  which 
is  traveling  at  an  intermediate  velocity  due  to  the  stream 
collision  process.   Furthermore,  as  Collard  and  Wolfe  (1974) 
have  shown,  the  peak-to-peak  velocity  dampens  and  converges 
on  the  average  velocity  as  the  stream  moves  outward.   During 
1964-65  the  sun  had  four  sectors,  which  means  the  distance 
between  fast  and  slow  streams  was  half  that  of  a  two  sector 
sun.   Though  velocity  data  are  not  available,  it  is  safe  to 
assume  that  the  amount  of  interaction  damping  was  much  higher 
at  1  AU  due  to  the  four  sectors  than  during  the  two-sector 
periods  of  the  present  study.   The  result  would  be  to  further 
insure  that  the  density  wave  at  the  earth  was  already  at  or 
near  the  average  wind  velocity. 

Barrow  (1972)  has  examined  Jovian  data  from  1961  to 
1968.   He  noted  that  due  to  the  earth-sun-Jupiter  geometry 
the  time-of-f light  delay  is  a  function  of  the  position  of  the 
earth  relative  to  Jupiter.   For  example,  after  opposition  par- 
ticles on  the  sun-Jupiter  radial  are  already  past  1  AU  by  the 
time  the  solar  rotation  presents  that  stream  to  the  earth. 
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Figure  18.   Superposed  epoch  analyses  of  sector 

boundaries  at  1  AU  and  Jovian  activity 
in  1964.  The  boundary  plot  is  delayed 
13.5  days  with  respect  to  the  activity 
plot.   (After  Kennedy  etal. ,  1973.) 


Unfortunately,  there  is  an  error  in  Barrow's  work. 
His  Figure  2  is  incorrect.   The  plot  should  have  been  made 
from  0   =-180°  to  ^  =  180°,  not  from  0°  to  360°.   This  error 
led  him  to  the  statement,  "...  immediately  before  opposi- 
tion .  .  .  values  may  be  such  that  Jupiter  can  be  active 
before  the  particles  reach  the  earth."  This  is  obviously 
incorrect  since,  before  opposition  the  particles  always 
arrive  at  the  earth  first.   This  is  compounded  by  the  fact 
that  he  does  not  differentiate  between  lo  and  non-Io  related 
events.   Since  he  has  selected  which  Jovian  storms  to  study 
in  a  way  that  emphasizes  major  storms,  much  of  the  data  he 
used  is  lo  related.   Barrow  reports  a  positive  correlation 
but  his  errors  cast  doubt  on  the  validity  of  the  result. 

Vladirairsky  and  Levitsky  (1976)  have  recently  reported 
a  correlation  between  Jovian  activity  and  solar  proton  flares. 
They  maintain  that  the  predominantly  more  dense  and  slower 
wind  found  in  the  solar  equatorial  plane  by  Hundhausen  etal. 
(1971)  forms  a  quiet  zone  which  shields  Jupiter  when  its 
heliocentric  latitude  is  between  2°  North  and  2°  South. 
They  report  that  when  Jupiter  is  outside  that  range  the  inci- 
dence of  Jovian  activity  1  day  after  a  proton  event  is  more 
than  twice  the  average  incidence. 

Solar  proton  events  are  caused  by  a  stream  of  very  high 
speed  particles  blown  away  from  the  sun  by  a  very  energetic  flare. 
The  velocity  of  these  streams  is  typically  10  to  5  x  10  km/s. 
Such  a  stream  is  different  from  the  main  body  of  flare  ejecta 
which  moves  slower  by  a  factor  of  10.   The  latter  are 
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responsible  for  SCs  while  the  proton  events  result  in  PCAs. 

The  sector  rotation  period  suggests  that  the  high 
speed  streams  are  emitted  at  solar  latitudes  of  about  20° 
North  and  South  rather  than  from  the  equator.   This  is  re- 
garded as  the  most  likely  explanation  of  Hundhausen ' s  finding 
of  higher  density  and  lower  speeds  in  the  equatorial  plane. 
Vladimirsky  and  Levitsky  believe  that  this  thin  band  acts  as 
an  imperfect  barrier  for  particles  that  try  to  cross  the 
equatorial  plane.   These  assertions  certainly  seem  plausible, 
as  does  the  proton-event- Jovian  relationship. 

Summary.   The  correlations  obtained  by  Conseil  etal . , 
Kennedy  et  al. ,  and  Vladimirsky  and  Levitsky  may  be  real. 
The  positive  results  of  solar-Jupiter  correlation  studies 
reported  by  other  workers  must  be  considered  less  than  con- 
vincing, if  only  because  of  the  uncertain  influence  of  the 
lo  related  events.   Sastry's  concern  that  the  periodic  lo 
effect  might  mimic  a  periodic  solar  relationship  is  well 
founded.   Most  of  the  early  studies  also  failed  to  take  the 
earth-sun-Jupiter  geometry  into  account.   Consequently, 
timing  errors  of  several  days  are  induced  in  the  projected 
arrival  of  M- region  events,  though  flare  correlations  would 
not  be  affected  by  this  problem. 

Analysis 

Chree  Analysis  and  Correlation 

Preceding  chapters  have  described  how  the  Jovian  deca- 
metric  data  were  filtered  and  segregated  to  provide  a  separate 


listing  of  each  source  at  each  frequency  for  both  the  lo  and 
non-Io  component.   Likewise  it  was  shown  how  solar  wind 
velocity  data  for  1  AU  were  used  to  generate  a  table  of  model 
particle  flux  densities  at  Jupiter.   Motivated  by  the  belief 
that  the  supposed  solar  influence  was  not  a  strong  effect, 
it  was  decided  to  use  a  special  case  of  Chree  analysis  to  mini- 
mize noise  in  the  data.   The  data  were  superposed  on  themselves 
with  a  period  of  about  25.7  days,  the  Jovicentric  synodic  solar 
rotation  period,  and  collected  in  25  channels  each  representing 
one  day. 

In  order  to  account  for  daily  sample  size  variations  the 
superposition  was  done  separately  for  the  observation  and 
activity  data.   In  this  fashion  the  zone  counts  for  observations 
for  each  channel  were  summed  over  the  long  apparition  while  the 
same  process  was  going  on  separately  with  the  activity  zone 
counts.  When  the  summations  were  complete,  the  ratio  of  the 
counts  for  each  pair  of  activity  and  associated  observation 
channels  represented  the  apparition-long  activity  index  for 
that  channel. 

The  table  of  solar  wind  flux  values  was  superposed  upon 
itself,  summed  and  then  averaged  to  give  the  solar  wind  infor- 
mation in  the  same  form  as  the  Jovian  data.   The  two  Chree 
analysis  tables  were  subjected  to  a  Banning  smoothing  proce- 
dure (see  Chapter  5)  and  then  shifted  with  respect  to  each 
other  in  one  day  (channel)  increments  with  the  cross-correla- 
tion coefficient  being  determined  at  each  shift.   A  shift  of 
zero  days  corresponds  to  the  projected  dates  of  arrival  of  the 
density  waves  at  Jupiter.   The  whole  process  of  Chree  analysis 
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and  cross-correlation  is  accomplished  by  a  computer  program 
designated  "Jupsol  V"  (See  Appendix) .   Figure  19  shows  an 
example  of  the  Chree  analysis  and  cross-correlation  functions. 
The  subject  of  the  development  of  confidence  limits  for 
the  Chree  analysis  tables  is  reserved  for  a  later  chapter. 
The  confidence  intervals  for  the  cross-correlation  procedure 
are  obtained  using  the  standard  "t  test"  (see  Hartnet,  1970). 
If  two  functions  are  both  randomly  distributed,  then  the 
function 

t  =  ^  '^  "  ~  2 

where  r  is  their  correlation  coefficient  and  n  is  the  number 
of  points  in  each  function,  has  a  Student  t  distribution  with 
n  -  2  degrees  of  freedom.  This  function  will  permit  testing 
a  null  hypothesis  that  the  functions  are  not  correlated.  The 
correlation  coefficient  limits  for  various  confidence  inter- 
vals for  the  25  point  cross-correlation  used  in  the  present 
study  are  given  in  Table  5. 


Table 

5 

Confidence 

Intervals 

Level 

Interval 

90% 

-0.337  <.  r  £  0.337 

95% 

-0.396  ±x   <_   0.396 

99% 

-0.505  ±x   <   0.505 
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Figure  19.   Superposed  epoch  analyses  of  Jovian  activity 
and  solar  wind  density  are  shown  in  the 
upper  two  plots.   The  lowest  plot  is  the 
plot  of  the  cross-correlation  function,  R, 
of  the  activity  and  density  functions. 
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The  Apparitions 

As  has  been  mentioned  the  two  Jovian  apparitions  to  be 
studied  here  are  those  of  1967-68  and  1974.   Table  6  has 
been  prepared  to  compare  and  contrast  the  various  conditions 
surrounding  these  measurement  periods.   It  is  noted  that  in 
many  ways  the  two  apparitions  are  antisymmetric.   The  dates 
of  opposition  are  about  6  months  apart,  the  radial  distances, 
Dg  and  sunspot  cycle  conditions  are  all  near  opposite  extremes, 
One  point  of  similarity  deserves  note.   The  heliocentric 
latitude  of  Jupiter  was  well  with  the  4°  quiet  zone  suggested 
by  Vladimir isky  and  Levitsky. 

Table  6 
Apparitional  Statistics 


Parameter 

1967-68 

1974 

Opposition 

23  February 

6  September 

Radial  Distance 

5.38  AU 

4.98  AU 

Synodic  Sector 
Period 

25.75d 

25.58d 

Heliocentric 
Latitide 

1° 

-1" 

Dg  Effect  on 
Activity 

Near  Minimum 

Near  Maximum 

Sunspot  Cycle 

Near  Maximum 

Near  Minimum 
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The  Cross-Correlations 

The  cross-correlation  functions  obtained  for  the  1967-68 
non-Io  component  of  source  A,  B,  and  C  are  presented  in 
Figures  20,  21,  and  22  respectively.   These  show  strong 
positive  correlations  at  or  near  zero  days  relative  shift  for 
sources  A  and  C  at  all  frequencies  except  for  source  A  at  18  MHz. 
The  plots  also  show  a  statistically  significant  anticorrelation 
for  source  B  at  18  MHz.   The  various  coefficients  and  the 
corresponding  time  shift  relative  to  the  predicted  wave-Jupiter 
encounter  are  summarized  in  Table  7.   It  should  be  noted  that 
day  zero  corresponds  to  the  encounter  day,  negative  days  are 
before  encounter  and  positive  days  are  after  encounter. 

The  cross-correlation  plots  for  the  1974  apparition  are 
shown  in  Figures  23  through  26.   They  show  a  similar  picture. 
Source  A  has  positive  peaks  at  26  and  27  MHz,  Source  B  anti- 
correlates  at  18,  22,  and  27  MHz,  while  source  C  has  positive 
correlation  peaks  at  18,  22,  and  27  MHz.   These  values  are  , 
summarized  in  Table  7. 

An  interesting  feature  of  the  1974  plots  is  the  fact 
that  the  features  of  interest  fall  not  in  the  -1  to  0  day 
range  but  they  occur  between  -4  and  -2  days. 

Conclusions 

During  both  apparitions  there  is  a  tendency  for  the 
extrema  to  precede  the  predicted  encounter  day.   However,  it 
has  been  noted  that  the  crests  of  the  density  waves  actually 
observed  at  1  AU  do  in  fact  precede  the  stream  collision  point 
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Table. 7 
Correlation  Extrema  Near  Encounter 


1967-68  Apparition 

Frequency 

Source 

(MHz) 

A               B 

C 

18 

-0.57*  (a  Od 

0.79  @ 

-Id 

22 

0.56  @  -Id          - 

0.87  @ 

Od 

27 

0.47  (a  -Id 

0.32  @ 

Od 

1974  Apparition 

18  -         -0.32*  @  -3.5d  0.54  @  -2d 

22  -         -0-51   @  -3.5d  0  g5  ^  _2^7d 

26  0.42  @  -3d  -  0.32  @  -2.7d 

27  0.40  @  3.5d    -0.48*  @  -3.5d        - 


*Dorainated  by  a  single  event, 
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by  about  one  day  and  the  half-amplitude  point  on  the  leading 
edge  of  the  wave  precedes  the  crest  by  still  another  day. 

Neither  of  the  two  models  took  this  fact  into  account, 
though  it  was  expected  that  the  problem  would  be  more  evident 
in  the  simple-field  model  due  to  the  very  thin  wave  fronts 
which  it  develops.   By  contrast,  the  field-free  model  permits 
some  diffusion  in  the  3.8  -  5.4  AU  region.   The  result  is  that 
the  wave  front  is  more  diffuse  and  the  leading  edge  of  the 
front  is  predominantly  composed  of  the  unimpeded  high  speed 
groups  which  reach  Jupiter  before  their  trapped  counterparts 
in  the  simple-field  model.   Consequently,  one  would  expect 
to  see  the  real  density  wave  interact  with  Jupiter  before  the 
day  predicted  by  either  model.  Moreover,  one  would  expect  to 
see  the  simple-field  model  predict  later  arrivals  than  the 
field-free  model.   Qualitatively  this  seems  to  be  the  case 
with  the  results  in  Table  7. 

Another  factor  which  may  account  for  the  early  arrival  of 
the  wave  is  the  acceleration  of  the  solar  wind  with  radial 
distance  predicted  by  many  quiet  wind  models.   This  accelera- 
tion was  not  observed  by  Mariner  2,  nor  was  it  seen  by 
Pioneers  X  and  XI.   However,  this  writer  questions  whether 
simple  time  averages  which  are  blind  to  the  density  varia- 
tions is  properly  descriptive  of  the  mass  transport  process. 

Finally,  it  has  been  suggested  that  the  triggering  agent 
may  not  be  the  density  wave  itself  but  the  shock  fronts 
produced  by  the  supermagnetosonic  collision  between  the  two 
streams.   This  front  would  propagate  forward  in  the  moving 
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wind  stream  at  the  fast-mode  velocity  and  thus  arrive  at 
Jupiter  ahead  of  the  density  wave  itself.   For  example,  a 
shock  which  left  the  wave  at  1  AU  in  a  wind  having  an  average 
bulk  speed  of  470  km/s  would  reach  Jupiter  about  2  days  before 
the  wave. 

The  ultimate  question  as  to  whether  or  not  the  corre- 
lations represent  a  physical  reality  is  reserved  to  the 
final  chapter.   There  are  other  items  of  evidence  yet  to  be 
presented. 


CHAPTER  5 
SPECTRAL  ANALYSIS 

Previous  Investigations 

In  general,  most  spectral  analysis  studies  of  the 
Jovian  decametric  emission  have  been  directed  toward  one  of 
two  goals:   the  measurement  of  the  Jovian  rotation  period 
or  the  measurement  of  frequencies  associated  with  the  lo 
effect.   The  first  of  these  studies  was  conducted  by  Shain 
(1956)  who  plotted  the  probability  of  emission  versus  the 
CML.   While  this  study  in  a  sense  was  not  really  a  frequency 
analysis,  it  did  indeed  reveal  the  presence  of  the  periodic 
source  structure.   Similar  studies  have  been  done  since  and 
improved  the  precision  of  this  result.   Using  essentially 
the  same  technique  Bigg  (1964)  plotted  the  probability  of 
emission  versus  ISC  to  discover  the  lo  effect.   Using  a 
rather  different  method  Duncan  (1965)  presented  strong  evi- 
dence confirming  Bigg's  discovery.   Duncan's  procedure  in- 
volved choosing  a  period  and  performing  a  superposed  epoch 
analysis  of  all  available  data.   He  then  defined  the  function 
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where  M  is  the  number  of  channels  used  in  the  analysis,  N  is 
the  total  number  of  total  events  in  the  study  and  N.  is  the 
number  of  events  per  channel.   Closer  inspection  will  show 
that  this  equation  is  nothing  more  than  the  equation  for 
the  variance,  a    ,   of  the  distribution  of  events  per  channel, 
normalized  so  that  the  average  event  per  channel  is  equal  to 
one.   In  other  words,  his  function  Z  is  the  square  of  the  RMS 
value  of  the  waves  that  would  be  seen  in  the  Chree  analysis. 
Except  for  the  square  root,  this  is  exactly  the  same  procedure 
used  in  Chapter  3  to  measure  the  growth  of  density  waves  in 
the  various  solar  wind  models.   Having  calculated  Z,  Duncan 
would  then  add  a  small  increment  to  the  period  used  in  the 
first  Chree  analysis  and  do  a  second  Chree  analysis  and 
calculate. Z  again.   By  many  iterations  of  this  procedure  he 
was  able  to  obtain  a  function  of  the  "waviness"  versus  period. 
As  the  period  of  Chree  analysis  approached  the  period  of  some 
component  in  the  data,  Z  would  climb  until  it  reached  a  peak 
when  the  period  was  equal  to  the  period  in  the  data.   Thus, 
his  function  Z(t)  represented  a  power  spectrum  analysis. 
Duncan  has  used  this  method  in  subsequent  papers  (1966,  1967, 
and  1971)  to  study  both  the  lo  effect  and  the  System  III  ro- 
tation period. 

The  method  has  one  serious  drawback,  however.   It  pro- 
vides no  means  of  distinguishing  between  the  fundamental, 
harmonics  or  sub-harmonics.   If,  for  instance,  the  period  at 
the  moment  was  identically  equal  to  the  fundamental  of  a  period 
in  the  data  there  would  be  one  full  cycle  of  that  wave  in  the 
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Chree  analysis  plot.   If,  on  the  other  hand,  the  period  were 
equal  to  twice  that  of  the  one  in  the  data,  there  would  be 
two  full  cycles  of  that  wave  in  the  plot.   The  calculation  of  Z, 
however,  would  give  the  same  value  in  both  cases. 

Bigg  (1969)  has  used  Duncan's  method  to  "discover"  two 
previously  undetected  Jovian  satellites.   Douglas  and  Bozyan 
(1970)  have  shown  that  both  of  these  periods  were  actually  due 
to  beats  between  harmonics  of  the  terrestrial,  Jovian,  and 
loian  periods,  not  to  "new  satellites." 

Kaiser  and  Alexander  (1972  and  1973)  used  a  somewhat 
different  technique  to  perform  their  spectral  analysis.   In 
an  effort  to  overcome  the  difficulties  imposed  by  the  un- 
equally spaced  observations  they  employed  a  two  step  process. 
They  would  first  pick  a  period  and  then  do  a  crude  first  order 
least  squares  fit  of  a  Fourier  series  to  their  data.   The  co- 
efficients of  the  series  provided  a  rough  estimate  of  the 
power  in  the  time  series  at  that  period.   They  would  adjust 
the  period  to  maximize  the  power  in  its  vicinity  and  then  do 
a  much  higher  order  fit  on  the  series  and  examine  its  spectrum. 
By  repeating  this  procedure  and  varying  the  period  under  con- 
sideration they  were  able  to  develop  a  plot  of  power  versus 
period.   In  their  1973  paper,  one  of  the  subjects  they  addressed 
was  the  question  of  a  period  in  the  Jovian  data  that  might 
indicate  solar  control.   Using  17  years  of  data  they  examined 
the  spectrum  in  the  region  near  25.5  days  for  evidence  of  a 
peak.   Their  result  was  negative  and  they  concluded, 
that,  therefore,  there  was  no  evidence  of  solar  control. 
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Their  conclusions,  however,  are  not  supported  for  two  reasons. 
First,  as  already  discussed,  the  principal  period  of  a  solar 
influence  would  be  expected  to  be  no  longer  than  about  13 
days,  not  25.   Secondly  any  such  effect  would  have  been 
washed  out  in  their  study  because  they  employed  17  years  of 
data,  whereas  the  solar  sector  structure  is  stable  for 
periods  on  the  order  of  several  months,  but  no  more. 

Analysis 

Autocorrelations  and  Fourier  Transforms 

It  has  been  long  known  that  a  time  series  is  related 
to  its  frequency  power  spectrum  by  a  Fourier  transform  process. 
There  are  two  general  approaches  which  may  be  employed  owing 
to  an  interesting  property  of  the  relationship.   If  S(t)  is  a 
time  series,  F{S)  is  its  Fourier  transform,  F*(S)  is  the 
transform's  complex  conjugate,  C(S)  is  the  autocorrelation 
function  of  the  time  series  and  P(f)  is  its  power  spectrum, 

then 

P(f)  =  aF*{S)P(S) 

and 

P{f)  =  3FIC{S)] 
where  a  and  g  are  constants.   This  is  to  say  that  one  may  obtain 
the  power  spectrum  of  a  time  series  by  either  obtaining  the 
square  of  the  modulus  of  the  Fourier  transform  of  the  series  or 
by  obtaining  the  Fourier  transform  of  the  autocorrelation  func- 
tion of  the  time  series.   (The  reader  is  referred  to  Blackman 
and  Tukey  (1958)  and  Bracewell  (1965)  for  a  detailed  discussion 
of  these  properties.)   The  procedures  are  entirely  equivalent  and  the 
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"transform  squared"  technique  is  often  the  choice  of  workers 
today.   However,  for  reasons  that  will  become  apparent  the 
"transformed-autocorrelation"  procedure  is  used  in  this  study. 

The  application  of  both  the  transform  and  the  auto- 
correlation has  limitations  when  dealing  with  real  data. 
Since  any  shift  other  than  zero  necessarily  neglects  points 
at  opposite  ends  of  the  time  series  and  its  image,  the  maximum 
shift  is  limited  by  the  finite  length  of  the  time  series.   The 
longest  period  one  might  study  with  the  autocorrelation  function 
is  on  the  order  of  the  maximum  shift  employed.   Ideally,  this 
maximum  shift  should  be  small  compared  to  the  length  of  the 
series.   Shifts  greater  than  20%  must  be  regarded  carefully. 
If  one  wishes  to  look  for  longer  periods,  the  best  method  is 
direct  examination  of  the  series. 

The  Fourier  transform  suffers  from  similar  problems. 
It  is  affected  by  the  finite  length  of  the  data  because  the 
transform  theoretically  has  integration  limits  at  positive  and 
negative  infinity.   Practically  speaking  the  transform  thinks 
it  is  operating  on  a  product  of  two  functions  when  it  deals 
with  finite  data  length.   These  two  functions  are  an  infinitely 
long  data  function  and  a  "window"  function  which  has  a  value 
of  1  only  where  the  actual  data  exist  and  a  value  of  0  every- 
where else  (this  particular  window  is  often  referred  to  as  a 
"boxcar"  function) . 

The  effect  of  the  window  function  is  that  the  resulting 
spectrum  is  the  convolution  of  the  data  spectrum  and  the 
window  spectrum.  The  simplest  window  one  might  choose  is  the 
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Figure  27.   Time  series  and  spectrograms  shov;ing  the 
effect  of  finite  boxcar  window.   The 
uppermost  plots  represent  an  infinitely 
long  cosine  wave  time  series  and  its 
spectrum.   The  center' plots  are  of  a  box- 
car time  series  and  its  spectrum.   The 
lov;est  plots  are  those  of  the  combination 
of  the  cosine  and  boxcar  time  series  and 
the  resultant  spectrum. 
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boxcar  function  mentioned  above.  The  transform  of  this  function 
is  of  the  form 

r:.  _    sin  X 
F  -  a 

Consequently,  if  the  data  spectrum  consisted  of  infinitesimally 
thin  (monochromatic)  lines,  the  effect  of  the  window  is  to 
widen  the  spectral  line  by  making  each  line  have  the  sinX/X 
shape  shown  in  Figure  27. 

The  most  striking  feature  of  Figure  27  is  that  in 
addition  to  widening  the  line  it  also  has  some  negative 
fringes!   Both  of  these  problems  can  be  controlled  somewhat 
by  choosing  a  more  subtle  window.   When  the  transformed  auto- 
correlation approach  is  used,  one  can  attack  the  window  problem 
either  in  the  time  domain  by  weighting  the  original  time  series 
by  something  other  than  the  boxcar  function,  or  in  the  inter- 
mediate (1/time)  domain  of  the  autocorrelation  function,  or 
in  the  frequency  domain  after  the  Fourier  transform-  By  any 
route,  one  is  effectively  introducing  a  smoothing  procedure 
which  acts  on  the  finally  derived  spectrum. 

There  are  a  number  of  window  functions  in  common  use. 
Perhaps  the  two  most  frequently  used  functions  are  the  Banning 
and  Hamming  functions.  The  forms  of  these  functions  vary  accord- 
ing to  which  step  in  the  time  series-spectrum  process  they  are 
introduced.   If  they  are  introduced  as  the  last  step  (after 
taking  a  boxcar  windowed  transform)  the  two  procedures  merely 
provide  a  3-point  weighted  smoothing  of  the  spectral  function. 
This  smoothing  corresponds  to 
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P^  =  0.25  P^_^  +  0.50  P.  +  0.25  P._^,    Banning 

P!  =  0.23  P.  T  +  0.54  P.  +  0.23  P.  ^    Hamming 
1         1-1         1         1+1  ^ 

where  P.  is  the  i   value  of  the  boxcar  windowed  power  spectrum. 
Figure  28  compares  the  spectra  of  several  windows. 

The  Banning  procedure  is  used  extensively  in  the 
present  study.   One  of  the  salient  features  of  Figure  28  is 
that  with  a  Banning  window,  the  half-power  width  of  the  win- 
dowed transform  is  two  channels.   Thus,  the  value  of  each 
point  in  the  power  spectrum  is  coupled  somewhat  to  each  ad- 
jacent channel  as  the  equation  above  indicates  and  thus  there 
are  (N-l)/2  independent  points  in  the  analysis.   One  might 
also  say  that  there  are  N  Banning  dependent  points. 

The  data  impose  a  second  restriction  on  the  transform 
in  that  they  are  not  continuous  but  consist  of  discrete 
samples.   Clearly  then,  if  the  samples  are  evenly  spaced, 
the  data  contain  a  strong  periodic  component — the  sampling 
frequency.   This  component  is  just  an  artifact  of  the  data 
collection  regime  and  should  not  be  allowed  to  show  up  in  the 
spectral  analysis.   The  standard  approach  to  this  problem  is 
to  limit  the  highest  frequency  analyzed  to  one-half  the  samp- 
ling frequency. 

Another  interesting  feature  of  the  analysis  is  the 
effect  the  process  has  on  the  precision  of  the  results.   The 
autocorrelation  function  in  this  study  has  values  for  each 
one-day  shift.   That  is  to  say,  the  function  consists  of 
evenly  spaced  values  at  one  day  intervals  in  the  time  domain. 
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Figure  28.  A  comparison  of  the  line  width  and  fringes 
with  boxcar,  Banning  and  Hamming  windows. 
The  Hamming  window  produces  a  narrower  main 
lobe  peak  than  the  Banning  window  but  its 
sidelobes,  though  initially  weaker,  damp 
out  much  slower  than  the  Manning  sidelobes. 
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The  application  of  the  Fourier  transform  yields  evenly  spaced 
values  in  the  frequency  domain.   As  a  result,  the  "period  band- 
width" associated  with  the  final  spectrum  is  directly  propor- 
tional to  the  period. 

This  effect  is  not  unreasonable  considering  that  the 
nature  of  the  transform  is  to  take  events  in  the  time  domain 
and  create  the  frequency  domain  image  of  those  events.   In 
the  process  the  evenly  spaced  points  in  the  time  domain  map 
into  evenly  spaced  points  in  the  frequency  domain.   Thus,  the 
higher  the  frequency  the  higher  the  precision  in  the  measurement 
of  period  associated  with  that  frequency.  The  spacing  of  the 
periods  in  Table  8  on  page  110  illustrates  the  feature. 

Basically,  the  question  here  is  one  of  recognizing  the 
independence  of  values.   One  can  compute  the  value  of  the 
power  density  for  any  frequency  within  the  limits  already  dis- 
cussed, but  one  should  only  be  interested  in  points  that  are 
marginally  coupled  to  each  other,  no  more  and  certainly  no 
less.   The  failure  to  recognize  this  requirement  led  Holmes 
(1965)  into  a  serious  error  in  his  study  of  Jovian  periods. 
He  followed  the  same  basic  procedure  as  that  employed  in  this 
study — the  autocorrelation-transform  approach.   However,  he 
chose  to  separate  his  calculations  at  equally  spaced  intervals 
in  time  (period).  As  a  result  his  "spectra"  contained  strongly 
coupled  or  dependent  points  in  the  longer  period  region  and 
omitted  a  significant  number  of  independent  points  in  the 
short  period  region. 
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In  the  present  study  the  maximum  autocorrelation  shift 
was  54  days.   Using  the  Fourier  transform  with  a  Banning  win- 
dow this  will  lead  to  55  Banning  dependent  points  in  the  fre- 
quency spectrum.   The  relationship  of  the  frequencies  and 
periods  is  shown  in  Table  8. 

Aliases 

While  it  is  desirable  to  avoid  problems  associated  with 
the  sampling  frequency  this  cannot  be  fully  achieved  unless — in 
addition  to  restricting  the  analysis  to  frequencies  less  than 
one-half  the  sampling  frequency — the  input  data  have  previously 
been  filtered  to  remove  all  frequencies  above  the  cutoff  fre- 
quency in  the  analysis.   Even  though  such  information  is  out- 
side the  range  considered  in  the  analysis,  the  frequency  of 
the  beat  or  heterodyne  between  high  frequency  information  and 
some  harmonic  of  the  sampling  period  will  fall  into  the  analysis 
range.   These  false  features  in  the  spectriom  are  known  as 
"aliases." 

Unfortunately,  the  filtering  of  the  Jovian  data  is 
impossible.   Both  the  rotation  of  Jupiter  and  the  revolution 
of  lo  occur  at  frequencies  above  the  0.5  cycles  per  day  (c/d) 
analysis  cutoff.   Consequently,  any  effort  to  filter  will 
cause  all  of  the  data  to  be  rejected!   As  a  result  aliases 
will  occur  at  e.   frequency  f  whenever 

^er\^^s   "  ^1  ^  i  ^s    ^  =  ^'  2,  3,... 
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where  f  is  the  sampling  frequency  and  f  is  the  "out -of -range" 

frequency.   In  the  case  of  the  present  study  this  relation 
becomes 

f^  =  |k'-  fl   <  0.5     k  =  1,  2,  3,... 

For  example,  the  System  III  synodic  rotation  frequency 
is  2.417  c/d.   Consequently,  we  may  calculate  the  alias  of 
the  fundamental  System  III  frequency 

f^  =  (2  -  2.417|  =  0.42  c/d 

which  corresponds  to  a  period  of 

p^  =  2.40d 

a. 

By  the  same  process,  the  alias  of  the  second  harmonic  of  the 
System  III  period  is  found  to  be 

f^  =  |5  -  4.8341  =  0.166  c/d 

Pg  =  6.02  d 

The  form  of  the  observation  function  0(t)  was  discussed  in 
Chapter  2.   There  it  was  noted  that  this  function  was  regarded 
as  the  product  of  three  periodic  rectangular  functions.   In 
the  context  of  this  discussion  the  earth  function  E(t)  repre- 
sents the  sampling  function  while  the  Jovian  and  loian  func- 
tions J(t),  I(t),  and  N(t)  represent  the  input  data. 
Therefore,  the  frequency  spectrum  of  the  products  J(t)I(t) 
and  J(t)N(t)  will  consist  of  the  fundamental  and  harmonic 
frequencies  of  each  component  in  the  product  plus  the  sums 
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Table  9 
Calculated  Aliases  for  J(t)I(t) 


Origin 

Frequency 
(c/d) 

Period 
(d) 

Frequency  Index 
(1/108  c/d) 

J 

0.417 

2.40 

45 

2J 

0.166 

6.02 

18 

I 

0.435 

2.30 

47 

21 

0.129 

7.77 

14 

J-I 

0.144 

6.94 

15-16 

J+I 

0.016 

62.50 

2 

J-2I 

0.285 

3.51 

30-31 

J+2I 

0.457 

2.19 

49 

I-2J 

0.263 

3.80 

28 

I+2J 

0.393 

2.55 

42 
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and  differences  of  those  frequencies.   Table  9  shows  the 
expected  frequencies  of  the  fundamental,  and  first  harmonic 
combinations  of  these  products-   Note  that  the  component 
frequencies  of  both  products  are  the  same  by  the  definition 
of  N(t)  . 

The  graph  in  Figure  29  displays  the  actual  spectrum 
of  the  observation  functions  for  source  B  at  18  MHz  during 
1967-68.   It  will  be  noted  that  peaks  at  virtually  all  of 
the  first  and  second  order  frequencies  in  Table  9  are  in 
fact  found  in  the  actual  functions. 

Timing  precision.   The  question  of  periods  naturally 
requires  a  statement  about  the  precision  of  their  measure- 
ment.  In  the  data  reduction  process  no  effort  has  been  made 
to  keep  track  of  the  exact  time  an  event  occurred. 
Consequently,  the  fundamental  unit  of  time  in  this  study 
is  one  day.   As  a  practical  matter  the  uncertainty  of  the 
time  of  a  storm  is  +6  hours  from  local  midnight.   Since  data 
from  the  entire  200-250  day  apparition  are  used  to  arrive  at 
the  various  periods  the  precision  is  on  the  order  of 

^  =  ±-§#-=±0.125% 

In  comparing  the  System  III  spectral  peaks  with 
theoretical  values  a  disturbing  discrepancy  arose.   The 
theoretical  alias  frequencies  had  been  calculated  based  on 
an  estimate  of  the  average  terrestrial  and  Jovian  synodic 
period  during  the  apparition.   In  making  the  estimate  the 
assumption  was  that  the  fluctuation  of  the  synodic  period 
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due  to  orbital  motion  was  approximately  sinusoidal.   The  average 
value  was  taken  to  be  the  RMS  value  over  the  apparition.   These 
values  are 

p^  =  23^55"*24^  +  l"'48^ 

p  =   g^ss'^Ol^  +  48^ 

vJ  — 

The  problem  was  that  the  fundamental  alias  of  the  Jovian 
rotation  frequency  is  predicted  by  these  periods  corresponds 
to  a  period  of  2.45d,  but  the  observed  peak  was  clearly  at 
2.40d  in  the  spectral  analysis.   Further  the  2.45d  figure  is 
significantly  out  of  the  precision  tolerance. 

It  was  recognized  that  if  the  peak  were  to  be  at  2.40d 
then  one  of  the  periods  above  was  wrong.   It  became  apparent 
that  one  reason  for  the  discrepancy  is  the  fact  that  in  the 
effort  to  avoid  daylight,  the  time  of  the  center  of  observing 
sessions  is  systematically  skewed.   At  the  beginning  of  an 
apparition  observations  are  made  for  a  short  time  in  the 
early  morning  hours.   As  the  apparition  progresses  the  obser- 
vations and  their  centroid  occur  earlier  and  earlier.   The 
net  effect  is  that  the  synodic  period  of  the  observatory  is 
lengthened  with  an  accumulated  difference  of  about  one  hour 
a  month.   If  one  assumes  that  the  terrestrial  period  is  wrong 
the  "correct"  value  can  be  computed  back  from  the  2.40d 
period  in  the  data.   Doing  so  yields 
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p  =  23  57  57 

TTI      O 

The  error  due  to  the  "station  lag"  is  1  58  .   There  was  a 

s 
second  error  of  19  due  to  an  incorrect  provision  for  Jupiter's 

orbital  motion  which  was  also  discovered,  leaving  a  total 

in  s 
correction  of  2  17  .   This  added  to  the  originally  calculated 

theoretical  value  gives 

=  23'^57"'41^  +  l"'48^ 
'^e  — 

This  figure  is  in  excellent  agreement  with  the  experimental 
value.   As  a  result,  the  experimental  value  is  used  in  all 
calculations  such  as  those  showing  the  expected  alias  fre- 
quencies in  Table  9. 

Amplitude  confidence.  The  establishment  of  confidence 
intervals  for  the  amplitude  of  spectral  density  measurements 
is  a  difficult  task.   If  the  data  are  that  of  a  random 

Gaussian  distribution,  the  sampling  distribution  is  approxi- 

2 

mately  a  x   distribution  of  order 

n  =  2N/M 

where  N  is  the  number  or  points  in  the  sample  and  M  is  the 
maximum  shift  or  lag  in  the  autocorrelation  (Blackman  and 
Tukey,  1958).   Current  statistical  theory  has  not  yet  re- 
solved the  problem  of  the  time  series  which  contains 
significant  non-random  components.   This 
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problem  is  exacerbated  by  the  relative  shortness  of  the 
time  series  (200-250)  points. 

If  the  time  series  did  have  a  random  Gaussian  distri- 
bution, the  confidence  interval  about  the  average  power 
density  P  is  given  by 


nP     ,  p  ,    nP 


-2  -  ^  ^  -2 

'^n,a/2         ^n,l-a/2 

Since  the  number  of  degrees  of  freedom  n  depends  directly  on 
the  number  of  points  in  the  series,  the  width  of  the  interval 
is  strongly  dependent  on  that  factor.   For  example,  the  95% 
confidence  interval  for  216  points  and  a  54-day  maximum  shift 
is  given  by 

n=  2  x-2|f-=  8 

0.46  P  1  P  <  3.67  P      95%  confidence. 

Thus  with  216  points,  the  interval  is  very  wide,  spanning  a 
range  of  8:1.   On  the  other  hand,  if  the  sample  contained 
2160  points 

n  =  2  X  -i|f-  =  80 
0.75  P  £  P  £  1.40  P 
which  results  in  a  range  of  2:1. 

In  extreme  cases  with  non-random  functions,  such  as 
those  whose  spectra  are  shown  in  Figure  29,  common  sense 
dictates  that  the  line  features  reflect  the  non-random  com- 
ponent.  Thus  one  can  recognize  the  "floor"  between  peaks 
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as  representing  the  random  component.   When  this  occurs,  one 
can  apply  the  confidence  interval  relationship  to  the  floor 
values  only.   However,  it  is  not  obvious  where  to  apply  this 
relationship  when  confronted  with  a  spectrum  like  that  in 
Figure  30  or  Figure  32.   This  is  indeed  unfortunate  because 
such  situations  are  precisely  those  for  which  the  intervals 
are  most  necessary. 

An  extremely  conservative  approach  would  be  to  apply 
the  confidence  interval  relationship  using  the  overall 
average — including  possible  non-random  peaks.  This  procedure 
will  always  overestimate  the  width  of  the  interval  with  non- 
random  data. 

A  more  liberal  approach  would  be  to  define  local  minima 
as  representing  floor  points.   Thus  one  would  consider  the 
average  of  the  random  component  to  be  the  average  of  the 
minima. 

The  Galileo  Club 

At  first  sight,  one  of  the  problems  that  spectral 
analysis  might  hope  to  address  is  that  of  the  existence  of 
additional  satellite  exciters.   Unfortunately,  this  will 
prove  a  difficult  task.   Almathea  and  the  Galilean  satellites 
are  the  only  ones  which  are  within  the  Jovian  magnetosphere 
more  or  less  continuously.   It  seems  reasonable,  therefore, 
to  restrict  attention  to  this  group.   Almathea  is  the  inner- 
most moon  and  is  quite  small  (250  km  in  diameter) .   The  four 
Galilean  satellites,  lo,  Europa,  Ganymede  and  Callisto  are 
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all  quite  a  bit  larger,  ranging  from  3000  to  5000  km  in 
diameter. 

The  Galilean  satellites  form  a  harmonically  locked 
orbital  system.   Apparently  after  eons  of  perturbations  they 
have  assumed  orbital  periods  which  are  related  by  rational 
numbers.   Table  10  demonstrates  this  relationship. 

Table  10 
Orbital  Periods  of  the  Galilean  Satellites 


Satellite 

Period 
(h) 

Ratio 

with  Preceding 
Period 

lo 

42.47 

_ 

Europa 

85.23 

2:1 

Ganymede 

171.71 

2:1 

Callisto 

400.53 

7:3 

As  might  be  expected,  this  situation  causes  a  number  of 
redundancies  in  the  various  periods  in  the  power  spectrum. 
Fortunately,  the  frequencies  of  all  but  lo  and  Almathea  are 
lower  than  the  0.5  c/d  cutoff  frequency  so  aliases  will  be 
less  of  a  problem.   Expeifience  has  shown  the  third  order 
alias  effects  are  negligibly  small  and  thus  they  will  be  ignored, 
Table  11  shows  the  expected  frequencies  for  the  various 
satellites.   In  scanning  the  table  the  redundancies  become 
apparent.   lo,  Europa,  and  Ganymede  all  have  predicted 
periods  in  the  range  of  frequency  index  14  to  16,  while 
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Europa  and  Ganymede  have  components  in  the  30  to  32  range. 

The  plots  in  Figures  30,  31,  and  32  depict  the  frequency 
spectra  of  the  daily  activity  index  Ag(t)  as  defined  in  Chapter 
2  for  1967-68.   In  principle,  this  function  is  free  of  the 
observation  spectra  previously  displayed  in  Figure  29.   In 
these  plots  the  lo  effect  has  deliberately  not  been  removed 
in  order  to  demonstrate  the  effect  of  a  periodic  exciter  on 
the  data.   Source  B  has  long  been  recognized  as  the  prototype 
for  lo  related  activity.   These  plots  show  clearly  that  this 
is  the  case.   The  lo  features  in  the  vicinity  of  indices  14-16, 
30-31  and  47  are  the  dominant  features  in  the  source  B  spec- 
trum.  All  but  the  strong  peak  at  index  17  can  be  attributed 
to  the  loian,  Jovian,  and  joint  interactions  outlined  in 
Table  9.  The  same  pattern  is  also  evident  in  the  source  A 
plot,  though  the  unexplained  index  feature  is  absent.   Most 
of  these  periods  are  absent  in  the  source  C  plot.   The 
strongest  peak  in  the  C  spectrum  corresponds  to  the  principal 
System  III  alias  at  2.40  days.   This  period  has  supposedly 
been  removed  from  the  data  by  the  compensation  procedure  which 
defines  As(t).   This  point  will  be  addressed  shortly. 

A  Matter  of  Priorities 

Tables  9  and  11  illustrate  that  the  number  of  true  and 
alias  frequencies  which  might  show  up  in  the  occurrence  spectrum 
is  large.   Nevertheless,  it  is  clear  that  there  is  a  descending 
order  of  importance  as  the  various  interactions  become 
more  complex.   In  this  regard  it  is  worthwhile  to 
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Table  11 
Sidereal  Frequencies  of  the  Five  Innermost  Satellites 


Satellite 

Origin 

Frequency 
(c/d) 

Period 
(d) 

Frequency  Index 
(1/108  c/d) 

Almathea 

A 

0.020 

50.92 

3 

2A 

0.040 

25.00 

4 

A-J 

0.082 

12.19 

9 

A+J 

0.090 

11.09 

9 

lo 

I 

0.435 

2.30 

47 

21 

0.129 

7.77 

14 

I-J 

0.144 

6.94 

15-16 

I+J 

0.016 

62.5 

2 

Europa 

E 

0.281 

3.57 

30-31 

2E 

0.438 

2.29 

47 

E-J 

0.139 

7.20 

15 

E+J 

0.299 

3.35 

32 

Ganymede 

G 

0.140 

7.17 

15 

2G 

0.279 

3.58 

30 

G-J 

0.280 

3.57 

30 

G+J 

0.440 

2.27 

48 

Callisto 

C 

0.060 

16.71 

7 

2C 

0.120 

8.36 

13 

C-J 

0.360 

2.78 

39 

C+J 

0.480 

2.08 

52 
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formally  recognize  this  hierarchy. 

All  fundamental  frequencies  will  be  classified  as  being 
of  "first  order."  Frequencies  arising  due  to  the  cross-product 
interaction  (sum  or  difference)  between  two  fundamentals,  in- 
cluding a  self- interaction,  are  classed  as  "second  order." 
Those  frequencies  which  occur  due  to  a  three-fold  fundamental 
interaction,  including  self-interactions,  will  be  termed  as 
"third  order"  effects. 

A  priori  one  would  expect  first  order  effects  to  dominate 
the  spectra  of  a  source.   These  influences  are  precisely  what 
the  compensation  procedure  leading  to  A  (t)  was  intended  to 
remove,  at  least  for  Jupiter  and  lo.   If  this  goal  has  been 
achieved  then  the  remaining  periods  would  be  of  second,  third 
and  higher  orders,  or  they  would  be  first  order  effects  of 
influences  not  yet  considered. 

To  this  point  in  the  discussion  third  order  effects 

have  been  neglected  as  insignificant.   However,  if  the  first 

order  effects  can  be  removed,  then  the  higher  order  influences 

obviously  descend  one  order  and  thus  become  more  important. 

In  the  examination  of  the  spectra  which  follows  the  hierarchy 

of  interactions  will  be  used  to  overcome  the  problems 

associated  with  such  things  as  redundant  frequencies  due  to 

the  harmonically  locked  Galileo  Club.   This  will  be  done  by 

the  invocation  of  three  "golden"  rules. 

1.   The  fundamental  frequencies  and  their  harmonics 
of  all  orders  for  E(t),  I(t),  and  N(t)  will  not 
be  seen,  nor  will  their  aliases,  due  to  the  A  (t) 
procedure.  ^ 
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2.  A  cross-product  (joint  interaction)  frequency  or 
its  alias  will  only  be  seen  if  at  least  one  of 
its  constituents  has  not  been  removed  by  the 
Ag(t)  procedure. 

3.  A  higher  order  frequency  or  its  alias  will  not 
be  seen  if  lower  order  frequencies  or  their 
aliases  composed  of  the  same  constituents 

or  their  subharmonics  are  not  seen. 

Figures  33  through  39  show  the  occurrence  spectrograms 
of  the  non-Io  component  of  each  of  the  sources  at  18,  22,  and 
27  MHz  for  1967-68  and  18,  22,  26.3,  and  27  MHz  for  1974. 
The  small  amount  of  data  available  in  some  of  these  source- 
frequency  combinations  has  rendered  their  particular  plots 
useless,   specifically,  all  of  the  source  B  plots  and  the 
1967  source  A  at  27  MHz  plot.   The  average  source  observation 
time  represented  on  these  plots  was  about  1400  zones  over  the 
apparition  and  yet,  due  to  the  small  activity  probability, 
fewer  than  25  zones  were  observed  storming.  Thus,  activity 
was  limited  to  only  a  few  events.   The  effect  of  this  condition 
is  visible  in  most  of  these  plots. 

The  26  MHz  plots  are  also  viewed  with  suspicion.   While 
the  activity  probability  was  much  higher,  the  large  array 
was  operational  for  only  about  620  zones  of  observations — less 
than  half  the  time  logged  by  the  low  sensitivity  systems;  thus 
there  are  many  gaps  in  the  time  series  where  there  are  no  data. 
Consequently,  the  bulk  of  attention  will  be  directed  to  the 
source  A  and  source  C  low  sensitivity  systems. 

The  spectrograms  have  been  plotted  against  an  arbitrary 
power  scale  of  from  0  to  5  units.   If  one  adopts  the  view  that 
the  minima  in  the  function  represent  the  floor  or  background 
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Figure  33.   Source  A,  1967-68  activity  spectra.   The 
27  MHz  spectrum  is  data  limited.   Its 
apparent  periodic  nature  is  an  artifact 
produced  by  the  autocorrelation  of  two 
events. 
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Figure  34.   Source  B,  1967-68  activity  spectra.   All 
frequencies  are  data  limited. 
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Figure  35.  Source  C,  1967-68  activity  spectra.  The  22 
and  27  MHz  plots  are  contaminated  by  System 
III;  see  text. 
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Figure  36.   Source  A,  1974  activity  spectra. 
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Figure  37.   Source  B,  1974  activity  spectra.   All 
frequencies  are  data  limited. 
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Figure  38.   Source  C,  1974  activity  spectra.   The  27 

MHz  spectrum  is  contaminated  by  System  III. 
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Figure  39.   26.3  MHz  activity  spectra  1974.   All  sources 
are  data  limited. 
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of  random  noise,  then  the  background  level  is  1  or  less  in 
the  source  A  and  source  C  plots.  Thus,  the  90%  confidence 
interval  is 

0.52  <^  P  <^  2.93. 

Contamination 

It  turns  out  that  the  first  golden  rule  can  be  violated 
by  at  least  one  odd  circumstance.   Examination  of  the  source 
C  spectrograms  reveals  statistically  significant  maxima  at 
index  45  (2.4  days)  in  1967-68  at  22  and  27  MHz  and  in  1974 
at  27  MHz.   Low  level  vestiges  of  this  feature  are  found  in 
all  the  remaining  source  C  plots. 

Inasmuch  as  index  45  is  the  alias  of  the  System  III 
rotation  fundamental  frequency,  its  discovery  in  the  spectro- 
gram caused  considerable  concern. over  the  validity  of  the 
removal  of  Jovian  and  loian  periods  by  the  Aj,(t)  procedure. 
The  success  of  the  procedure  with  source  A  allayed  the 
fear  but  did  not  explain  the  presence  of  the  spectral  line. 

Examination  of  the  input  data  to  the  autocorrelation- 
spectrum  routine  ("Jupcor  III",  see  Appendix)  revealed  the 
cause  of  the  problem.   Despite  the  fact  that  the  non-Io 
component  of  source  C  has  a  relatively  small  emission  proba- 
bility in  terms  of  the  percentage  of  observed  time  spent 
storming,  it  storms  quite  frequently.   For  the  source 
frequency  combinations  showing  the  intense  line,  about  6  days 
out  of  10  produced  a  storm,  perhaps  better  a  "spurt",  since 
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each  was  only  one   zone  long!   As  a  result,  the  time  series 
consisted  largely  of  l/Og(t)  which  naturally  has  the  same 
period  as  0- (t) . 

The  explanation  for  the  state  of  the  data  is  uncertain. 
The  phenomenon  occurs  only  in  source  C  and  it  becomes  more 
prominent  with  increasing  frequency,   it  occurred  in  both 
1967-68  and  1974.   It  is  tempting  to  conclude  that  the  source 
is  active  frequently  but  only  the  peak  of  a  storm  is  powerful 
enough  to  be  detected  by  low  sensitivity  (yagi)  antenna 
systems.   If  this  is  the  case,  the  storms  themselves  must, 
nevertheless,  be  short  because  comparing  the  1974  26  MHz  high 
sensitivity  activity  index  to  the  27  MHz  low  sensitivity  index 
(Table  3)  shows  an  increase  in  the  activity  index  by  a  factor 
of  only  2.76.  However,  the  26  MHz  spectrogram  has  only  the 
slightest  hint  of  the  frequency  index  45  spectral  line. 

In  any  case  the  situation  is  some  kind  of  data  limiting 
problem  and  must  be  regarded  as  an  artifact.   Thus,  caution 
must  be  used  in  the  interpretation  of  any  other  features  in 
the  spectrograms  of  source  C  22  and  27  MHz  in  1967-68  and 
27  MHz  in  1974. 

Solar  Control 

It  was  recognized  early  in  this  work  that  the  auto- 
correlation and  power  spectral  functions  view  periods  in  data 
from  different  perspectives.   A  clear  example  of  this  is  seen 
by  comparing  Figures  40  and  41.   In  Figure  40  both  the  auto- 
correlation and  power  spectral  functions  describe  a  fairly 
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monochromatic  13-day  period  in  the  solar  wind  particle  flux 
density  as  predicted  by  the  field-free  model  for  1967-68.   On 
the  other  hand,  the  autocorrelation  function  of  the  1974  simple- 
field  flux  density  shown  in  Figure  41  has  its  largest  peak  at 
25  days,  while  the  spectral  function  shows  the  principal  power 
peak  near  13  days. 

The  difference  between  the  1974  autocorrelogram  and 
the  spectrogram  is  not  a  contradiction.   The  velocities  of 
the  two  high  speed  streams  in  1974  were  noticeably  different. 
Due  to  this  asymmetry  (refer  back  to  Figure  8)  the  best  auto- 
correlation occurred  at  25  days  even  though  the  principal 
period  was  13  days.   Nevertheless,  the  spectrogram  shows  very 
little  power  at  25  days.   The  transformed-autocorrelation 
technique  for  obtaining  power  spectra  was  used  in  order  to 
have  the  advantage  of  both  viewpoints. 

No  assijunptions  have  been  made  about  the  nature  of  the 
density  wave  or  shock  front  interaction  with  the  Jovian  magne- 
tcsphere  nor  about  the  workings  of  the  radiation  mechanism(s) . 
As  a  result  it  is  not  known  whether  the  asymmetry  in  initial 
wind  velocities  is  reflected  in  the  observed  Jovian  radiation. 
It  is  possible,  for  instance,  that  some  threshold  process  might 
"accept"  one  wave  and  emit  while  "rejecting"  the  other 
and  remaining  silent.   Alternatively,  one  wave  might  bend  the 
beam  of  radiation  away  from  the  earth  while  the  other  might 
direct  the  beam  at  the  earth.   Therefore,  the  possibility  of 
a  25-  or  26-day  period  in  the  Jovian  emission  when  the 
velocity  is  highly  asymmetric  should  not  be  ignored. 
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In  the  1967-68  Jupiter  spectrograms  none  of  the  plots 
show  statistically  significant  levels  at  either  13.5  or  27 
days  (the  nearest  points  to  13  and  26  days) ,   though  several 
show  higher  than  average  power  at  13  days.   The  1974  functions 
are  somewhat  different.   The  cross-correlations  in  Chapter  4 
showed  that  source  A  correlated  with  the  wind  at  26  and  27  MHz. 
The  spectra  for  both  these  source-frequency  combinations  have 
high  power  levels  at  13  days  though  they  are  not  outside  the 
confidence  interval.  The  source  C  cross-correlations  for 
1974  had  positive  peaks  on  18,  22,  and  26  MHz.   In  all  three 
cases,  the  spectrograms  show  statistically  significant  maxima 
at  27  days.   Furthermore,  at  27  MHz  where  there  was  no  cross- 
correlation,  the  27-day  maximum  is  present  but  much  weaker. 

Satellite  Control 

Excluding  the  3  peaks  caused  by  System  III  contamination, 
there  are  24  features  of  statistical  significance  in  the  11 
source- frequency  combinations  which  are  being  examined  at 
this  point — sources  A  and  C  at  18  and  22  MHz  in  1967-68  and 
18,  22,  and  27  MHz  in  1974.   In  order  to  visualize  these 
features  a  histogram.  Figure  42,  has  been  constructed  showing 
the  distribution  of  spectral  index  numbers  which  exceed  90% 
confidence  for  all  source  A  functions,  all  source  C  functions, 
and  the  sum  of  the  A  and  C  functions.   The  source  A  diagram 
shows  two  prominent  features,  one  at  about  indices  17  and  18 
and  the  other  at  about  index  40.   The  source  C  histogram  shows 
only  the  index  17  peak.   The  composite  diagram  shows  three 
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Figure  42.   Histograms  of  90%  confidence  spectral  features. 
1967-68  and  1974  values  are  merged  for  each 
source . 
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dominant  features,  one  at  index  17,  one  at  index  40  and  one 
at  index  4.   The  index  4  maximum  corresponds  to  the  27-day 
period  mentioned  in  connection  with  the  solar  correlation. 

The  index  17  peak  is  rather  interesting.   It  is  near, 
but  essentially  uncoupled,  from  the  second  order  lo  peak  (J-I) 
that  straddles  indices  15  and  16.   Though  the  entire  line 
profile  about  17  does  include  power  at  15  and  16,  it  is 
weaker  and  there  is  no  contribution  at  14 — the  site  of  another 
equally  strong  second  order  lo  peak  (21).   Since  the  Jovian 
frequencies  are  well  suppressed  not  only  in  the  composite 
histogram  but  in  the  individual  spectrograms  which  provide 
the  indices  15  and  16  contribution,  rule  2  says  that  if  15 
and  16  are  Jupiter-Io  contaminants  (J-I)  the  source  must  be 
unsuppressed  lo  information.   However,  the  21  peak  is  missing 
as  it  should  be  by  rule  3  because,  more  importantly,  the 
fundamental's  alias  at  index  47  is  missing.   The  only  infor- 
mation at  47  in  the  composite  histogram  was  contributed  by 
only  1  of  the  7  spectograms  which  contributed  to  index  17, 
and  it  is  probably  not  the  lo  fundamental.   As  a  result,  we 
can  rule  out  lo  contamination  as  the  source  of  index  17,  or 
the  whole  feature  from  15  to  20  for  that  matter. 

Another  curious  aspect  of  index  17  is  that  mentioned 
in  relation  with  the  activity  spectograms,  Figures  30  and 
31,  where  all  the  lo  activity  was  left  in.   Specifically, 
the  source  B  spectrum  has  a  sharp  independent  peak  to  the 
side  of  the  lo  indices  14,  15,  and  16  maxima  but  this  peak 
is  not  found  in  the  corresponding  source  A  lo  activity  line 
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profile.   This  strongly  suggests,  first  that  the  index  17 
line  is  capable  of  being  quite  strong — rivaling  some  lo  peaks — 
and  second,  in  addition  to  all  other  arguments,  it  is  not  an 
lo  feature  because  it  goes  away  in  source  A. 

Almathea  and  the  other  Galilean  satellites  have  been 
examined  to  see  if  they  can  account  for  the  peak.   The  only 
possible  candidate  for  a  simple  satellite-System  III  process 
analogous  to  the  lo  effect  is  Ganymede.   But  the  fundamental 
from  Ganymede  falls  not  on  index  17  but  index  15  and  thus 
it  seems  to  be  an  unlikely  source  of  the  "17  effect."  Joint 
interactions  between  two  and  three  satellites  have  been  ex- 
plored through  the  third  order  with  only  one  promising  candi- 
date, lo-Europa. 

Table  12  displays  the  expected  frequencies  of  the  lo- 
Europa  system.   It  is  seen  that  one  of  the  two  lowest  order 
frequencies  corresponds  to  index  17.   The  remaining  second 
and  third  order  frequencies  correspond  reasonably  well  to 
the  secondary  features  in  the  vicinities  of  indices  14,  30, 
and  47.  ~ 

The  index  40  line  is  not  associated  with  any  Jovian 
or  lo  line  except  for  a  third  order  line  (2J+I)  at  index  42. 
This  third  order  line  is  not  even  visible  in  the  A  and  B 
activity  spectra  and  is  too  far  removed  from  index  40  to  have 
any  connection  with  the  index  40  peak.   In  fact,  no  combination 
of  satellite  and  Jovian  interactions  studied  in  the  search  out- 
lined above  can  account  for  the  line.  Index  40  corresponds  to 
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Table  12 
Frequencies  of  the  lo-Europa  System 


Source 

Frequency 
(c/d) 

Period 
(d) 

Index 
(1/108  c/d) 

I+E 

0.155 

6.46 

17 

I-E 

0.283 

3.53 

31 

J+I+E 

0.263 

3.80 

29 

J+I-E 

0.300 

3.33 

32 

J- I+E 

0.430 

2.32 

46 

J-I-E 

0.133 

7.54 

14-15 

2.70  days.   The  line  is  seen  only  in  source  A  but  it  is  seen 
in  both  1967-68  and  1974.   In  both  cases  it  was  seen  at  all 
frequencies  though  it  was  not  plotted  for  1967-68  MHz  because 
those  data  are  suspect  and  it  was  not  plotted  for  1974  to 
18  MHz  because  it  just  failed  to  reach  the  90%  confidence 
limit.   For  now,  the  source  of  the  index  40  line  is  a  mystery. 


CHAPTER  6 
OCCURRENCE  DISTRIBUTIONS 

Statistical  Approach 

For  many  populations  it  is  either  impossible  or  im- 
practical to  measure  desired  parameters  of  every  member  of  the 
population.  One  of  the  basic  objectives  of  modern  statistics 
is  to  enable  one  to  infer  from  a  sample  certain  characteristics 
about  the  population  from  which  the  sample  is  dravm. 

In  the  case  of  the  Jovian  decametric  radiation,  the 
stroboscopic  effect  of  the  three  rotating  and  revolving 
bodies,  the  earth,  lo,  and  Jupiter,  conceals  much  of  the 
emission  population  from  the  observer.  One  must  be  content  with 
the  samples  which  are  available.   The  preceding  chapters  have 
dealt  basically  with  the  question  of  periodic  or  quasi- 
periodic  properties  of  the  radio  emission.  Periodicity  is 
one  form  of  non-random  process  and  evidences  of  this  property 
have  been  offered.  The  focus  in  this  segment  of  the  work  is 
the  attack  on  the  general  question  of  randomness. 

The  Service  Problem 

If  one  assumes  that  the  Jovian  emissions  stem  from  some 
truly  random  process,  then  the  characteristics  of  that  process 
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can  be  modeled  in  a  fashion  analogous  to  many  terrestrial 
situations.   For  example,  the  likelihood  of  random  Jovian 
emission  during  any  given  period  of  time  would  behave  the 
same  as  the  likelihood  of  a  bank  customer  approaching  a 
teller  for  service  during  any  given  period  of  time.   In 
both  situations  the  radio  astronomer  and  the  teller  are 
passive  observers  viewing  presumably  random  events. 

If  time  may  be  subdivided  into  small  periods  so  that 
the  probability  of  an  event  occurring  in  that  subinterval  is 
small,  if  the  probability  of  occurrence  is  constant  in  time, 
if  occurrences  are  independent  and  almost  never  occur  twice  in 
one  subinterval,  the  probability  distribution  of  these  events 
is  a  Poisson  distribution.   To  the  extent  the  experimental 
distribution  departs  from  the  Poisson  one  or  more  of  these 
conditions  are  not  being  met. 

Thus  if  Jupiter's  emissions  are  random  independent  events 
with  a  constant  probability  of  occurrence,  the  distribution  of 
the  number  of  storms  detected  per  unit  time  (say  a  few  5° 
CML  zones)  should  obey  a" Poisson  distribution. 

Furthermore,  if  the  requests  for  service  are  random  and 
the  time  of  completion  is  short,  the  lengths  of  the  service 
times — the  storm  lengths  in  this  analogy — obey  an  exponential 
distribution  (see  for  instance  Hartnet,  1970).   Consequently, 
one  expects  that  the  distribution  of  occurrences  and  the 
distribution  of  storm  lengths  will  yield  information  about 
the  underlying  processes. 
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Unfortunately,  there  are  some  potential  problems  in  pur- 
suing this  course.   (As  an  aside,  it  must  be  said  that  this 
particular  effort — distributions —  is  the  newest  and  hence  the 
least  developed  part  of  this  study.)   The  key  feature  is  the 
testing  of  the  concept  of  constant  probability  (i.e. ,  not  time 
dependent).   To  do  this  adequately  one  must  be  certain  that 
there  are  not  other  time  dependent  influences  on  the  probability 
of  detection  of  Jovian  radio  storms.   Of  course,  there  are. 
Due  to  interference  and  degraded  ionospheric  transparency,  the 
detection  of  activity  toward  the  end  of  an  apparition  is  more 
difficult.   Lebo  (1964)  has  considered  this  and  several  related 
effects  and  has  shown  that  they  all  have  measurable  influences. 
Equipped  with  this  knowledge  the  next  logical  question  is:   Can 
the  influence  of  these  effects  be  removed?  Heuristically  it 
seems  that  they  can,  but  it  appears  that  the  burgeoning  science 
of  statistics  has  yet  to  deal  fully  with  this  problem.   A 
further  discussion  of  this  particular  situation  and  its  possible 
solutions  will  be  taken  up  in  the  final  chapter.   For  the 
present,  a  less  sensitive,  more  robust  aspect  of  this  area  will 
be  pursued. 

Storm  Lengths 

Due  to  the  various  data  limiting  problems  mentioned  in 
the  preceding  chapter  the  following  presentation  is  limited  to 
a  merge  of  the  source  A  data  for  18,  22,  and  27  MHz  in  1967-68 
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and,  separately,  in  1974.   Figures  43  and  44  show  the  distribu- 
tion of  observation  lengths  (in  zones)  and  storm  lengths  (in 
zones).   The  observation  spectra  in  both  cases  are  entirely 
reasonable.   They  show  that  the  distribution  of  zones  observed 
is  essentially  constant  from  1  through  18.   But,  since  there 
are  19  zones  in  the  source,  corresponding  to  about  2.6  hours 
of  observing  time,  about  30%  of  the  time  the  observing  station 
"saw"  all  19  source  zones.   As  a  result  there  is  a  strong  peak 
at  19. 

Inspection  of  the  observation  length  and  storm  length 
distributions  for  both  1967-68  and  1974  reveals  that  80%  of 
the  observations  lasted  6  zones  or  more  while  80%  of  the  storm 
lengths  were  6  zones  or  less.   Thus  there  was  a  20%  probability 
that  the  20%  of  the  storms  which  exceeded  6  zones  were  trun- 
cated by  short  observations.   Put  another  way,  the  overall 
probability  of  observed  storms  being  truncated  by  the  beginning 
or  end  of  the  observation  of  the  source  is  only  4%.   As  a 
consequence,  the  storm  length  distributions  in  Figures  43  and 
44  are  good  estimates  of  the  true  distribution  of  storm  lengths 
at  observable  flux  levels. 

If  the  storms  are  produced  by  a  random  process  (that  is, 
one  assumes  they  have  a  Poisson  distribution)  and  the  storm 
lengths  are  relatively  short  and  of  random  length,  the  storm 
length  should  have  an  exponential  distribution.   Specifically, 

f(n)  =  Ae"'^^ 
where  \   is  the  reciprocal  of  the  average  value  of  the  storm 
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lengths  and  n  is  the  zone  number.   However,  in  the  present 
problem  it  has  been  defined  that  the  shortest  possible  storm 
is  one  zone  long.   The  equation  above  assumes  that  the  minimum 
value  is  zero  zones.   This  is  rectified  by  a  simple  translation 
of  axes  so  that  one  obtains 

f(n)  =  Ae-^("-l^ 

while  the  average  value  of  the  storm  length  must  be  reduced  to 
n-1.   These  values  are  calculated  directly  from  the  experimental 
distributions  and  are  found  to  be  2.16  in  1967-68  and  2.07  in 
1974.   Thus,  the  values  for  A  are 

0.463   in   1967-68 
and 

0.483   in  1974. 

The  theoretical  and  actual  storm  length  distributions 
are  shown  in  Figure  45.   The  actual  storm  lengths  differ  from 
those  in  Figures  43  and  44  only  in  that  they  have  been  un- 
normalized  to  show  the  actual  storm  counts  for  each  length. 
It  will  be  noted  that  the  number  of  single  zone  storms  falls 
short  of  the  predicted  value  in  both  cases.   Generally  speaking 
the  actual  number  of  longer  storms  exceeds  the  predicted  value. 
If  apparition-phase  edge  effects  were  significant  in  this  plot, 
just  the  opposite  would  be  expected  since  the  storms  would  be 
artificially  shortened  as  ionosphere  or  interference  per- 
mitted only  the  strongest  part  of  a  storm  to  be  observed. 
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Confidence  intervals.   If  one  has  a  table  of  observed 
and  theoretical  values  for  a  frequency  distribution,  the  samp- 
ling distribution  of  the  statistics 

~   M    (X.-T.)^ 


1^  T. 

where  X^  are  the  observed  values,  T.  are  the  theoretical  values 
and  M  is  the  number  of  values  observed,  is  a  chi-square  distri- 
bution with  M-1  degrees  of  freedom.   This  statistic  can  be  used 
to  determine  confidence  intervals  for  the  "goodness  of  fit" 
between  the  observed  and  theoretical  values.   It  is  worth  re- 
marking that  it  is  imperative  that  the  frequency  distribution 
must  be  in  the  un-normalized  form  as  in  Figure  45  because  the 
confidence  interval  in  the  chi-square  test  is  very  dependent 
on  the  actual  number  of  events  observed.   Normalization  of  the 
frequency  distribution  as  in  Figures  43  and  44  obscures  this 
information  and  leads  to  a  nonsense  result  in  the  test. 

Application  of  the  chi-square  test  to  the  storm  length 
distribution  leads  to 

2        _ 
^  1967-68  ~  ^^'^^ 

X^1974    =  110.35. 

Both  values  far  exceed  the  99.5%  confidence  interval  and  thus 
both  observed  storm  length  distributions  are  considered  non- 
random. 
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Chree  Analysis 

In  Chapter  4 ,  Chree  analysis  was  used  extensively  to 
compare  Jovian  events  with  modeled  solar  wind  events.   In  the 
process  many  superposed  epoch  diagrams  of  the  Jovian  data 
similar  in  character  to  those  shown  in  Figure  46  were  gene- 
rated. A  persistent  feature  seen  in  most  of  these  plots  was 
the  four  peaked  structure  seen  there. 

The  universality  of  these  features  has  led  to  an  investi- 
gation of  the  probability  that  these  modulations  could  occur 
by  chance.  Each  point  in  the  analysis  plot  represents  the 
mean  of  a  sample  drawn  from  the  population  of  emissions.   If 
one  assumed  a  null  hypothesis  that  the  population  had  a 
Gaussian  distribution,  the  population  of  the  means  of  all 
samples  of  that  population  would  be  a  t  distribution.   Thus, 
one  could  use  that  distribution  to  check  the  variation  from 
the  long-term  sample  mean.   However,  the  assumption  of  a 
Gaussian  distribution  can  be  eliminated  by  again  resorting 
to  the  chi-square  test. 

In  applying  the  chi-square  procedure  to  the  Chree 
analysis  there  is  a  problem  in  that  the  number  of  zones  ob- 
served is  not  the  same  for  each  day  in  the  epoch.   However, 
this  will  be  sidestepped  by  assuming  that  all  days  have  the 
same  sample  size.   The  goodness  of  fit  test  will  then  be  done 
on  the  resulting  distribution  of  activity  zones  versus  days. 
Again  attention  will  be  directed  to  source  A  at  18,  22,  and 
27  MHz  during  1967-68  and  1974. 
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Figure  46.      1967-68   source  A  activity   index  Chree 
analyses. 
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Under  ordinary  circtunstances  the  25-day  period  of  the 
analysis  would  result  in  24  degrees  of  freedom.   However,  in 
the  cross-correlation  procedure  the  Chree  analysis  tables 
were  Banning  smoothed  before  correlating.   The  smoothed 
tables  will  be  used  here  as  well,  thus  leaving  12  degrees  of 
freedom. 

The  results  of  the  chi-square  test  are  summarized  in 
Table  13. 

Table  13 
Chree  Analysis  Confidence  Intervals 


Frequency  x^         Confidence 

(%) 


1967-68 
18 
22 
27 

1974 
18 
22 
27 


44.2 

>99.5 

25.4 

>97.5 

32.7 

>99.5 

89.7 

>99.5 

55.7 

>99.5 

121.2 

>99.5 

Since  all  values  of  the  statistic  exceed  97.5%  confidence  and 
all  but  one  exceeds  99.5%  confidence,  it  is  concluded  that  the 
Chree  analysis  variations  are  also  not  the  effect  of  a  random 
process. 


CHAPTER  7 
OVERVIEW 

Suimnary  of  Results 

Wind  Models 

Two  models  for  solar  wind  velocity  stream  interactions 
have  been  developed.   These  models  describe  the  transformation 
of  the  velocity  wave  into  a  density  wave.  The  magnetic-field- 
free  model  appears  to  give  acceptable  results  when  the  peak- 
to-peak  velocity  amplitude  is  not  too  great  and  the  radial 
distance  is  also  not  too  great.   A  test  has  been  devised  to 
measure  the  propriety  of  the  model  in  a  given  situation.  The 
simple-field  model  solves  the  diffusion  problem  that  restricts 
the  field-free  model.   However,  it  permits  the  density  wave  to 
build  up  as  a  infinitesimally  thin  sheet. 

Both  models  are  adequate  for  a  coarse  description  of  the 
wave  exchange  process.   Nevertheless,  future  work  might  well  be 
directed  at  a  more  detailed  model  of  the  density  wave,  includ- 
ing the  effects  of  magnetic  and  thermal  pressure  which  will 
limit  the  compression  of  the  plasma  to  some  physically 
realistic  size. 

Jupiter- Solar  Correlation 

Employing  the  two  solar  wind  models,  the  1967-68  and  1974 
apparitions,  with  stable  solar  magnetic  sector  structures, 
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have  been  studied.   Evidence  of  a  correlation  between  the  non- 
lo  related  component  of  Jovian  decametric  emission  in  the  18  to 
27  MHz  range  and  the  solar  wind  flux  has  been  offered  for  both 
apparitions.   Positive  correlations  exceeding  95%  confidence 
occur  in  source  A  at  22,  26,  and  27  MHz,  and  source  C  at  18, 
22,  26,  and  27  MHz.   The  validity  of  the  source  A  22,  26,  and 
27  MHz  correlations  is  clouded  by  the  predominance  of  single 
zone  storms.   Evidence  of  an  anticorrelation  between  the  lo 
component  of  source  B  at  18,  22,  and  27  MHz  is  provided. 

The  solar  wind  velocity  data  which  serve  as  the  input  to 
the  simple-field  wind  model  used  for  the  1974  apparition  are 
limited  to  81  days  near  opposition.  A  re-examination  of  that 
apparition  when  more  complete  data  are  available  should  be 
undertaken. 

The  correlation  maximum  precedes  the  predicted  density 
wave  arrival  at  Jupiter.  This  raises  the  question  of  whether 
the  density  wave  itself  or  shocks,  caused  by  stream  collisions 
moving  ahead  of  the  wave,  is  the  exciting  agent.   Other  possi- 
bilities include  the  build-up  of  the  density  wave  1  to  3  days 
ahead  of  the  peak  of  the  velocity  wave,  and  the  acceleration  of 
the  solar  wind  with  radial  distance  as  predicted  by  so  many 
models.   The  possibility  of  an  accidental  correlation  also 
remains. 

Occurrence  spectral  analysis  shows  a  period  in  the 
decametric  emission  for  1974  near  the  26-day  solar  rotation 
period  giving  further  weight  to  the  reality  of  the  correlation. 
The  1967-68  spectral  data  are  inconclusive. 
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In  connection  with  the  wind  acceleration  the  author 
believes  that  the  method  of  reporting  solar  wind  velocities  for 
Mariner  2  and  Pioneers  X  and  XI  would  conceal  the  acceleration 
if  it  exists.   It  is  suggested  that  both  velocity  and  density 
or  flux  data  from  the  Pioneer  X  and  XI  projects  can  be  used 
with  a  mass  weighted  average  technique  to  resolve  this  problem. 

The  proposition  of  Vladimirsky  and  Levitsky  (1976)  that 
there  is  a  4°  wide  quiet  zone  centered  on  the  solar  equatorial 
plane  can  be  neither  confirmed  or  denied  since  Jupiter  was 
well  within  that  zone  during  both  apparitions.   Their  sugges- 
tion that  ejecta  or  shocks  from  solar  proton  events  produce 
Jovian  emissions  is  compatible  with  the  results  of  the  present 
study.   The  energy  density  of  such  events  exceeds  that  of  common 
flare  ejecta  by  several  orders  of  magnitude.   Though  it  is  not 
possible  to  predict  the  energy  density  of  the  stream  generated 
mass  density  waves  without  knowing  how  compressed  they  can  be- 
come, it  is  believed  that  the  proton  events  exceed  the  stream 
wave's  energy  density,  Thus  if  the  stream  generated  wave  can 
excite  decametric  emission  surely  proton  events  may  do  so. 

Spectral  Analysis 

A  persistent  spectral  line  at  frequency  index  17  (6.35d) 
is  seen  in  sources  A  and  C  during  both  apparitions  at  most 
frequencies.   This  and  other,  secondary,  features  can  be  ex- 
plained by  a  joint  lo-Europa  effect.   Miller  (1970)  has  pre- 
sented and  discounted  evidence  of  such  an  effect  by  showing 
activity  on  a  plot  of  lo  versus  Europa  locations  from  geocentric 
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superior  conjunction.   This  matter  definitely  bears  further  in- 
vestigation employing  more  than  the  two  apparition's  worth  of 
data  used  here. 

When  one  does  a  spectral  analysis  of  the  source  B 
activity  data  without  removing  the  lo  data  a  pronounced  activity 
spike  is  found  at  index  17  during  1974.   However,  the  same 
treatment  of  source  A,  with  lo,  yields  no  visible  line.   This 
intensifies  the  argument  for  a  joint  effect  since  the  index  17 
peak  cannot  be  explained  by  lo  alone. 

A  second  line  at  index  40  (2.7d)  has  been  discovered  in 
source  A  during  both  apparitions.   It  has  thus  far  defied 
identification. 

Statistics 

The  lengths  of  non-Io  related  storms  in  source  A  have 
been  shown  to  be  non-random  as  have  the  apparently  periodic 
fluctuations  in  the  solar-Jupiter  Chree  analysis  plots.   In 
source  A  at  18  MHz  these  Chree  analysis  peaks  are  seen,  though 
shifted  in  position,  even  though  there  seem  to  be  no  solar 
correlations  at  that  frequency.  These  four  peaks  correspond 
to  a  period  of  about  6 . 4d  and  may  be  connected  with  the  lo-Europa 
effect. 

Future  work  in  the  probability  distribution  area  may  be 
directed  at  removing  the  known  long-term  time  dependence  from 
the  activity  distribution.   In  principle,  this  might  be  done 
by  fitting  a  low  order  polynomial,  such  as  a  parabola,  to  the 
data.   The  inversion  of  this  function  can  then  be  used  to 
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amplify  the  activity  probability  near  the  extremes  of  the 
apparition  and  thus  remove  the  effect. 

Conclusions 

It  appears  that  below  the  strata  of  the  terrestrial, 
Jovian  and  loian  effects  on  the  decametric  emission  there  is 
a  set  of  discernable  secondary  and  tertiary  exciters  partici- 
pating in  the  "underlying"  processes.   The  current  examination 
has  revealed  evidence  of  three  such  exciters,  one  correlated 
with  solar  wind  interactions  with  the  Jovian  magnetosphere, 
one  identified  with  an  lo-Europa  effect,  and  the  third  at 
index  40  whose  possible  identify  is  not  known.   It  is  hoped 
that  further  studies  will  resolve  these  effects. 
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APPENDIX 
Data  Reduction  Programs 

Four  basic  computer  programs  form  the  nucleus  of  this 
study.   These  programs  are  written  in  FORTRAN  IV  and  should 
be  compatible  with  virtually  any  compiler  with  few,  if  any, 
modifications.   The  principal  functions  of  these  routines 
are  outlined  below.   The  programs  are  internally  documented 
with  niimerous  comment  statements.   Source  listings,  including 
these  comments,  follow  the  outline. 

S nop low.   The  solar  wind  interaction  models  are  gen- 
erated by  this  program.   The  input  data  consist  of  a  table 
of  daily  solar  wind  velocities  on  some  solar  radial  line  (the 
sun- Jupiter  radial  in  the  present  study).   This  table  is 
entered  in  card  form  using  a  10(15,13)  format  where  the  15 
represents  the  five  digit  Julian  date  and  the  13  represents 
the  corresponding  velocity.   In  this  study  the  table  has  been 
in  chronological  order  with  no  missing  dates  though  the 
program  will  handle  both  scrambed  and  missing  dates  (the 
latter  are  interpolated) .   The  output  is  both  printed  and 
punched. 

Jupzone .   This  program  is  the  interface  between  the 
"intermediate  deck"  and  the  two  Jupiter  data  reduction 
programs,  Jupcor3  and  JupsolS.   This  program  contains  the 
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filter  which  separates  and  counts  zones  for  a  given  source 
into  activity  and  observations  for  the  lo  and  non-Io  related 
components.   Jupzone  accepts  its  input  from  either  tape 
or  cards  and  generates  both  printed  and  punched  output. 

Jupcor3 .   Both  solar  wind  and  Jovian  data  are  the  input 
to  this  routine.   JupcorS  performs  the  autocorrelation  and 
Fourier  transform  procedure  on  these  data  and  prints  both 
values  and  plots. 

JupsolS.   Solar  wind  and  Jupiter  data  are  fed  to  this 
program  which  begins  by  performing  the  Chree  analysis  based 
on  a  user  specified  period.   The  solar  wind  Chree  analysis 
table  is  the  cross-correlated  with  each  of  the  Jovian  tables. 
The  routine  prints  values  and  plots  for  both  the  Chree 
analyses  and  the  cross-correlation  functions. 
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WliMU    MOOcL.     ScVL-r^AL    CPTIONS     A.-!L    AVAlLAoLii.     THE    U5=i.-i    i-'.AY    UjTAIM 
A    (  MAG.-icTIC)     FItLJ    Fi^dt    MUiJtL  t     A    3IMr>Lt£     FItLO    MUOt-L    JiiilCH    COM'J  tRVES 
KINsTlC    .-luMaNrUM,     oJT    ALLuWb    MA33    TU    J.JILD    UP    A.'i     I  ,MF  I.'H  TiTLY    THIM 
SHiii^TS.     Qifi    FINALLY.     A    KOjbL    .ViilCH    ASSUMinS    A.m    INVliKot    PuWEK    LA.V 
FOkCc     IS    CXtrtTED    JV    cACH    uROUr'    GIJ    Eve:JY     OTHcR    GHUUP. 


OIMt.NSICN    SUNST(  J,:3^5)  .  JL)ATA{  2i  300  )  , '.VNJ.^kMI  a,  160  J  .P.-<M:iX(  0,60) 
UiMtiNiiCN    B(30J  .WI^^Jl5.  JOG) 

COMMON    SUN3T,  J0ATA,JN0PSrt,LA5T.P.^MSX,KCJUNT,3,nlIiJ0.  I  STAR  i  .  I3TUP 
COMMJN    A  INC 

W.MDPRM    WILL    KEPriiiicNT    iMt    AHi^AY     OF    6    HJUH    GrtOUn'S    UF    PARTICLES 

WNDPrtMl  I  .N)=  OATc 

/;nDPR.-1(  ^,N)  =  VELOCITY     IN    KM/S 

iilNDPHMl3,M-  DISTANCE     IN    AU 

WNDPRM(4,N)=  WA3S  IN  ARTIRARY  UNITS 

WN0PR.M(3.N>=  SPAHc 

iKNt)PRM(6>N)=  bUNCH  tf 


INITIALIZE    WNOPHM 

00     10 0     1=1.160 
DO    100    J=l,6 
100    WNDPRMI J,!)=0 

PRINT     THE    SCALE    FOR    'PLOTS'. 

WRITE(6.9C1 )    AINC 

L=0 
AIC=1 

•  103  LQCP*  STEPS  GROUPS  THRJOGH  WNOPRM  AT  6  HOJR  INTERVALS 

NC=4*LASr«-160 
00  103  M=1,NC 

•101     LOOP'     t^EPOSITIONS    EACH    GriOUP    IN    WNOPKM    EACH    6    HOURS 

DU    101     1=1. 16C 

J=lol-I 

WNOPRM(  1, J)=SUNST( 1  ,1J  +  (0.25*L1 

IF(J.Ea.l)     GO     TO    102 

WNt)Pf(,v|(2.  J)=WNDPKM{^,(J-1)  J 

«lNDPft.M(3.  J)=««NUP/<M(3.(J-l))+(ai(NOPHM(2,(  J-U  )/b925.83) 

MN0PRM(4.  J}  =  WN0PR.v.i4,(J-I)) 

WNOPRMI  5.  J)=rtNOPHM(i.(  J-m 

iKNDPRMI  b.  J}=WNDPRM(6.  J-1  ) 

GO    TO     101 

EMIT    MUST    RECENT    GROUP    JF    PARTICLES    FHJM    THE    5JN 

102    K=AK 

AK.=  AK4-0.25 

IFIK-GT  .LAST)     GO    TO     110 

L=L«-1 

WNOPHMt  2.H  =SUNST(2.K) 

WNDPRMIJ. 11=0 

lKN0PRM(4.1i=l 

wNOPkMC  5.  1J=1 

WNDPRM(6>  1)=4. 

GO    TO    1 01 
110    MN0Pk.V4(  2.  11=0 

WNDPHM(  3,  H=0 

WNOPR.M(4.1)=0 

«kOPRM(b. 11=0 

tllN3PaM(6.1>=0 

L=L+1 
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FOrtTi<AN 
0040 


0042 
0043 
0044 
0045 
0046 
0047 
0048 
0049 
0050 
0051 
00  52 
0053 
00  54 
0055 
00  56 
0057 
0C58 


0059 


0060 
0061 
CO  62 


0063 
0064 
0C65 
0066 
0067 
0066 
0069 
0070 
0071 
0072 
00  73 
00  74 
0075 
0076 
0077 
0073 
0079 
0080 
OOdl 


LiiveL 
101 

c 

c 
c 

c 


21 
CUNTINUe 


UATr.    =    /aJJ. 


150 


105 
104 


FOR    FIfc-LD    FRc£    MODuL     IN^^dhiT:     ..U     TO     ll^i. 
FOR     SIMPLt    FitiUU     .MJOcL     lh,^e.ir:     00    TO     ll.j. 
G0~TQ'TT5 


THc    FULLOlalNG    KCUTlNa    CALCULAfiia    THE    FOHCH    £Xi£;^TcO    JM    HACH    oHaU? 

of    tACH    UTHE.H    GKuUh"    ANJ     AuJJSTS    TH;;     SUtJJcCT    GRUUP"b    V-LOCI TY 

ACCOkUI.Nii    TO    THt    .N^r    ACCdLciiAT  t C.N    Pi-fOJUCiiC 

THc     AbSUMPTIUN    .MAOi     IS     THAT     THt    FJH>_c     IS     Ai-J     INVdr^SL    POw^-R    LA;< 

FOhCc.     FuRTHck     ir    13    ASSUMuK    THAT     THE    rUKCc.    ACTS    UN    A    UulT    .-lASS 

OVtR    A    PcHIDU    OF    SI/C    hOjRS.     CONScCiUENfL Y     THC    Ut.T    CHA.Ni.-      IN     SPXCD 

IS    PROPORTIONAL     TO    A(X)**(-P). 

"A"    lo    THE    FukCE    CU:JSTANT 

"p"    IS     THfc    POmr.H    OF     Th£    AjSOM-iO    POu/tiR    LAJ    FORCc 

CONTINUE 

A=3 

P=2    •. 

DO    104     1=1,160 

IF{ *N0PHM(6tI ) .EQ.O)    GO     TO     104 

IF(  I.GT.l  J     lllNOPRM^^,I-l  J=aN0PrtM(2.I-l>*»^ 

F=0 

00    105    J=1.160 

If  (  *NL)PRMt6.J).EQ.O)    GO    TO    105 

IF(  I.cQ  .J)GO    TO    105 

X  =  l  *NOPl<;J|(3  ,1  )-wNDPRM(3.  J)  ) 

IF(X.ca.O)    GO    TO    105 

G=A»<<ABS    X)*»(-P)J*(Ad5    X»/X 

F=F+a 

CONTINUE 
CONTINUE 
•(NOPR.M(  2  .  1 60)  =.«llNL)RrtM(a  ,  1 60  ) +F 


lis    CONTINUE 


FOR    SIMPLE    FIELO    MQi>EL    LOCK    ONCE    A    DAY,     NOT    EVERY     SIX    HOURS. 

K=WN0PRM{1. 1) 
C=WNDPRM( 1 ,1) -K 
iFtC.Nc.O.r    OQ    TO    103 

11  =  1 

116  JI-1 

IX  =  H+1 

IF(li.GT.160)     GO    TO    114 

Z=l 

S=kN0PRM(2»II ) 

IF(S.cO.O)     GO    la     116 

117  IF(iiMJPRM(3  t  Jl)  .LI.  J.NaPHM{3,II)  )     GO    TO    118 
WNOPR-MC  3,  JI  >  =  WN0PRM(3>1I  ) 

Z=Z+1 
S=S+»(NOPRM(  2  ,  JI  ) 

118  JI=Ji+l 

IFt Jl  .LT.II)     GO    TO    117 
Jl=l 

119  lFt*NDPRM(3,JI)  .n£.*N3PR.>1(3,II)  )     GO    TO    120 
«(N0PRM(2»J1  )  =  3/Z 

120  JI  =  JIi-l 

IF( Jl .LE.II J    GO    TO     119 
GO    TO    116 


114    CONT/NUE 


169 


FURTHAN    IV    G    LcVEL       21 


0064 
00d5 


AT    THIS    PCINI     wNJPriM    HAS    0£i£N    RtOP.oA.-J  t  ZLTJ    FGk    ThlS    6    iiuJR    PbrilOD 
IT    MAY     NUw     BE    EXAMliNcLy    FO^     -vnATiIVtr<    P^'iOMtHT  I£S    A?.r.    •JF     INTtRcif 

CALL    PARM 


103    CONTINUE 
RETURN 

901  FORMAT!  •  •,//.  26X, 'SOLAR  ..  1  NO  PARMi-liiTER  "DS.NSITY"  VS.  JULIAN  DATE*. 
1  ,///,•  JULIAN*  ,1-iX,  •  INCREM-iJif  Ifj  bPACi£  ib'.F4.1,'  A^  r.<  J.N  JM  I  CAU  vJNl 
2TS.»./,'       OATc' ,oA. 251  'i,. •),•!•  ) 
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0001  SUb.iOuriNE    PARM 

C 

c 
c 

f  PAoL^r'^c-^Tif-^^'     ^'^V':^'^    "^^     •SOL.viN'.     iATKACTS    Fl^OM    THE    MOUi£L    WINO 

V,  HKtbuM     OP  I  IONS    Akc;  : 

C  1.     THt    FLUX    OtIiMSlTltS    AT    X    AU. 

Jr  rv.O     VALUco     MAY     jii     OoTAIMiiO. 


THtY     ARc     AT     X  =  A£riC*I     AJ  . 


^  Tmc    RcSULTS    Ai.iH     GTOrico    AS 

>  wINl)(2,  ICuUr4T)  =  j(  U       AMl) 

^  d.     THh    PArtTlCL;:.    0£:4SITIES 


c 


LOCATtO    At40    STU.iEO    AS    AJOVd. 


^  J.     THt    VtLGCITIiiS    AT    i^ADIAL 

J:  DISTA.-JCES    AS     %Pi.LiriEO 

^  IN     STATtM£NT5    IjO    AiNO     Ibl. 

c 

loll       '     o!;:^;:i{^;^^^^^[!^j^^^^:j^s}'''-'"^'-^^''«'^'^'^''^''^'^"^^*^'^"' 

•iO^*  DIMENSION  A(o,  loul 

OOot  cSS  |y|2!^'-"'*^^""'^'''"''-^''-^^'''«'''^''-'^'^-"^'''^'-I-^D'IS^'**^T.IST0P 

C  . 

<"^07  DATA  A,  ICOUNT/960*a.,0/ 

C 

C       LOOK  ONLY  UNCE  EACH  OAY  -  NOT  EVERY  SIX  HOURS 

C 

0008  K=WNDPkM( 1,1  J 

0009  C=*N3PHM(l  ►D-K 

OwlO  .      IF(C.N£.0.1  GO  TO  105 

C 

c 

C  INITIALIZE  a. 

C 

C 

0011  DO  10.3  1  =  1,30 

0012  103  d(I)=0. 
C 

c 

C       FOR  FLUX  I.MSckT:  NOTHING. 

C      FDR  DENSITY  INSERT;  GO  TO  Ufa. 


0013  "  ,15  CONTINUE  "*-*-**'****-**?i5S*^*,**»**?— 
C       FIND  FLUX  PAST  POINT  X  IN  STEPS  OF  AINC  IN  AJ 

0014  DO  100  1=1,30 

0015  Ai=I 

0016  X=A1»AINC 
C 

C  Th"plInI    OF**x!    '^^^'^    ELEMENTS    OF    WNDPrfM    TO    SEE    IF     THEY    HAVE    PASSED 

C  * 

0017  DO    106    K=l,4 

0018  IF(WNL)Pf.Ml3,K)  .GE.X)     GO    TO     107 

0019  GO    TO    1C6 

0020  107    6(n  =  ti(  IH-I 

0021  lOfa    CONTINUE 

c 

C  NOW    CHECK    THE    REMAINING    ELEMENTS    UF    liif.Mi>PRM    AND    A     (THE     'OLO*     VALUES 

C  OF    WNDPHM)     TO    SEE    IF    ANY    HAVE    PASSti>    THJUGH    TH:^    PLANc    OF     X.  '"    " 

0022  DO     101     J=l,  ISti 

0023  IF(  ACS,  J)  .LT.X.AND.'.Vn:)PRM{3iJ«-*)  .J=  .XI     GO     TO    102 

0024  GO    TO    101 

0025  102    3(IJ  =  d(  U+WN0PRM(4,.J  +  4» 
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FORTHAH 
0001 


IV     G    Li£Vc:L       dl 


OAr£     -     /ojl  4 


io/3i:/i>2 


oooa 

0003 
0004 

0005 
0006 

0007 


oooa 

0009 

coio 


0011 
0012 


0013 


0014 
0015 
0016 


0017 
0018 
0019 
0020 
0021 


0022 
0023 
0024 
0025 


bUB.<UUriN£  PARM 


THIS     SUBRUL>iINc.     CALLcO    dV     'SCiL^IN*.     cATHACTS    rWJf.l     THE     MOOcL    WIND 
PARAKcTckS     THGSc    rfrtlCH    ARc    Of-'     iNTcRtSr    FUR    LAfc:.^    DISPLAY. 
PRtScNT    OPflUNS    are: 

1.  THt    FLUX    Di-riSlTItS    AT    X    AU. 
rv,0    VALUco     MAY     jc     OyrAlMr;{J. 
THcY    ARt     AT     X=AINC*I     AJ. 
Ttif.    RcGULTS    Ar<c     STORciJ    AS 

wiNij(2,  icuurir)=j(  I)      ANO 

Vvi.NO(3.iCUU,-4T)-J<  I)  . 

2.  THt    PAftTlCLt;    Oc^NlSlTIES 
LuCATtO    AriO    jfUiJtU    A3    AJOVc. 

J.     THt    VtLUCITIcL;    AT    RADIAL 
Oi  STANCES    AS    SPECIF- 1  cfO 
IN     STATEMENTS     laO    Anl)     Ibl. 


UI.MENiIGN    SUNSrC  J,3abJ  ,  JDATA(2,  300  >  .  //.MOPft;-»{  6,  loO)  .PRMSX(  6,60) 
DIMENSION    dliQ)  ,iilHOli>,A00) 
DIMENSION    A (6. 160) 

CUMMUN     SUNSI,  JDATA.  WNDPHM,LASTtFHMSX,KCJUNT.d.rfINL),I5TAf;r,  I  STOP 
COMMON    AINC 

DATA    A,     ICOUNr/g60*0..0/ 

LOOK    ONLY    QNCe    EACH    OAY    -    NOT    EVERY    SIX    HOURS 

K=WN0PkM(1 rl) 
C=«NDPRMC1  ..IJ  -K 
IFCC.NE.O.)    go    TO    105 

INITIALISE    B. 


00    103    1=1.30 
D(I}=0. 


FOR  FLUX   insert:   nothing. 

FOR    DENSITY     INSERT:    GO    TO    116. 


107 
lOb 


102 


m^iff  trf^  *******  ff'r****!^^**^!^'^**:^ 
CONTINUE 
FIND    FL'X    PAST    POINT    X     IN    STEPS    OF     AINC     IN    AU 

DO    100     I=U30 
AI  =  I 
X=A1»AINC 

CHECK.    THE    FIRST    FOUR    ELEMENTS    OF    *NOPrJM    TO    SEE    IF    THEY    HAVE   PASSED 
THE    PLANE    OF     X. 

DO     106    K=1.4 

IF{*Ni)PRMl3,K)  .GE.XJ    GO    TO    107 

G0    TC     1 06 

d(I)  =  B(I)4-l 

CONTINUE 

NO*    CHECK    THE    REMAINING    ELEMENTS    OF     tfNDPRM    AND    A     (THE     'ULO'     VALUES 
OF    WNDPHHJ     TO    SEE    IF    ANY    HAVE    PASSED    THJUGH    TH:.    PLANc    3F     X. 

00     101     J=l« 156 

IF<  A(  3.  J)  .LT.X.AND.WN0PRM{3iJ*4)  .lit.X)    OU    TO    102 

Go    TO    101 

ti(I)  =  D(  l)'*'MNDPRM(4..J-f4) 
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0001 


0002 
0003 
0004 

0005 
0006 

0007 

0008 
0009 
OCIO 
0011 
0012 
0013 
0014 
0015 
0016 
0017 
0013 
0019 
00  2C 


0022 
0C23 
0024 
0025 

0026 

0027 
0028 
0029 


iV     &    LtVEL        21 

SUc.^OUTlNE    PLUra 


0A1£    =    /oJl-4- 


THIS     oUiRuUTlNE     15    CALLED    dY     'FLUX*.     IT    PhJOUC^S    A    DISPLAY    OF 
JULIAN    OAlc     VS.     Sl.LA«    W  I  Nb     UtlNSlTY     VS.     iiAUlAL    OlS.TACt.     IN     A.U.     FROM 
I  He.     &UN  • 

JIMdNSlCN    SU.MST(3,3.i5)  .  JOAT  A(  2  .  iOO  )  .  W,  JJPr<M(  o,  loC  )  .;^r<MSX(  o  ,  60) 
UlMtNSiLN    B(  JO).WlNO(ij,:jOC) 
DIM=NSII_N     lOUTdl  J.KUuKyO)  ,     LULir(Jj) 

COMMJN     SUNST.JDATA.wNuPWM.LASr.PHMSX.KCJUNT.B.WlNO.IirA.Vr.  ISTJP 
COMMON     AINC  .1-1-*"    i.iaiJf 


Data    Ijut/? 


DO     100     1=1,83,3 

J-(li-2)/3 

IA-(o( J )+.5)/10. 

A  =  Ia 

JA=(t-(  J  )  +  .5)-10.*A 

XUUr( I) =IuUT(l J 

Ir(lA.dG.C)     KCUT(  H-l)=Iuur(  1  ) 

1F(  lA  .tO.O.AND.JA.EO.O     KfjU  F  t  I +2)  =  1  l)UT(  I) 

IF{  lA  .Nc.C.ANO.JA.cO.O  )     KauT{  l+il)=lOur(  1  1  J 

00     10  1     K=l,y 

IF(lA.cO.K)     KJUTC  lH)=IJUT(Ki-l  J 

101  IF(JA.EG.K)     KUUri  I+.a)  =  iauTiK<-l  J 
100    CONTiNUE 

MAP 'jK  I  NO    DENSITY. 

W KITE (6, 900 )     WNDPrfM( 1,  l),(KOUT(  I). 1  =  1,90) 

PUNCH    OUT    WIND    MAP. 

00  102  1=1.30 

102  LOUr( I ) =5(1 ) 

IOAT£=tfN0PRM( 1,1) 

«li^lTd(7,901  )  lOAIt,  (LOUrtl),I  =  l,30) 

HETUkN 

900  FOKMAK  •  •  ,F6.0,3X,90(A1)) 

901  FOrtMATt I5,3CI2) 
END 
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FORTRAN  IV 
0001 


0002 
0003 
000* 

0005 
0006 
0007 

oooa 

0009 
0010 

ocu 

0012 
0013 
CO  14 
00  IS 
0016 
0017 
0018 


G  LEVEL 


0019 
0020 


0021 

0022 
0023 
0024 


21  OUTPUT 

subrtouriNE  auTPor 


DATE    =    7oJl'f 


PUNCH    UUT     "WIND" 


DIMcNSICN    SUNbTl  J,3ii3)  ,  JCATAC  2,300)  ,  .il.JiJPfiMC  6  .  1  bO)  .PrSMSXl  o  .60) 
OlMcNSI  L,N    8(30)  .  WlNOCb  ,  JCO) 
OIMcNSIUN    Iyv(2.l0) 

CGM.IGi-J     SUNSr,  JD  ATA,  WNDPi-<M.  LAST.  Pi?M6X.<CaUNr,b,.vl;-IO.  ISTARf.ISrjP 

COMMON    A INC 

WfiITE(o,902) 

DO     1    N=2.5 
J  =  C 
K  =  l 

00  I 
IF(K 
K  =  I-J 

1  dr(  I. 
I W  (  2  . 
IF(K. 
IF(I. 
CUNT  I 


I-l .KCCUNT 
cO.lO)     J=J+10 

<)=.VIND(  1.1) 

K)=WINO(N, I) 

eu.lC)     riKi  TECy.-iCO)  (  (IW(L.M)  .L  =  l  .2)  .H=l,  10) 

£0.k;CUUNT)     wKl  Tut  7,500)  {  (  la<L,M)  .L-l  .2)  .M  =  l  ,K) 

Nue 


DO    2     I-l, KCCUNT 

WRITE  (6 .901)     MINJ(I.I).     UINL)(4,I)< 


RETURN 


««IN0(5.IJ 


900 
901 
902 


FORMAT! 10(15, 13)) 
FOftMATi'  •,3Fa.O) 
FORMAT!  •!•) 
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ei-J'JF^CMF     -FORT  -SYf;PqiNT 

FORTPAN    IV    G    LEVEL       21  MAIN  DATE     =     76314  16/02/17 

c 
c 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

Of>01  0!M(-N=.ICN    j;CNP?(l.ft),     KOATA  (  )?,2C0':') 

0002.  CCM^CN     INDEX,     .S-OAT/i,     JZOM''S»     I  FLUX 

c 

0003  CALL     INPUT 

c 

0004  DO     999     INDEX=1,4 
C 

0005  CALL    ZCKES 
ft0"6  899    CONTINUE 
0007  STOP 
000  3                                         FND 
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PCRTPAN     !V    G    LFVFL       ?A  INPUT  OATr    =     7631ft  l'S/C2/17 

000  1  SUiBPCUTINF     INt'UT 

C. 

C 
000  2  nlM=-.N^.  rCN    J?CNt"?(4,4)  ,     KDAT  C  (  1  ■=.  ,  PTC  0  ) 

C003  C0-»V(:N     IN(>EX,     KPATs,     J^^MTS.     ii-Llix 

c  .  ■ 

C  PEAO     IN     ALL     DATA 

C  V 

0014  c;FAP(5.'--rf.  )  (  (  JjrN=^S(  !.J),I  =  l,a).J  =  1.4) 

OQ05  Pf:Ar,(?,60?  )      IMINtlMAX 

c 
c 

0006  OC    21    J-1.20r0 

0007  1    R'^Af^C.^C  1  ,FNC -^1  )  (KOaTA{  I  .  J).  I-l  .15)  .IFLUX 
O'^OS  l"(  KOiTACa.  J)  .LT.  IMIN)     GO    TO     1 

OOOC  I^^I  KHATAI?,  J  J.GT.IMAX  )     GC    TO    21 

00 IC  GT     TO    23 

0011  21     KOATA(l.J)=0 

00  12  KnATA(5.J)=0 

C 

c 
c 

001  3  D'^    l*"    1  =  1  .500 

0014  \r\    WPITF(6.6C?)     KDATAt?!.!  ),KPATA(  12.1  )  ,KOATA<  13.1)  ,K0ATA(14,I)  , 

2KDATA(15.I) 

0015  eC5    P0?5MAT(  1X.5I7) 
C 

0016  M^TUt'N 
C 

C       ALL  DATA  MATf^lCES  APE  NOW  LOADED 
C 

0017  23    CONTINUE 
C 

C 

001  e  6CO  c^^^r^•lT{  ie<  1 K.  I  2) ) 

rOlQ  eCl    FmMAT(  IX,  I2.1X.  312.  I5,f-I2il£X,2!2.  12X,2I2.3X,I7)  5D    Z0^ 

0020  ':02    C0P^'lA  r<  I?,l  X.I5> 

0021  cNT 
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FQCITPAN     IV     G    LFVFL         " 

'*'''-"  S  DATt     =     7f3!4  16/^^2/17 

C 

c 


000  5 
one 
0007 
0008 

000  g 


001  1 


CO  13 


001  £ 


0022 


0023 


C  THIS     PCUTINF    CrUNT^     JUO'rro      t-im-c     ^r-,...r.^ 

C  •>J':N-rj    Fcp    c-ACH    r^Y.       ray-     ,^-:'-'":^,.S'^V    ;7n^    ^■•'^    Ofif-f-v-T    -nr,-     TO    AMD 

c     CUTPUT  AP-  Pr-vro^o     ^  r  v.-.r- vATicrj-  days.  p.:,-inted  ano  punch^'d 


„«i?  JM5nr?rjN;pc:i  +  7 


SOUf^Cli      INITIALIZATION    STOUFMC:. 

J35nci-J/rNFS(  1  ,  IMD'^X) 
J'^JS'^C!  =J7rNFS(?,  lUtjc-X) 

t-i"?:  =J:fCNF.s(  3,  iNn'--x) 

JHSPC2  =  Jt!Sf-CH-72 


Il'ja^ILVl  <-72 


NJZS-0 


00X1  NJZ-0 


T  J  Z  5=0 


POlt  IJ/=r 

22 II  JUL=6 

«"»e  JNJB1=0 

2?'^  ir3UNT=0 

00|0  ICAKC=1 

S.,fi  IF-:=KnATA(l.l> 


ILINES^O 

IP    THEF^F    AFE    NC    MOPF    SOU^CPS    SPtC.IF  IFD— fxI  T. 

IF(JZ0NFS{l,lND-X).F0.o.ANO.JZ0NeS{2.INPFx).E0.0)     GO     TO    399 


c 

c 

c 
c 
c 

c 

c 

C       WRITE  FIRST  PACE  HEADINGS. 
cols  ^       W«lT^{t:t;lJ  ''''■''''     •'^^'^^^-  ^^^».  INZI.  IFF 

c 

C       PUNCH  DATA  eLCCK  HSADER 

°^^*  ^  «'f^ITF(7.6C4)(J2CNFS(I.lNDEX),I  =  1.4) 

C 

C  FND    CF     SOURCE     INIT I AL IZAT lOM. 

C 

c 
c 

C  SEGIN    DAILY    COMDUTATIDNS. 

C 
loll  »     IF=KOATA(I.ICAR0) 

oSfl  JY=KDATA{2.ICAPD) 
Oolo  ^'  =  KnATA(3.1CAFD) 

0^,Y  IO=KDArA(A.ICAR0) 

00,1  J'JLIAN=KDATA(c;.icACD) 

loll  JNC.9=KDATA(7,ICAFn) 

003*  J3ACl=KDATA(o.ICASO) 

OO^fi  JNACl  =  KrATA<C,ICA''D) 

■^^  leOB=KnATA(  10,  ICAPD) 
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f>^7 


I>n5'-Kr/.TA(  1 1 ,  icX'^o) 


CMS  ICA'"-t;  =  ICA'^D+l 


c 

C  DATA    rCF    riexT    CAY'S    VALUr.S     A^l     STOnO    PENDING    OUTPUT    OF    CURP^MT 

C  OAY'S    VA1.U=:S. 

c 
c 

C  IF    THIS     IS     r\^=     LAST    OAfft    CAOO    —    P':TUPM. 

C 

c 
C03-5  iF(  ir.F.-ii.->  .AND.  JULIAN. ira.r)   go   tc   409 

c 

c 

C  I"^    THIS     TS     ThF    "^IPST    PAY     FC     THIS     SOU'^xCE  .     TMfc'RE     IS    NOTHING    TO 

C  OUTPUT    YFT    —    r-HGiN    COMPU  r  AT  IQNS. 

c 
c 

Of)ACi                                             IF(  JUL.FQ.O  )    GO    TO     7 
C 
C 
C  IF    This     13    AM    AOPITION'AL    CARP    FCR     ThF     CUf-KENT    OAY CONTINUE 

C  ccmputatign. 

c 

c 

0041  I^t JULIAN. rC. JUL)     GQ    TG     18 

c 
c 

C  IF  CAPD  CUT  CF  CROER.  PP I  NT  MESSAGE... 

C 

0P42  IF  (  JUL  IA!^.LT.JUL)  GO  TO  20 

00  43  Iri:)UNT=!CCUNT+l 

0044  GO  TO  19 

c 

0045  20  PRINT  903 

0C4  6  PniNT6»M  , IF, lY, IM,  ID . JULI AN . JQC3. JN09. J3AC1 . JNACl .180  3. I  NOB 

C 

0047  le  !F( JNOP.EO. JNOEl )  GC  TO  1 

0048  IF{  JCOB.!=:O.JNCBl  )  JBr.3=JB0B  +  l 
004O                  GO  Tr  7 

c 
c 

C       IF  FOPTV  LI^b"S  CF  OUTPUT  PPINTEO  GD  ON  TO  NEXT  PAGE,  LA-3EL  A?»0  GO. 
C 

c 

0050  19    ILINES=ILINF£+1 

0051  IF(  ILINFS.EC.4.1  )    GO    TC    40 

0052  GO    TO    22 

C 

c 

0r53  40    WPITEI^.'^P*;) 

C054  WPI  TF(6,fc'>3) 

C055  ILIM'^S=1 
C 

c 

C       PRINT  AND  PUNCH  OUTPUT....... 

c 
c 

cose  22  IF(NJ2S.GT.oi  GO  To  25 

0057  /(PITE(e,602)  JUL,  IJ  2S  .  IJZ.NJZS  .NJJ 

OOSe  GO  TC  21 

c 

c 

0059  25  <ilKITE(6,f'C2)  JLL  ,  I  J  ZS  ,1  J  ?  .NJZS  .NJZ 

0060  WRITE(7,e07)  JUL , I  J Z5 . 1 JZ .N J2S .NJZ . IF 

C 

c 

C       IF  THERE  IS  A  DAY  MISSING  FROM  THE  INTEPMEDIATC  DECK SHOW  IT  AS 

C       A  OAY  WITH  KG  CF3EP VAT ICnS . 

c 

c 

0051  21     1F( JULIAN. EC. (JUL*1 ))     GO    TO    23 

0062  JUL=JUL+1 

C 
r 

C  IF    FORTY    LIKES    OF     OUTPUT    PRINTED    GO    ON    TO    NEXT    PAGE,     LA13EL    AND    -50. 

C 

c 
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i-.VfrU        ^1  ZCNSS                                           nATF     "     76311                               16/02/17 

OOf'2  42     ILIN!^S=rL  INf  S4-1 

i??^^  r-^C  ILlNrS.L-  C.'il  >     GC     TC    4  1                                                                         . 

Of>r5  rjo    to    a3 

C 
C 

C067  W^I  TF  (  f' ,,  (5^  ) 

<>''6e  !LINFS=1 

c 

c 

C  P'-^INT    ANP.    PUMCH    MIS^.  rrjG    OAV. 

c 

;'';^^  ''^    voiTf-(f  .ff.e  )juL 

^23?  V(mTF(7,-;U')      JUL,       IF 

0071  GC    TC    ^l 

c 

CO 7  2  J>3    CHNTINUP 
C 

c 

c 

c 

c 
c 

C073  Ijr=o 

'^'^75  IJ'S^O 

0077  7    M=i 

0"7S  vi=^ 

^J^T^  IFf  Jf?np.r,T.jrG5)     JNCr,=  jN0G+7? 

^2S?  IF{  iBOP.GT.INOe)     INC'3    =     I NCB    «■     72 

"^i  IP(  JRACl  .GT.JNACn     JNACl  =  JNACl  +  72 

""''S  JQAC2  =  JDtri +72 

''O^-  JNAC2  =  JNACI  ^72 

00P4  Jf?0=JBCi< 

ooee  GD   TC  ^ 

'^^''^  4     JS^C(N)  =  J609 

'''?''l  t    JT^W=jr3CBH 

0002  N=N<-1 

92^:  I!^(  JBCO.LP.JKCP)    GO    TO    3 

0"^*  jN'3Bi=jNna 

fO^S  8    nci    9    J  =  N,59 

0096  9  JSPC{J)=0 
C 

C 

C  IF  NO  SOURCE  ZCNFS  CBSERVFD  M=0 .  GC  READ  ANOTHER  DATA  CAPO 

C 

0097  IF(M.EO.O)  GC  TO  17 

C 

c 

C  CALCULATE  JUP  ZONES  COPPESPONO ING  TO  10  POUNDAKIES  FOR  T  1*0  CYCLES 
C 

0098  Ii30  =  IB0e 

2'>2^  K=2548,e/5?5.495 

OJ29  Jnzi  =  Jeci+x*(iRzi-ien> 

">•*  J!NZI  =  jar.  +  x»nNZl-I  E0» 

21!?|  Jl9?2  =  J°CtX*(lEZ2-IeC) 

"^03  JINZ2=JSC+X*(INZ2-IQ0) 

c 

r  TE^'^nrl^r,'^?";    ^^^    ^'^TICES.CNE    WITH    ALL  '  OBSERVFD    IC    RELATED    ZCN^S 

C  THE     CTHER     WITH    NON-IO     ZONES                                                                                                «-<-i-j 

C 

2^'"»  DO     12    J=1.99 

"J^f  T-^f  JSRC(J).EC.O)    GO    TO    11 

oJpt  IF(JSPC(J).Gr.JIP7i.ANP.JSP.C(J).LE.JINZl)    GO    TO    10 

oiee  Nspcfj>=ji^c{:^'''"'*'"°-'^''''-"-'---'''''^*  ^-  ^^  '° 

Cl'^9  1S=C<J)=0 
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rO«TRAN     IV    r,     Lr  Vv  L        ?.\  (^TN'TS                                            OAfc     =     /f^Jl*                               1<=-/C2/17 

Oil''  GT    Tf-     1  2 

0111  10     !G^C(  J)=J'^?C(  J  ) 

0112  NS7C<J)=n 
C  U  3  GT     f     1  2 
«^114  1  1     r  -•>rAJ)=C 
0115  N'-^xC(J)=0 

one  12  coNriNUE 

r 

c 

C  COUNT    STCPVING    ANO    NtJN-STOr^M  I(NG     10    ANO    NCM-IC     lOHES 

c 
c 

0117  on     irr     J  =  l  ,^0 

0118  IF(  IS=C(  J)  .FiCO)     GO    TC     la 

0115  rj^-ij?+i 

0120  tP(  I-.(-C(J)  .G*::.  J-'ACl  .AND.  IS>'C(J  )  .LH.JNACl  )     GO     TO     13 

0121  IF  (  l;iFC(  J)  .G5.  Jr;AC2  .AMn.  ISkC(  J)  .Lf:.JMiC2)     GC    TO     13 
ryiPZ  GO    TO     1» 

C123  13     IJ^?=IJ?S+! 

0124  1'*    IPCM^r^cc  J)  .EQ.O)    GO    TO    16 

0125  NJ^=NJ?+1 

0126  IF( -IS'^CtJ)  .Gf^.  JTACl  .ANr-.N>:  =  C(  J)  .Lf^.JNACl  )     GO    TO     15 

0127  I-^CNSRCCJ)  .G?.  JKAC2.Ar;n.NSf;C(  J)  .LE.  JNAC2)     GO    TO    15 

0128  GO    TP     16 

0129 

C130 

C131 

0132 

C13  3 
013A 

ftl3£  fCl    Fn:?MAT    IIX,     12.    IX,    312.     15.    612) 

0136  602    FOPMAT     (•  '.     I'S.     3X,     4110) 

0137  603    =nPMAT(////,lCX,«  JULIAN' ,°X.  •  10    PFLAT<Er>     ZCN!;S  •  . 'iX  ,  •NGN-IO     ZONrSV./ 

IX,/,  1  IX.'OATF*  .'^X,  •  ACTI  v;^>  ,5X.'  TOT/>L*  .AX  ,  ■  ACT  I  V£«  ,  rX  ,•  TOTAL'  ,// ) 

013e  604    =0RMATC'4(  I2,2X  )  ) 

0139  605    F07MAT(  •  1  «.////, oSX.'JUP^CNE' ,//,ir>X. 'SCU^^C?    L1.=S    FPCM    /'ir-JF  •  ,  I  3,  ' 

ITO"  .  13.  ■  .  '  .  3X.  '  IC     ZCNFS    APF    FROM', 13.'     TC  .  I  3  ••  .  •  ,/,  1  OX.  '  F;:.FOUJ?NCY 
2'. 13.'     VH2.        ZG'JFS     AHic    triviT    DPGKE-S    OF    ci'L    ■tlDE.') 

OlAO  606    FOPMATC     '  ,cx  ,I';.3X.-i(  cx.'O'  )) 

0141  eC7    "^O'-MATf         •  ,  15  .3X.4I  It  .  2X.  12) 

0142  'SO'?    ="RMAT('         •  ,  I5.3X,4(OX,  •(■>•  ),2X,  12.  •  NC    DATA    FCP     THIS    DATE.') 

0143  «i09    FnOMaT(  •  I  •  .////) 

014  4  90 ?  FORMAT  (//.5X.'DATA  CAPO  OUT  OF  OROhR  F0LLPV3 • ./ ) 

01*5  FNr) 


1  5 

NJ7<;=NJ£S+1 

!<=• 

C  JNTINOe 

17 

JUL=JULIAN 
GO     TC     1 

c 

*. 

c 

39  fl 

CALL    EXIT 

c 

40  -i 

K'-TUf<N 

c 

c 
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FORTflAN    IV    G    UtVEL       21 


DArt    =    76  31'% 


16/0o/4fi 


0001 
0002 


>»»»>>»>»>»»>»        JUPCUKJ3       <««<<<<<«<<««<<« 


THIS    R'JUTlNt    CALCULATES    AUTUCURRELATl  LlN    AMJ    POvldU    Sr>oiCTilA    OF 
JUPiTEK     ACTIVITY     AHJ    J  JoiiK  VAT  IJNS .     IT    OOti     FHc.    oAM;;    F  Clri     THE    SOLA.^ 
IIIWU    .MULicL    OcNSIIY* 

DIMtNoI  CM    AdP0CH(3.  JOu)  .     rt  I  i>iD  (  2  .  3C  0  J  .     IJAIA(b.30O).     JZuNtSt*) 
CaciKQN     AiiPGCH.  lat^.  Iiiriii,     iVi.'^Dt     lOATA.    JZONuSf     XrMcCJl 


CALL     INPUTl 


00C4 
0005 
0006 
0007 


ooca 

0009 


0010 
0011 
0012 
0013 
0014 


0015 
0016 
0017 
0018 
0019 
0020 
0021 


0022 
0C23 
0024 
0025 
0026 


0027 
002d 
0029 
0030 
0031 
0032 


0033 
0034 


OG    aOO     1-1,300 
AtPuCHI  1,  I)=.VIND{  I.I) 
A£PaCH(2,I)  -WIND  (2.  I  J 
A£PaCri( J,Ii--l 


CALCULATE    i«lND    D^INSITY    AUI GCOKkELAT I ONS    AND    iPtCTRA. 


201 
202 


403 
402 


wRITc(&,9001 
CAUL    CQHcL 


bcGlN    CALCULATIONS    FCR    THE    JUPITER    DATA. 


CALL     INPur2 
00    4C  1     1=1.300 
AEPOCH( l.I)=IOATACL.IJ 
AEPUCHl 2.1J=-1 
AEPOCrtl3. 1J=-1 


CALCULATE     10    ACTIVITY     iNDICtS. 


00    402    1  =  1.300 

C=It)ATA(l>I) 

A=I(}ATA(2.t  } 

B=iOATA(i>I ) 

IF(  b.  £ii.O)d=-l 

0=A/ii 

IFia.dU  .-1 JD=-1 


CALCULATE    NtJN-ta    ACTIVITY    INOICcS. 


A=10ATA(4.I 1 
d=IDATA(5.I) 
IF<3.£vi.0)a=-l 
£=A/'B 

IF(  o»Ea.-l}£=-J 


LOAD  INDICES  INTO  AcPJCH. 


DO    40  2    J-  1  .  30  0 

IFIAEPUCHd  .Jt.EU.C)     GU    TO    403 

GJ    TO     402 
AEP0CM(2,JI=D 
AdPOCHI 3.Ji=e 
CUNTINJE 

CALCULATE    ACTIVITY    AUT UCUKKEUAT I QNS    AND    SPECTRA 

WRITE(6.-«C1  }(  JZONESl  1}  f  1  =  1.4).      IFREQl 
CALL    CJkEL 
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FCrtTRAN    iv     t.    LhVfcL 
C 


MAIN 


0035 
0035 


iJATE    =    /oJl-f  .  16/03/48 


DO    200     1=1  , JOO 

0030  ^^        A-IJArA<L.I> 

"^'•^^  200    AcHi}C.H(  J,  I)=A 


C 


C  CALCULATE    a3ScRVAriCN    mUTOCo^KELATI  CNS    «NC    SPECf^^A. 

004?  C^LL'^cSRiL^''-'^^''"^' "•'=!'*»  '      '^«^^^ 


"^O^S  GO    TO    201 

*  90X     FuH.MAr(  •  1  •  ./y /•      ,.    v     .    ,.,,•-,■,.      ,, 

r?;:'f\-     ••^^••'^     '--^^^^    ^-^UM-J^^li-^e^L"^  ;^0M    ZU.C..13.. 

3D£.'J  >JUiM^.r  tc  ,ij,    .   zj.Ni^^  Ar<c  ^  I  vt  0c6R-iSS  OF  C.-lLJl 

C 

0046  STOP 
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FORTRAN     !V    G    LclV^L       21  ^,,,. 

°''°*  SUBROUTINE     INPUT  1 

C 
I  T»,S    R.or.ME    .«.s    1.    ,«    SU..„    „„„    ,,,,    ,,„    ,,^,,^     ,^     ,^    ^,^_^_ 

DO     1     1=1,291.10 
°<^*°  1     CONTINUE 

c 

^  IF    -INO     IS    CUT    OF    DATA    -    DEFINE     IBEG    AND     I  END . 

tin  ^    iBEG=wlNO(l.l| 

"^'^  ^  ItNO^WlNOd  ,K-li 

^"^•^  RETUHN 

C 
0014  ^    901     FUftkAT(lO(F5.0.F3.C)) 


0004 
000b 
0006 
0007 
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FORrRAN  IV 
0001 


0002 
OOOJ 


OATt  =  /6J1* 


0004 
OOOb 


0007 

oooa 

0009 


0010 


0011 


0012 
0013 
0C14 
0015 
0016 
0017 


0013 
0019 


0020 
0021 
0022 


0023 


G  LtVEU   21  MAIN 

SOBkOUTINc.  INPUT2 

c 

C  THIS    HOuriNt    LOAOS     THL     JUPITER    OATA    BLOCKS. 

c 

C  LOAO    JUPITER    DATA    BLOCK. 

C 

1     Rc:AU(5.'J01>     I  OPT 

;    RcAL»<6,9Ci:.£N0  =  <J)(J£uncS(I)  .1=1  ,4i 

SKIP    BLANK    OR     IGPT    CAROS. 

IF(  J2uN£S(2).i£0.0.A,NJ.JZ0NcS(o).ea.0)    GJ    TO   2 

1    1=0 
•    1=1+1 

riEADl5.903,.£NO=6)(I0ATA{  J,I).J=1,5),IFREQ 

IF  IpATA  FULL  KcTORN 

IF( I.EQ.300)  GO  TO  6 

IF  OUT  OF  INPUT  OATA  —  ENTER  ZEROES  FOR  NEXT  DAY. 

IFdDATAd.D.EQ.C)  GO  TO  6 

IFREQl=IFR£Q 

GO    TO    4 

00    7    K= 1,5 

IOATA(I,K)=0 

iaES;:IOATA(  l,li 

IENJ=IDATAM,I-1) 

RETURN 

CALL    EXIT 

FORMATI II ) 

rOSMAT(4C  12. 2X)) 

FORMAT! 2X. 15, 3X, 41 10,2X,I2) 

END 


io/C8//ra 


901 
902 
903 
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FURTaAN  IV 
0001 


0002 
0003 


0004 
0005 


0006 
000  7 
OOOd 
0009 
0010 
0011 
0012 
0013 
0014 


0015 
0016 
0017 
0018 
0019 
C020 
0021 
0022 
0023 
0024 
0025 


0026 


0027 


0028 
0029 
0030 


0031 
0032 


0033 
0034 


21  MAIN  DATE  =    76314 

suoaauTiNE  cuHEc 

THIS    RCuriNE    CALCULATES     THE    AUTUCuk^ELAT  I  CN    CUii.-r-IC  IciN  f  S  . 

L>IM£NSiaN    AEPCCHCJi.JOO)  ,     C0RCaH(J,55) 
CUMMuN    AcPQCH.IdiiC,  IcNU 

THE    MAXIMUM    LAG     ,»ILL    Dc    55    DAYS. 


JLIM=0 
ILIM=b5 


16/0a/48 


INITIALIZE  VALUES. 


00    304     1-2,3 

00    30  4     K=;l,  ILIM 

SUM1=0 

SUM2=0 

SQK1=0 

S(i32=0 

CRas=o 

JDIV=0 
JLlM=IENU-IbEG-K.+2 


CALCULATE    THE    COEFFICIENTS. 


DO    303     J=U  JLIM 

L=J+H-l 

ALPHA^AEPOCrid.Ji 

dETA=AEPUCH( I.L) 

IF( ALPHA. LT.O, art. dETA.LT.O)     GU    TO    303 

JOIV=JOlV+l 

SUMl=SUMl+ALPHA 

i.UM2=5u.vi2i-3ETA 

SGkl=S0Rl+ALPHA*42 

S0fi2=SQk2+3ErA**2 

CftOS=CftQS<-ALPHA»tJETA 


CONTINUE 

AVOID    DIVISION    oY    ZErtu. 

If ( JDIV.EQ.O)     GO    TU    3G2 


A=CR0S-SUM1 *SUM2/JDI V 
E=SOft.Tl  SQR1-SIjM1**2/JD1V) 
C  =  SC»RT(  SaH2-5UM2»«2/JDIVJ 


AVOID    DIVISION    aY    ZEriU. 
IFIB.ECl  .0  .OH.C.EQ.O)    GO    TO    302 
COEF=A/(B*CJ 


GO    TO    305 
C0EF=-99 


305    CONTINUE 


185 

C  LOAD     VALUti;>     INTO    CUWCUH  . 

C 

00  36  COKCQr( 1 ,K) =K-1 

0037  JOi»  CGHCJF(  1  ,K>=CUt£F 
C 

c  ■ 

C  QEGIN    PLOrriNG    iidfJUtNCd. 

c 

003a  |_=2 

0039  M=3 

00*0  N=55 

0041  IHISlf=l 
C 

^'^''^  CALL  3PCTkM(CJ«Cu>-,c,.v|.IHM) 

0044  L=3 

0047  HcruSN 

0045  tNo 
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-ORTrJAN    IV 


0002 
0003 


3004 
000b 
0006 


0007 

oooa 


0009 
00  10 
00  11 
0012 
0013 
0014. 
0015 
0016 
0017 


0018 
0019 
0020 
0021 
0022 


0023 
0024 
0025 
0026 


;U       21  :'.AI.M 

SUaHOUTINE    SPCIHM    <K. LI, HI. NX) 


JATS    =    7o31* 


THIS     K!O.JTIN£    CALCLILATdJ     THc    PUaER    SPcCrRUM    JF     Trli    fJ-'MCriON 
AOTuCOrih^LATION     CudFFlCI^NT     lo    STO!?l;l)     I  .-i    H     USI.MU    A    "riAIMNiNG 


aiMtiNSIGM    H(Ml.Nl) 
DIMcNSIC!<4    G{i£.55) 


THERE     WILL     da    54     INCR^^MtNTS     IN    FRtQUEMCY. 


L=2 
M  =  54 

MA=M+1 


INITIALIZE    G. 


00    3    I-l.Nl 
3    G(2.1 )=0 


CALCULATE    THE    SPECTRUM. 


DO    2    K^l.MA 
DO     1     N=2iM 
AM=M 
AN=,N-  1 
AK=K-1 

a-J.14159*AN/AM 

A=R<Ll,N)>i'(COS(d>l-l)*LOS(B»AK) 
IF(StLl  .N>  .EQ.-'J9.)     A=0 
1    U(2.K)=G(2>K)4^A 


WHEN    K=l     THE    PERIOD     IS    INFINITE    — FORCE    A    FINITE    GVERFLOW. 

IF(K.EQ.l )     GO    TO    4 

S(1.<)  =  108../(AK) 

GO    TO    2 
4    G(l  «K.)  =  100C00. 
2    G(2.K  )  =  (G(2.K.  )^R(L1.  1))«2 

PLOT    AS    A    HISTOGRAM. 

1HIST=1 

CALL    PLGT(G.L.L>N1.IHIST) 

RETUR.N 

END 


ivHOSc 
tVlUDOW. 
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FOKTrtAN  IV  G  LcVtL   21 


DATE  =  76J14 


0001 


0002 
0003 


0004 
0005 
0006 
0007 


oooa 

0009 
0010 


0011 
0012 
0013 
0014 
0015 
0016 
0017 


OOld 
0019 
0020 
0021 
0022 
0023 
0024 
0025 
0026 
0027 


002d 
C024 

0030 

0C31 
0032 
0033 
0034 
0035 
0036 
0037 
0038 
0039 
0040 
0041 
0042 
0043 
0044 
0045 
0046 
0047 


SUukUUriNt     PLOT     (A,L.W,N,IHISr) 

OIMcNSIUM    AtM  ,fj)  .3CALi(  1  1)  .UUrPUT(  101> 
OATA     A^jF  ,(3LNi<,oA.</lH*,  IH     ,1H|/ 

FINO     THE    MAXIMUM    ANJ    MIMIMUM    VALUi^    3F     THd    DATA     I  ,>»     Thii    AHHAY 

10     AMAX=A(L,1) 
AMIi\=0 
DO     1     I-1,N 
IFlAi L. 1) .GT.AMAXJ     AMAX-ACL, U 

THE    VALUE    -99     15    GIVEN     TO    'VJu    DATA"    PUlMTi—    TMESE    Arte     IGNO.JEO 

1F<A(L,  I  J  .E0.-9.^)    GJ    ro     1 

IF(  A(L,  I)  .LT.AMIN)     AMIN::A(L,I) 

1  CONTINUE 

SELECT     THE     PRQP^A    SCALE    FO^    THE    PLOT 

ISCALE=l 

JSCALE=0 

DO    2    1=1 ,50 

J  =  -1»I*-1 

IF(  AMAX  .GT.  I.ANO.AMAX.LE.I  I  +  l)  )     ISCALc=H-l 

IF(  AMIN.LT*  J.ANl>.AMIN.Gc.(  J-l)>     JSCALE=J-1 

2  CONTiNUE 

-PT*?^u-''^r!ii^"r5^    Gr^EATEa    THc.    MAXIMUM    OR     THE    AoSOLUTE  OF    THE    Ml^^IMUM 
?hL     .^^Vr-*!^*^.^"^'*'-    ^"    ^"^    ^<^AL-    MAXIMUM    VALJE.  .  .  MINIMUM 

I  MEN    JcCILc     IF    TriERw    ARE    NEGATIVE    VAL  "3     I '"J    THi=    rti  ta  _        tc     c<i    •r.j.r 

PLOT     «IUL    GO    FROM    THE    NEGATIVE    O?    ThI    MAXIMUM    3CAL^  POINT    fo    TH^ 

PusinvE  VALUE—  qthchwioE  the  plot  will   st"r t  at  2er J 

KSCALE=-1*JSCALE 

IF( IBCALt.LT.KSCALE)     I SCALE=KSCALE 

IF( JSCALE.EQ.O)  AJ=0 

IF( JSCALE.NE.OJ  AJ=ISCALE 

00  3  J=  1,11 

1  =  J-1 

SCALc( JJ=I»ISCALE 
SCALE ( J  J =SCALc( J )/l 0 

3  IF(JS'CALE.NE.a)     oCALE(  J)-5CALE(J»*2-AJ 
a(RITE{6.900J     tSCALE(J),J=l,UJ 

PLOT  THE  POINT 

aMAX=ISCALE 
DO  4  K=1,N 

iGO*To'4*'^**^'^*°'*^'^^*'*'^**V"^"^***'^°*°'*°^*'**^*"^*»''-Ea'0"> 

1=(  ACL,K)/BMAX}«100-)-.5 

IF{  JSCALc.Nc.O.)  I  =  (tA(L,K)/oMAX)«-l)/2*100.+-.5 

00  5  J=l,101 

IF(lHISr.EQ.l>  GO  TO  8 

GO  TO  9 

8  IF( I .GE.( J-1) )  GO  TO  6 
GO  TO  7 

9  IF(  t.ECi.I  J-1)  )  GO  TO  6 
GO  TO  7 

6  OUTPUT! JJ=AST 
GO  TO  5 

7  OUTPUT! J)=8LNK 
5  wONTlNUE 

IF( JSCALE.NE. 0)  OUTPUT! 51 i=aAR 

MRITE(6.:>C1  )  AIL, K)  ,A<1,K)  .(OUTPUT!  J>  .J=l  ,101) 
4  CONTINUE 

WRI TEI 6.902 )! SCALEIJ).J=1,11J 

12  RETURN 
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0051  ^    V02    KuRMATC     •.l5x..i'.10C I  •  )  ./ .oX  .  1 1  (  5X  .K  o  .  I  )  ,/ ,  •  1  •  ) 

C 
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FOSrSA.M    IV    G    LEVEL       21 


0001 
0002 

0003 
0004 

0005 
0006 
0007 
0008 
0009 
0010 
0011 
0012 
0013 


V«:**«-*-*'!t4**  *******  JUPoUL     V 


*-*****.*v  =  **»«*^** 


THIS    PrtOCRAM    ACCcrTi    AS     InPJT    OATA     rH£    SOLAH    jflflU    L)c:.\:aITY,     JUi'ITCfi 
IG    ANJ    NON-IU    OdStRVAT  lUN    i.-Ju>    ACTIVITy    ZjNi    CuUNTS.     ALL     V£RSUS 
JUL  I  A.N    DATE. 

ir     THEN    PKOOuCEi    bU>^cSPU3'iO    iJPOCHi    FUR    aJLAfi    ftINO    i)tNSITY    ANO 

THE    lU    AND    NON-iU    ACTlVi'Y     1  .Nii  1  Cr.-:. .     »•■*     THE     t^yMilErt    UH     CHANNELS     IN 

THc    cPjCH    (      idlNi     )     MUbT    u£    AH    OUij     iJirifGeH.     *»* 

FINALLY.     IT     CROoi    Cur.'rt  cTu  AT£S    THc    V/INJ    AND    JUPI  Tc:k    ACTIVITY    SUPER- 
POSED   EPOCHS    ANO    PLJTj    THt    ^£3ULTS.        THc     oUPt-.-;PL;S£tJ    EP.JCn:i    ARE 
THEN    SJdJECTEO     TiJ    THE    (AimImING    SMDuTHIiJC    PRJCIZDJi^E    AND    CRC5S 
CORRELATED    A.NU    PLQTTE')    AGAir4. 


>»»»   MAIN    <<«<.« 

DIMENSION  rflNDJ 2,3CJ) . 
DIMENSION  BEPOCH'.  .'.ej)  , 


IZ0NE3(4»,  1DATA(5,100C),  AEPGCHC 3 . 31 ) 
CuRCOFt  J,31>.  IFRc.Cl(*),  lEPGCHtS.Jl) 


C0r4M0N  INDEX.  PEHtOO,     IblNSi 
COMMON  CORCOF.  X jPT .  IFREO. 

INDEX=1 
CALL  INPUT 
CALL  EPCCHA 
CALL  EPOCHB 
CALL  COREL 
CALL  OUTPUT 
lNDtX=lNDcX+l 
GO  TO  1 
END 


WIND. 
INS MO I 


I  ZONES, 
I  COUNT  I 


lOATA,  AEPOCH,  etPOCH 
lEPJCH 
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FOHT.?AN     IV     G    I.eviiL        2J  ir.PUT  jaT£    :.    ZbJii 


0001  SIBkG'JJINE     INPUT 


lo/37/3^ 


lEAO     INPUT    DATA PERIUO.     lalNS,     ..IfW.     ZO;,es    ANO     IDATA 


PEiKiOO.     IBiNS    ANO    rtXND    ARt     RtAD     IN    CUuY    CNC£.     ThAT     IS    WH-N     I,NDi:x=l 

C  PtfilOO     IS    THd    PiHiUO    uF     Th£    SUPiiRPQSE.J    cPOCH. 

C  laiNS     IS    THE    NUMBER    Cr     ^INS     (CHANr.dL3)     THE    iPoLH     IS     -.Q^TcD    OVeP. 

C  WINO     IS    THc     SQL^H    WINu    DhNblTY    VS.     DATtd     AT     JUPir-^R 

C  WIND    MAY    HAVE     JCO    CCNScUUc  Fi  Vd     bNT^^IfSr  -^"PT-R. 

C  lOATA     MAY    HAVE    1000    £NTta£S    FJR     £ACH    :,LOCK     OF    JUPITEK    DATA. 

C  FkISSenOeII^     '-^^"    -"^'^''^^     '^^^^'^    "'^    "'^^'^^•^'    '^'^TA     ON    UP    TO    FOUR 


c 

0002  CONTINUE 


C 

C  >»»  DATA  DECK  STACKING  ««< 

C       CARD  l:  PtRIOO.  IBINS 

C       CARD  2:  FIRST  CARU  OF  SOLAR  «IND  DATA 

r  ^A-v,     ,.     V^?!    '-'^^^    '"'^^^     ""^    ^    i3LANK     CARD 

C  CA«U     I.     lOPT.....    ciUALS     I     IF     JUPITER    DATA     IS     TO    Bt    PRINTED, 

^  eGUALS    o^ANK     IF     IT    IS    NCT    TO    sE    PRINTdU 

C  CARD    J-     ZaiMF^  "^-r-^^T-^    if.  ?''^^?     *^    ''"'^    ^'"''^    ^F    T""^    OVERALL    DATA. 

I  CO .  i  cS52^orr;s'.i^j;if  .=.5f.^;s=-i  r^  '""■  -  - 

C  CARD  L:  IOPT  equals  2  IF  THcRE  ARE  NO  .MORE  JUPITER  BLOCKS 

C  AFTER    ALL    DATA.../*,.. 

C 

c 

ooS!        ^  ^g-s::  ^s^t::F/!S^??-i^^i^f'x;i^;j?:  iiss^i-  ^=^^^-  — •  ■^^^^-ch 

0007  ~  IFIIMOex.GT.lJ    GO    -.0    4 

OOOS  REAOC5,900)     PcfilQ-,,     IdlNS 

C  HEAD    IN    SCLAR    >lli,'i)    CATES    AND    OcNSITlcS. 

00?o  °°     *     1=1.291.10 

0010  M=I-^9. 

loll  R|AD(5.901J     <(K.N0(J.L).J=1.2>.L^I.M> 

00}j  IFIWlNOd. !).£). O.i    GO    TO    2 
XSH                                    1    CONTINUE 

001*  GO    TQ    3 

C 

C  IF    MIND    IS    our    CF   DATA.     FILL    THE    REMAINING    SPACES    WITH    ZEROS. 

0015  2    DO    3    I=M.300 

0016  DO    3    J=1.2 

0017  WINDC J.I)=0 

0018  3    CONTINUE 
C 

r  i-.i°T^    =    1. IDATA    AlLL     NOT     b£    PRINTED 

C  RcAO     lUPT.      IZONfcS    A.MC     lUATA     {JUPITER    DATA). 

0019  4    REAOIS,S04i     IGPT 
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G  t.tVfcL   21 


0020 
0021 


0022 
002J 
0024 


0026 
0026 
0027 


00  30 


0031 
0032 
0033 
0034 
0035 
0036 

0037 
0033 
0039 
0040 
0041 

0042 

0043 
0044 
0045 
0046 
0047 

0046 


UiPor 


J-\T'. 


7h  i  1  't 


16/J//C2 


IF    uUr     OF    DATA,     P«a^»A.M    E^ECutx     -.     i  r       ■ 


CALL    cXIT, 


5 
830 


•iOO 
901 
902 
903 
904 


IF<I0PT.ecj-2)      CALL      -.XI  T 
RtAL.<6.9C2)     (  IZO.^doC'n.I.i,^, 

INiriALIZa     IFREvJ,     A.NO     lA. 

00    9     1=1 ,4 

lA-1 

i^tAO  JUPITiiR  OATA. 

r=o 

t=l*l 
lCOUNT= I 

«c.AO(5,903)  (  IDArA(  J,l>,j=i,5)  ,ipj^£ji 

IF    lOATA    fs  Full,   return. 

IF( I  .uQ.looc)     GU    TO     7 

IF     UUT    OF     INPUT    DATA.    FIL'.    ^£.v,AIN,NG    SPACES     IN    lO^TA    .  ITH    ZtP.US. 

iF(IOATA(lKll .EQ.Oi     GO     Tl    6 

iNu'SAlt^N^EN^I^^if  I.^.^i,?'^"""^'''    "^'-^     ^"^    P^^VIOUS    ^NT.V    GO    TO 
iCfJ^^*^*^*      Ii=HcO(IA)  =  Ur(cQl 

GO  la  5°**'*^'''**^"**^"*'^'**'^°*"^'''^'^'*^^*''   ^^  *"o  a 

IA=IA«-l 

IFKtOi I A)=£FRc01 
GO    Ta    5 

1C0UNT=I-1 
00     7    J=  I, 1000 
00    7    <=!,£) 
lOATAiK.  J}?G 

CONTINUE 

Rt  TUftW 

FaRMAT<F4.2,I2> 

FtJRMAT(  IO{R5.0,i  3.0)i 

F0RMAT<4( 12, 2X)) 

FORMAT  (2X,t  5.  J»  .'f(7X  ,  13  J  .2X  ,  IPl 

FORMAT     Ul) 

END  -     — 
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FOHifKAN  iV 
0001 


^.tVtL   21 


cPGCHA 


0002 
O0C3 


G004 
0005 


0006 
0007 
OCOd 
0009 


0010 
00  u 
0012 
0013 
0GI4 
CClS 
0016 


0017 
0018 
0019 
0020 
0021 
0022 
0023 


0024 
0025 
0026 
0027 


0028 
0029 
0030 


0031 
0032 
0033 
0034 


SUoftOUTINE    tPUCHA 


JATl 


76314 


16/37/0  2 


400 


401 


403 
402 


411 


412 
4  04 


410 


THIS     SUBHUuriNc     FAKcj     The     JUPlTiZ9     /'i-j-    rri..T     ,.,.-'     -. 

COUNTS  P.H  t'o^L:  ^.tu:-::,h  :r^^L''i^'ll:^.'^'l^^rr^l^!.^^ 

^^^r^^O^^    ^^-^^-11    {-3^?^^^^^I.{^-?/-    —    ?-4C.A^ 

^^^  ^S^^^A/T^^??-i^---,--:  l^^:  ii^it:.  ^^'^-^  — ^ 

THIrtD    fiOrfSi  ^*    *^''    ^'^'<^-S     INTO    THc    ScCuNO    AND 

DO    400    J=l,31 
IcPGCHt  1,  J)::j 
DO    400     I-i.S 
IfcPCCH{I,J)=0 

COMPUTE    THE    EPOCH    CAY    CCHANN^L)    EACH    i^EAL    DAY    CORRESPONDS    TO. 

DO    40  3     Ul.ICOUNT 
JULIAN=     liJATAtl.IJ 
D A  Y- J  UL  I  AN/PE  R 1 UO 
I OA  r=  DA  Y 
OAY=OAY-IDAY 
dAY=DAy*It(INS 
IDAY=DAY+1 

|p3criA?^^^;;L'NELy°    OBSERVATION    .ONES    FUR     10    ANO    NON-IO    FG«    EACH 

00    40  3    J=2,5 

Oo''404'i;Kt8iNi"''''"'''''°*^'"'''*^*«-''" 

A£POCH{  1  tli-IEPOCHd  ,I» 

DO    404     J=2,4.2 

L=2+lJ-2)/2 

<=J+1 

IF     1E,»0CH    CONTAINS    NO    ENTRY    FOR    A    GIVEN    DAY    SET     ITS    VALUE     TO     --99. 

iF( iePOCH{K»I ).£U.O)     GO    TO     411 
GO    TG    412 
ItPOCH(KiIJ=l 
lEPGCMIJ,  I)=-99. 

A=IEPOCH{  J,U 
B8=IEPUCH(K,I J 
AEPUCH{LtI)=A/a8 

being^^SpuIed!  lf''trTs'''^o'i'il%^'ctV  T^^^ctivity  indices  A«e 

ZONES    A«E    BEING    TALL^Eot  AEPQCHCC.  1  )=Ba.     THE    •JB^.ERV ATION 

FILL    ALL    UNUSEO    SPACES    IN    AEPOCH    WITH    ZEROES. 
JaiNS  =  lBINS*-l 

00  410  j=jaiNS,3i 

DO    410     1=1,3 
AEPQCHC I,JJ=0 
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C  AEPUCH    N'JiV    CONTAIiXNi    :     THd    £.PJCri    JAY,     lu    ACTI VI TY     I  .NOl  0:15 ,     N  JN-1 0 

C  ACTIVITY     iNOlCES. 

C 

0035  RETURN 

003O  END 
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FOHTHAN 
0001 


IV  G  LtVEU  ^l 


dr^actm 


JATb  =  76J14 


16/J7/02 


SUoKOUriNt  EPOCH; 


0002 
OOOJ 


0004 
0005 


0006 
0007 
COOS 
0009 
0010 
0011 
0012 
00 13 
0014 
0015 
0016 
0017 


0019 


0020 
0021 
0022 
0023 
0024 
0025 
0026 


0027 
0028 


0029 
0030 
0031 


0032 


THIS    bUoHOUTlNii    TAKc3    THL    SJl.AR    V<  INO    J':;r;SITIES    Ar<.J    GfcNr.rtATco    A 

SUP^'iPUiiiO    cPuCn     .VITH    'IjImd'     CHANtldL^     AND    A\i     d.^UCH    PhriiOO 

or-     'i'-lilQO*.     »■!=**-    :NOTc:      frit    IJUMO^H     Or     I^AY'i     IN     TH,^     dPOCH     (l.HA.>I.NLLS 

).     THAI      lb,     'IBINS',     MUiT     oti     AN    000    PJolTlVil     I:>i Ti^oc^  .     "PhrtUJ*. 

HJttt.Vt£.<     MAY     dc     A.\Y    PUci  T  I  Vc     N'JM  jfcH  . ':^<  »* «  _ 

THt    Di-NSITY     rUk     -ACH     u^Y     (  CHANNcL)      IS    UJMALlZuJ     JY     TH-      TOTAL 

NUMLacH    CF    ACTUAL.    UAY:>    (.ONTH  I  oUT  U40    T.J    cACH    LPOCH    'JAY     (CHANNEL)  . 

THIS     INFOKMATION     IS     SlUrtiC.     IN    TWJ    CYCLuS,     IN     THi    AR«AY     -JiPJCH'. 

0IMFN3I0N    •.VINO(2.JO.'>.      IZlJNtS{4>,     I  OATAt  3  .  1  000  )  .     AiPOCHI  3  .  31  ) 
OlMcNSICN    3t;P0CH(.3,  i-J  .     CUr^CCF  (  J,  3  l)  .     Il->(t0(4) 

COMWUN     INDHX,     Pc.-ilOi,     I61NS.     .V I ND .     IZON£S,     I  DATA.     AtPCCH,     b-IPJCH 
CUMMtJN    COkCGH  ,     lUPTt     Ir^hJf     INSM3.     ICOUNT 

BfcPuCH     IS    Fl.-iiT    USV  J     TO     STOKd    THt    cPl>Cri    OAY     (  CHANNdL    #).     Tri:£    SUM 
□p     XHc    FLUX     OFNSIIIdS    rC«     THAT    DAY    ANU    TMu    NUM-Jr.-*    UF     hi£AL    OAYo 
FOLDEO    ONTO     THAT    ,.riA.N.\EL.     THUS     IT     ISA     THKEE    f<GW    ArtRAY. 

INITIALIZE    ocPOCr.    BY    COADING    THE    EPOCH    OAY     INTJ    AJ:i    ONII    AMO    ZEROES 
INTlI    rtOwS    Two    A,!..     TrthEc. 

708    N=300 

JblkS=2*IbINS 

00    70  0    J=1.63 

BEPOCHl 1 , J1=0 

b>£P0CH(2.  J)=0 

cJEPOCHi  3.  J)  =0 

IF{ J.Lc  .IBINS >    GO    TO    704 

K=J-I3INS 

IF( J.OT.Jdl N3)     K=0 

GO    TO    700 
704    K=J 

700  b£POCH(  l,J)-iV. 

COMPUTE    THc    ^POCH    DAY     (CHANNEL)    EACH    REAL    OAY    CORRESPONDS    TO. 

701  DO    703    1=1, ». 

IF     •i*IND«     tUT    UF    DATA    GO    TO    703. 

IF(WIN0(1.I l.EQ.O)    GO    TO    703 

CALCULATE    THE    EPOCH    DAY    FOR    THE    FIRST    CYCLE (lOAY)     AND    FQS    THc 
SECOND    '^CLc     IJDAYI. 

JULIAN=  aINDC 1 .1 ) 

OAY=JU'.iAN/P£RiDO 

IDAY=DAY  — 

DAV=DAY-IDAY 

DAY=DA\*IOlNS 

IOAY--0AY+1 

J0AY=.DAY+ISINS 

COUJI     THE    NUMdcft    OF    ^EAL    DAYS    FOLDbD    ON    EACH    EPOCH    DAY    OF    BOTH 
CYCLlS. 

fa£PLCH(  3.JDAY)=lJ=P0CH(3,  JDAYJ  +  1 
BT.P.,CH(3.  IDAYJ=t.cPOCH(3.  JDAY) 


702 
703 


SUk     THE    FLUX    IN    EACH    EPOCH    OAY    FOR    3CTri    CYCLES. 

'itPOCH  (  2  .  JO  AY  )  =bcPOCH  (  2  ,  JO  AY  )  ♦*  I  NO  <  2  .  1  J 
cj£30CH(  2,I0AY  )jid£P0CH(2,  iL>AY)+  »INO{  a.I> 
CONTINUE 

IF    NO    ENTRY     HAS    tiEEN    MADE    FOR    A    GIVEN    EPOCH    DAY    SET     BEP0CH(3.I) 
TO     I     TO    AVOID    DIVISION    BY    ZER*. . 

DO    709     l=l.jaiNS 

NORMALIZE  THE  FLUX  JY  THE  NUMBER  OF  REAL  DAYS  IN  EACH  CHANNEL  AND 
STORE  THESE  VALUES  IN  BEPDwHt 2 . I ) . 
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^''°^"^  OATe     =    7t,314  16/37/02 

C  8EP0CH<C...     .C.    CC.TM.S    ..P..P.SL.    ..r.S    A^D    SO...    .  U.O    P.UXS 


0036 
0036 


HETURN 
END 
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FOhTKAN 
0001 


0002 

000  i 

0004 
0006 

0006 
000/ 

oooa 

0009 


0010 
0011 
0012 
0013 
0014 
0015 


0016 
0017 
0018 
0019 


0020 
0021 
0022 
0023 
0024 
0025 
0026 
0027 


002a 


0029 
0030 
0031 


0032 
0033 
0034 
0035 


IV    G    LtVEL      21 


CCHEL 


SUDHOUTlNc    COHEL 


DATE    =    76314 


16/37/02 


:  THIS    CAUCUU.TES     THE    CO««ELAr£ON    CCEF     .Er.EEN    AEPUCH    AND     SEOPCH. 

[  ---l^n^^S^^;^-^;J----:3{---.-oo..    Ae.acH.3.3n 

300  JLIM=laiNS 
lUlMs  JL  IH 

301  DO      10*     1  =  2,3 

OO     >04     K-1,CLIM 

INIUALI^E    THE    VARIOUS    SUMMATIONS    TO    ZESO. 
SU..«  1  =  0 
SU  »i:sO 
SO t 1=0 

so; 2=0 

CRl  5;:0 

JDi-y=o 

CALCULATE    TH£    VAUIUUS    SUMMATIONS    FO^    A    SHIFT    CF    (K-,,     DAYS. 

00    303    J=1,JLI.M 

L  =  J+K-1 

ALPHA =A EPOCH (I , J > 

SETA=utPQCH(2,L) 

iDiv^Jof:i:[*'-°''-''^^^-'-T-='    GO    TO    303 
SUM1=SUM1+ALPHA 
SU.M2=iUM2  +  dETA 
SQ«  1  =  SCJR  1  +ALPHA**2 
SQft2=SGR2*dErA«*2 
CWOSstKtiS+ALPHAvUETA 
303    CONTiNUE 

AVOiO    DIVISION    bV    ZERO. 

IF< JOIV.EO.O)    GO    TO    3Ci 

CALCULATE    THE    VALUE    0^     THE    CORRELATION    COEFFIENT. 

A=CR0S-SUM1*SUM2/JDI\ 
D=SQ«.T(  SaRl-SOMl*»2/JDIV) 
C=SQR  T  (  SQK2-SUM<i  »:f2,  JUI V  ) 

AVOID    DIVISION    BY    >ERO. 

cCIf="(-S;c'*^*"*°'  ''°  ^°  ^02 

GO    TO    309 
302   CQEF5-99 

HAVING    NEGATIVE     SHIFTS     IN    THE    CrS^I    r  ^o^Jf^.f^''    ^'^'^    RcGARDED    AS 
WHICH    OCCUR    «IJH    HlGh-^R    juP^TtR    SIy    vj^ll^''^^  "i^X.'^°-'^^^^SE'- Y.    EVENTS 
AND    ARE    REGASOED    TO     HAVE    PoilTIVE    In^FVlf  "^PPEN    AFTE«    ENCOUNTER 

ARRANGE    THE    DATA     Iti    THE    ARRAY    Cn«ri-P    t-i     Tui*-r    ^ 
VALUH     IS    IN    THE    CENTER    OF    THE    ARRAY        -trwr^I    ^"^    ^-^°    ^^^Pf 
SHIFT    AT    THE    BEGlNNlNo    GUlNG    oCw.w    70    ztLn    T-J^^Mf^^iMUN    NcGATIVE 
MAXIMUM    PCSITIVTE    SniFT    AT    THE     ENoI  ^    ^^^    CENTER    AND    UP    TO 

309    A=lLlM 
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FOKIWAN     IV     &    LdVcL        21  COHEL  JATt    =    7ti31*  :fi/j7/C 

C  THE    y;*KlADLii     'iA*     is    TH:i    WALUc    OK     Trirl     l.-i-JcA    OK     Fdir.    ChNTLH    Or     1  Ht 

C  ARkAY    CURCOr. 


C 


I A=A/2.+.3 


C 

G  THc     WARIALBc     Mn"     HAS    THfc    VALUc    OK     IHc;     jHIFT    AOOjT     THE    CErJTc.a    OK 

C  THE    AK«AY    COHCCr. 
C 

0033  KA=-tK-l) 

0039  Io-IA  +  .<A 

0040  IK(K..Gr.IA)    GO    TO    305 

0041  GO    TO    306 

0042  305    <A=ILIM+KA 

0043  Id=ILIM  +  I8 

0044  306    CONTINUE 

004b  CJKCOF{ l,ia)=KA 

0046  304  COKCOFi  l.ia)  =  COcF 
C 

C  FILL  UNUbED  SPACES  IN  CGKCCF  WITH  ZEROES. 
C 

0047  ILIM=ILIH+1 

0048  DO    307    I=ILIM,3I 

0049  DO     307     J=l,3 

0050  307    CORCUFC J,I)=0 
C 

C  CORCOF    NOW    contains:     the     shift    A4NO     THE    COftHELATiON    COcIFF  ICIENT. 

0051  RETURN 
C 

0052  END 
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0001 


0002 
0003 
0004 

0005 
0006 

000? 


0008 
0009 
0010 
0011 
0012 


0013 
0014 
0015 


0016 
0017 
0018 
0019 
0020 
0021 
0022 
0023 
0024 
0025 
0026 
0027 
0028 
0029 
0030 


0031 
0032 
0033 
0034 
JO  35 
0036 
0037 


0039 
0040 


0041 
0042 


:VEL       21  PLOr 

SUBRUUTlNb     PLOT    (  A  .U  .M  ,N,  l-n  1  Sf  » 


OATh    =     763 1'» 


lb/37/02 


L)lMtNS;aN     *HNJ(c?,3C0)  ,      UONcSfi)  ,     1uATA(  5,1000  ,     AEt^UCHC  J.31) 
0IMi£N5ICN    liLPGCrilJ  ,o3)  ,     COr^CiJl^  {  3  .:>  1  )  ,     ti-Kdij(4) 
Dl.MiiNSION    A(M,N).     SCAL-llDi     JUrSvJT(lOl) 


IZO.Nii:i.     IUaTA.     AcPOCH,     dcPJCH 


COMMON     INDEX,    PcrJIOj,     ^iitNS,     WlNO, 
COMMON    COkCOF.     iuPT.     l^KcO.     INSMO 

DATA    AST.BLNK.BAk/IH»,1H     ,1H|/ 

FINU    THfc    MAXIMUM    AMU     MINIMUM    VALUt5    OF    THt    DATA     IN    THt    AiiUAY 

I  COUNTED 
10     AMAX=A(L,1) 
AMIN=0 

00  I     1=1, N 

IFIAIL,  IJ  .&r.AMA> )     AMAX-A(L.li 

THE    VALUE    -99    >K^5    GIVdN    TQ    "NO    DATA"    POINTS —    THdSE    ARE     iGNWtO 

IF(A(L,  I) ,EQ.-?9)    GO    TO    1 
IF(A(U,  I)  .LT.Ak  lit)     AMIN=A(L.  I) 

1  CONTINUE 

SELECT    THE    P'iUPHH    SCALE    FO.-^    THE    PLOT 

ISCALE=1 

JSCALE=0 

IF(AMAX.GT.i00.U«.AMIN.LT.-99)     GO    TO     13 

LA  =  1 

MA=100 

NA=l 

GO    TO     14 

13  LA=lu0 
MA-IOOO 
NA=100 

14  SO  2  1=LA,MA.NA 
J--l*I*» A 

Ii-(AMA<.GT.I.AND.AMAX.La.(  I+NA)  )  1SCALE=I+NA 
iF(AMI.\.LT.J.ANO.AMiN.Gt.(  J-NAJ  )  JSCALE^J-NA 

2  CONTI  Itc 

OtClDc.  WHICH  IS  GRcATER  THE  MAXI;^UM  OR  THE  AbSCLUTE  OF  THE  MINIMUK 
SET  THE  GkEATcR  ciGUAL  TO  The  SCALE  MAXIMUM  VALUE... 

THEN  DtCIUt  IF  THERi  ABE  NEGATIVE  VALUEi  IN  THE  OATA  —  IF  SO  THE 
PLOT  WILL  GU  F«UM  TmE  NEGATIVE  OF  THE  MAXIMUM  SCALE  POINT  TO  THE 
POSITIVE  VALUE —  OTHERWISE  THE  PLOT  WILL  START  AT  ZERO 

KS:aLE=-1*J SCALE 

IFaSCALE.LT.KSCAL6i     I  SCALE=KSCAkE 

I =(JSCALE.EQ.O)    AJ=0 

IK  JSCALE.Nc.O)    AJ=ISCALE 

to    3    J=l.ll 

1  =  J-1 

1 CALE C J ) =1» ISCALE 

SET    SCALE    VALUES    FOR    ALL    POSITIVE    DATA 

SCALE ( J  J  =SC ALEI J  J/1  0 

SET    SCALE    VALUES    FOR    POSITIVE    AND    NEGA.TIVE    DATA 

3  IF(JSCALE.NE.O>    SCALE! J1=SCALE (Ji* 2-AJ 
WRITE{b,9C0)     (SCALE(J)*J=1 ,11) 

PLOT    Trie    POINT 

BMaX=ISCALE 
DO    4    K=1,N 

IF    OUT    OF    DATA    GU    JU    4 
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FOHIHAN 
0043 


IV  G  LEVEL 
C 


PLOT 


UATE  -  7<:,j1'+ 


i  i/J7/02 


0047 
0048 
0049 
0050 
0051 
0052 
0053 
0054 
0055 
0056 


oosa 

0059 


0060 
0061 


0062 


0063 


0064 


0065 
0066 
0067 
0068 
0069 
00  70 
0071 


SCALE  THE  DATA  FRDM  0  TO  ICO  IF  ALL  t'QSITIVJ 

I=(  ACL. KJ/BMAX)*! 00+. 5 

SCALE     THE    DATA    FRCM    -60     TJ    +50     IF    POSITIVE    ANO    NcGATIVE 

IF{J5CALE.NE.O)     I  =  (<A{L,K)/l>MAX)+1)/2t10G.  +  .5 


BLANKS,     OR,     BAHS 


InD    S'^o:?-^'StJ^^^'     •*"-'-    ^^    FILLED    WITH    ASTESlXcS, 
AND    HjzP rtCLbt^UT S    UNd    LINE     OF    PRINT    FOrt     THE    PLOT. 

DO    5    J-l,10l 

IF  A  histogkam  display  is  desired,  'iHisr*  -  1. 

IF<  IHIST.EQ.l  )    GU    TO    a 
GO    TO    9 

8  IF(  I.GE.I  J-l)  )     GO    TO    6 
GO    TO     7 

9  IF(I.£Q.(J-1))     (io   TO    6 
GO    TO     7 

6  OUTPUT! J)=AST 
GO    TO    5 

7  OUTPUTC J)=BLNK 
5    CONTINUE 

iTLZ^'lEToAi^    l*,"Crr-'?    r5''-''r^,"?^S"'     ^"^    -^^^^tR    ELEMENT     JF    OUTPUT 

1F(JSCALE.NE.0>    OUTPUT  (  5n  =  aAR 

P.«INr     THE    LINE    UF    ASTERIXdS,    BARS    AND    BcANKS, 

4    CONTliud'*^*'     *"-''^' **''•''*  •''^^''^"^'■'*-'=^'^"'* 
PRINT     THE    CLOSING    SCALES    AND    HEADINGS. 

MiRITE{6,902)(SCALE(J),J=l,ll) 
ICauNT=lC&UNT+l 

IF    SMOOTHING     IS    TO    dE    I.NHIdITED.     •  INSMO*     =    1. 

IF( INSMO.EG.l )    GO    TO    11 

IF    THIS    IS    THE    FIRST    TIME    THIS    TAOLE    HAS    3EEN    PLOTTED.    THEN    SMOOTH 

IFClCOUNT.EQ.n     CALL    SMOOTH!  A  ,L.  M.N,  ICOUNTJ 

IF     THE     SMOOTHING    HAS    OCCUREO    GO    SACK    AND    PRINT    OUT    THE    PLOT. 

IF(  ICOUNT.EQ.n     GO    TO    10 

IF    RETURNING    AFTER     SMOOTHING,    ADVANCE    PAGE    ON    PRINTER. 


12  WRITE (6.903) 

11  RETORN 

90C  FuKMAT(lX,7X,ll(5X.F5.1>,/.l6X,«l«,10f ?•)! 

901  FORMAT(lX,F7.3,3X.F*.G.lX.ioiAii         ' 

902  FORMAT!  16X,'i», IOC I  •  )  .•,  8X.  1 1  (  5X.F5  II) 

903  FORMAT(»l«>  i  .>- .  ox,  l  H  bX.F5.  1  J  J 

END 
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CCOl 


oooa 


OOOi 
0004 
OOOb 
0006 
0007 
0003 
0009 
0010 
001  1 

ooia 

0Ci3 
0014 
0015 
0016 
0017 
OCIS 
0019 

Qo^o 

0021 
0022 
0Q^3 
0024 
0025 


0026 
0027 
0023 
0029 
0030 

0031 
00J2 


:veL  ii  cMOoiM 

SUBl^UUr  iNti    SMUOrn     <  A.LtMiN,  icjumd 


OATE    -    7b314 


THIS    SUBrtOUTI:Jc     IS    CALLcO    6Y     'PLOT'.     IF     Pd.-;FORMS-A    HANGING    SMOJTH- 
ING    UN     WrlAThVcK     AhicAY     'PLOT*     HAS    JUjT    Pi.urriiD     Q.M     T Hi.     P.iINri£R. 

DIMcHSION    A(.M.N) 

jaiNS=N 

00    10     1=1 .N 

10  lF(At  1  ,1).GT.  j.AND.A(l.(I  +  l)),i:Cl.O)     JiilNS-l 
Ab  =  A(L.  J31.MS) 

35=A(L.l) 

DU     11     I=l,JbINS 

J  =  I+1 

IF(  I.tQ.  Ji3I  NS)     J=:l 

C3=AtL. J) 

A0=2 

B0  =  2 

CD=2 

300=0 

IFfAS.Lr.C)     A0=1 

IF<CS.LT.O)     CD=1 

I^-(  (A0  +  CU).t:Q.2)     dDi)  =  -l 

D=(AU*CO) 

IF(bS.LT.OI     0=0/2 

IF(bS.LT.O)     60=1 

A(L|i  )=  (ciO-1  )*a3/(2ta00)+l  (AO-1  )*AS+(CD- n*CS  )/0 

IF(  (a04-r30D>  .tC.o)     A(L.I)=-1 

AS=BS 

11  dS=CS 

IF  raa  plots  this  size   wont  fit  on  one  pace,  advance  tj  next  page. 

IF(NiLT.41)    GO    TO    14 
mRITE(6.901) 
GO    TB     12 
14     WRITe(6,902) 

12  KETUfiLN 

901  FOftMATI  •  1»  K*  THREE    DAY    WEIGHTED    :>M00THING«) 

902  FOhHATC '0' ,••,•     THrtcE    DAY    WEIGHTED    SMJGTHiNG'i 
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rOVi  THAU 
000  1 


0002 
UOOJ 


IV  &  LEVEL   21 


DAi'c  -  7oJl4 


16/3//32 


0004 
0005 


0006 


000  7 

0008 


0009 

0010 
0011 
0012 
0013 
0014 
0015 
0016 
0017 
0013 
0019 


0C20 
0021 


0022 
0023 
0024 
0025 
0026 
0027 
0028 
U029 
0030 
0031 
0032 
0033 
0034 
0035 
0036 
0037 
0036 
0039 
0040 
0041 
0042 
0043 
0044 
0045 
0046 


SUBrt3UTiNE    OUTPUT 

OUTPUT    ROUTINES    FOLLUlm 

OIMF.Ni,lG;j    .Vl.NCX^.JCd),      IZUHESl'n  .     I  DATA  (  S  . :  000  )  ,     AEP0CHf-',31) 

COMMON    INDEX,     PERIOj,     loINS.     WIND,     IZONliS.     IDATA,     AEPOCH,     OEPOCH- 
COMMUN    CUkCOF,     lOPT.     IF^iLO.      INS.v.O,      ICJUNT.     lEPJCH    ''-'^"^"'     OEPOCH 

PRINT     ZONE    AND    F^EQUtNCY    riEAOINGS. 

WKire(6,605)     (I2CNcS(  I)  ,1-1  .-i)  ,     {IFkEC(J)  ,J  =  1,4) 

INITIALIZE     IPLOT    ANJ     IN5MO. 

IPLOT=0 
I NSMO=0 

IF     •40PT«=1,    PHINT    OUT     THE    JUPITER    DATA 

IFCIOPT.NE.U     GO    TO    303 

800  WRITe(6,603J 

DO    801     J=l,ICOUNr 

IFdDATAI  I,  JJ  .EQ.C)    00    TO    801 

.<RITE(-J,602JC  IDA  TAJ  I  ,  J)  ,1  =  1  ,5) 

801  CONUNUt 
mRITE(6,901) 
DO    806     1=1 ,31 
IFdiGT.lSI  MS)     GO    TO    809 

aOb    WRITE(o,602J     (iEPOCH(J.I),J=l,S) 
809    CONTINUE 

^M./^^o^Jib'-"!?"^'^    SEOUENCE    'PLOT*     IS    CALLED    TO    PLOT    AEPCCH,    B-iPOCH 
^Mnnr^rM^'^AK/!:":!^-,''''-'^^'     '^'^'"'-^    'SMOOTH*     WHICH    PtSFO,<MS    AN    MANNING 
?^    tI^.^    °r^^*'^'''^^"    '^'•^'^    BEPOCH.     THE    ORI-ilNAL    VALUES    ARE    REPLACED 

l^rl^'^   akhays  ay  the  new  values,  'plot*  then  plots  the  smoothed 

th^S^.*«^o,^^S°*^"A  '^I'^A'-LY  -COREL'  IS  CALLED  AGAIN  TO  CORRELATE 
DATA     A^POCH  AND  dEPOCH.  THEN  CORCOr  IS  PLOTTED  FOR  THE  SMOOTHED 

WRITE<6,905) 
GO  TO  802 

TM.r^=L:^"oiII^  °*^-*  '^  NOT  PRINTED  OUT  THEN  BEGIN  THE  PLOTS  ON 
THE  oAMc  PAGE» 

303  WRITE'(  6,904) 

802  L=2 
M=3 
N=31 
IHIST=1 

CALL  PLOT(AEPOCH,L,M,.N,  IHIST) 

WRITE(6,9C7J 

L=3 

CALL  PLGT{A£POCH,L,M,N,IHISrj 
805  WRITE (6,909) 

L=2 

M=3 

N=63 

IHIST=1 

CALL  PLOTJqEPOCH.L.M.N.IHISD 
807    KRITE(6.911) 

INSMU=1 

L=2 

M=3 

IHISr=0 

CALL    PLCT(C0RC0F,L,M,15INS, IHIST) 

WRIT£(6,913) 

L=3 

CALL    PLCT(CORCOF ,L,M.iaiN3. IHIST) 

I  PL OT=I  PLOT  i-1 

IPLOT  =  1  SAYS  THAT  THIS  IS  THE  FIRbT  TIME  THE  PHQGRAM  POINTER 
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2^4^  IFdMLOT.Nt.l)     GG    TO    a04 

00*^  CALL    COREL 

0"*9  WRIT<£{6,910J 

JO^o  GJ   ro   ao7 

"^OSl  tJ04    WETJHN 

0052  602  FOKMAT  (•   •,  is,  3<   4110I 

Ooll  90i  ^"SKMAfl'.o'.'f^ir^-'-'^^"  ^'^^  TOTALS., 

0002  ^11  p?l/,:;^:;-  ^^^.^^^'^i^VrfSw^^^SL^^V^/^^s^^-- — -, 


I. 
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